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Since in x-ray diffraction both experimental results and 
theoretical considerations are in agreement as to the 
presence of structure in liquids, an x-ray study is made of 
concentrated solutions of LiCl in the three normal alcohols, 
ethyl, propyl and butyl. A definite but tentative conclusion 
is that in these concentrated solutions of a strong elec- 
trolyte, both Lit and Cl- participate with the molecules of 


the solvent in forming a common liquid cybotactic or 
quasi-crystalline structure. This is a new conception which 
may be made to include what have been described as 
“onic clouds’ and as “association.’’ A conclusion of 
secondary importance is that the solutions are in fact 
solutions of the alcoholates of LiCl, although above the 
corresponding transition temperatures. 


HE problem of the structure of concentrated 

electrolytic solutions is both interesting 
and difficult. Great advances have been made, 
but so far removed is a satisfactory understand- 
ing that additional evidence secured by every 
method of approach possible is highly desirable. 
The study of the structure of liquids by the 
diffraction of x-rays has shown beyond reasonable 
doubt that liquids generally have an orderliness 
in the arrangement of their molecules. One 
would expect that the x-ray method of investi- 
gation, though crude, could shed valuable light 
for, after all, it is a very direct approach to the 
nature of the structure of a liquid. 

Since the x-ray results are only statistical 
time averages, the nature of the structure of a 
liquid cannot be obtained with the uniqueness 
possible in crystals. At the present time either 
one of two conceptions of the liquid structure 
seem to be plausible. One is that immediately 
about any point of observation, for a distance of 
from 10A to 20A, the arrangement of scattering 
centers approximates very closely to that of a 
crystal. If the liquid is near the freezing tempera- 


ture, the crystal structure referred to is that of 
the solid. This approximately regular structure 
may change materially with temperature. An 
important feature of the picture is that the point 
of observation specified is any scattering center. 
The liquid is therefore assumed to be homogene- 
ous. This conception could be only approximately 
correct, for this homogeneity would demand a 
repetition of like structure or a crystalline 
structure in fact, which obviously is not the case. 
It is in the approximation that allowance is 
made for temperature effects and the relative 
movement of the molecules. So the conception 
may be regarded as representing really a 
statistical mean rather than an instantaneous 
condition. The other plausible conception is not 
so clearly stated but is an instantaneous picture. 
It assumes that, at any instant, the liquid is 
spotted throughout with regions which are quite 
orderly, as assumed in the other view, fading 
away into regions with little orderliness. There is 
a constant shifting of position with the result that 
a given molecule does not continue to remain in a 
well-defined group. At any instant, then, part of 
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the molecules participate in one or the other of 
these regions. An orderly region may be called a 
“group,” though of course without sharp 
boundaries. The name ‘‘cybotactic group” has 
been given to it for convenience. The liquid is not 
homogeneous in detail. The latter conception is 
probably nearer the truth though the x-ray 
observations do not distinguish between the two 
views. But the former one has the advantage of 
being easily applied to computation. This has 
been done and reasonable structures have been 
found! to give approximately the observed 
diffraction curves. It is fair to say that in view of 
the large number of consistent experimental 
results? and of the satisfactory comparison with 
theoretical considerations, liquids generally do 
have a structure and that these results, which 
give a clearer conception of solutions, may be- 
come of quantitative importance. To this 
structure the name ‘“‘cybotaxis’”’ has been given 
because it is neither fluid nor crystalline but yet a 
crude and unstable imitation of the latter. At 
times reference has been and will be herein made 
to cybotactic ‘“‘groups,” by using thereby the 
second conception given above of the liquid 
structure. 

By interpreting the results now presented 
through the use of previous x-ray studies of 
liquids, the following tentative conclusion seems 
evident. In a concentrated solution of a strong 
electrolyte, LiCl, both the ions and the molecules 
of the solvent, n-alcohol, participate in the same 
liquid (cybotactic) structure and there is no 
significant change in the liquid structure as the 
unsaturated liquid passes through the solute’s 
temperature of transition to an alcoholate. The 
importance of the conclusion, if true, is evident, 
for it gives a new relationship of ions with ions 
and ions with solvent. 


1 For water, Bernal and Fowler, J. Chem. Phys. 1, 518 
(1933); for para-azoxyanisol, Buchwald, Ann. d. Physik 10, 
558 (1931); for n-paraffins, Warren, Phys. Rev. 44, 969 
(1933); for ethyl ether near the critical point, Warren, 
Phys. Rev. 44, 969 (1933). These computations follow the 
method of Zernike and Prins, Zeits. f. Physik 41, 184 (1927) 
and Debye and Menke, Phys. Zeits. 31, 797 (1930) and 
Ergebnisse der Technischen Réntgenkunde, Vol. II, p. 1, 
Akad. Verlagsgesell., Leipzig, 1931. 

2 For the lines of experimental evidence see Stewart, 
Indian J. Physics 7, Pb. 6, 603 (1933). 


PREVIOUS CONCLUSIONS 


Prins*® investigated water solutions of LiBr, 
Lil and KI. He applied the Debye-Hiickel ionic 
cloud theory to ascertain the probability distri- 
bution of the anions and found certain qualitative 
agreement but also quantitative disagreement. 
Later he* investigated water solutions of KCl, 
NaBr, SrBre, BaCle, Ag NO; and 
He concluded (1) that the ions are arranged in 
the solutions in a more orderly manner than 
would be expected according to the geometrical 
Spielraum, and (2) that the ionic cloud compu- 
tations do not agree with experiment. Prins’ also 
remarks that one may account for the absence of 
the typical gas diffraction picture by assuming 
that the ions have about them solid water shells 
which are independent of concentration and 
approximately monomolecular. 

Meyer’ has investigated the water solutions of 
NaCl and LiCl. The Cl- ions were sufficient in 
number to give a definite effect if arranged at 
corners of a cube 6 A on a side or scattered at 
random throughout the liquid. The diffraction 
curve gave no evidence either of this structure or 
of gaseous scattering from Cl-. 

Good’ has used aqueous solutions of LiCl, NH; 
SO,, CaCl, KCl, KB, KI, and finds in the case of 
the sulphate a new inner ring which he suggests 
might be formed by complexes of water molecules 
and NH,j*. 


EXPERIMENTAL METHOD 


Mo Ka radiation, Soller slits, methyl bromide 
ionization chamber and differential readings with 
SrCO3 and ZrO: screens* 0.036 g/cm? and 0.0025 
g/cm?,® respectively, were used. By this method 
the elimination of all but the Mo Ka radiation 
was effected. This is shown by the diffraction 
curve of an aluminum rod 3 mm in diameter for 


3 Prins, Zeits. f. Physik 56, 617 (1929). 

4 Prins, Zeits. f. Physik 71, 445 (1931). 

5 Prins, Naturwiss. 21, 435 (1931). 

6 Meyer, Ann. d. Physik 5, 701 (1930). 

7 Good, Physica Acta 4, 205 (1930). 

8 Ross, Phys. Rev. 28, 425 (1926); J.O.S.A. 16, 443 
(1928). 

® The number of g/cm? were chosen so that the trans- 
mission of wave-lengths other than the Mo Ka would be 
equal for the two screens and yet the difference in trans- 
mission for the Mo Ka would be an optimum. 
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STRUCTURE OF N-ALCOHOL SOLUTIONS OF LiCl 


which the readings with the filters were alike in 
all but the Ka region. The Soller slits were those 
previously described!® and gave a peak width at 
half height of approximately 26 seconds with 
crystal diffraction. 

The solutions were those of LiCl in n-ethyl, 
n-propyl and n-butyl alcohols. The percentage of 
saturation was obtained by dissolving the correct 
weight in the solvent. The saturated solution was 
secured by heating the base of a test tube 
containing the solvent and an excess of LiCl for 
two weeks, the temperatures of the base being 
maintained at 50°C. The solubility of LiCl in 
grams per 100 grams of the solvent was de- 
termined at 25°C and found to be 12 percent for 
n-butyl. This is in slight disagreement with other 
recorded values'! for the same temperatures 
which are respectively 42.3 at 50°C and 10.6 at 
25°C. The values of 24.7 percent for n-ethyl and 
16.0 for n-propyl, quoted in Seidell, were verified 
and adopted. The solubility of the salt in n-butyl 
was redetermined after the solution had been 
standing in contact with excess LiCl for more 
than seven months, and the result obtained was 
12.2 grams, an increase of less than two percent. 

The solutions were placed in the path of the 
x-ray beam in a very thin-walled cylindrical glass 
container, having a diameter about 0.8 of the 
“optimum” thickness,” although the optimum 
thickness precaution is not so important with the 
differential filter readings as when using a 
zirconium filter only. 

The temperature regulation of the solution in 
the path of the x-ray beam was obtained by the 
circulation of air in a heat insulated enclosure 
having very thin mica windows. For tempera- 
tures in the neighborhood of 0°C, the air was 
passed through a copper tube immersed in a 
bath of salt, ice, and water. For much lower 
temperatures, a current of cold air was obtained 
by evaporating liquid air from a special vessel 
containing an electric heater. The flow of the 
air and thus the temperature of the air bath was 


Stewart, Phys. Rev. 30, 232 (1927). 

"Seidell’s Solubilities of Inorganic and Organic Com- 
pounds, 2nd Ed., D. Van Nostrand, 1919. 

“This optimum thickness was computed from the 
following value of u. 


= (1.36 X N442.32N)/ZA 


wherein N is the atomic number and A the atomic weight. 
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controlled by the supply of heat through current 
regulation. 

The following are possible corrections of the 
readings: 


1. The effective difference in thickness of the fluid in the 
container at different angles of diffraction. 

2. The variation in the absorption of the container and 
mica windows and in their scattering at different angles 
of diffraction. 

3. The lack of parallelism of the x-rays and the penetration 
of the Soller slits. 

4, The incoherent radiation. 


For the small angles of diffraction used in this 
report, less than 15°, corrections 1 and 4 are not 
important, for as will be shown, no use is made of 
the exact position of diffraction peaks and great 
accuracy is not required. Corrections 2 and 3 
were approximated as follows. The readings were 
taken throughout the range with the container 
empty. These were reduced because of absorption 
of Mo Ka when the liquid is in place, and then 
used as corrections for the actual observations. 
In taking readings the current in the ionization 
chamber was always read first with one filter in 
position and then the other. Then the spectrome- 
ter setting was altered to an equal angle on the 
other side of 0°. 


RESULTS 


Figs. 1, 2 and 3 show the experimental curves 
for the x-ray ionization current with n-ethyl, 
n-propyl and n-butyl alcohols after correction is 
made for slit penetration, container and mica 
window scattering. Each point is the average of 
two, one on each side of 6=0°. The percentage of 
saturation is shown for each curve, all being 
taken at 24°C. Each curve has a different zero 
ordinate, as indicated on the left-hand margin, 
and all are altered to have the same maximum 
ordinate. 

Fig. 4 shows the chief lithium chloride alco- 
holate crystal diffraction peaks obtained with 
samples secured by plunging a saturated solution 
of each into liquid air, the purpose being a very 
rapid formation of the crystals and thereby a 
lack of preferred orientation. Measurements 
were made on two solidifications in each instance 
and the readings were in sufficient agreement to 
assert that the chief diffraction peaks of each 
alcoholate have approximately the relative values 
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Fic. 1. Relative diffraction intensities, Mo Ka radiation, n-ethyl alcohol solutions of lithium chloride. 
Fic. 2. Relative diffraction intensities, Mo Ka radiation, n-propyl! alcohol solutions of lithium chloride. 
Fic. 3. Relative diffraction intensities, Mo Ka radiation, n-butyl alcohol solutions of lithium chloride. 


indicated. The complete curve in each case was 
not obtained because, as will be seen, the dis- 
cussion involves only the chief intensity peaks. 
The temperature of observation for the alco- 
holates was kept at — 15°C which was well below 
the temperature for the solubility curve in each 
case. The temperature of the transition point or 
of the junction of the alcoholate and salt 
solubility curve is 17.4°C for the n-ethyl,"™ 
+3°C for n-propyl and —3°C for n-butyl as 
determined within 1° or 2° in these experiments. 

In Fig. 5 a diffraction curve for each alcohol is 
taken several degrees below its transition temper- 
ature. Again zero ordinates are shifted for 
convenience and the maxima are reduced to the 
same scale. The solutions were 75 percent 
saturated at 24°C, and then cooled. The solid 
alcoholates would appear with these solutions at 
10°C with the n-ethyl, at —9°C to —10°C with 
the n-propyl and at —20°C to 25°C with the 
n-butyl alcohol. The temperatures of observation 
are respectively 12°C, —7°C and —10°C and are 


13 Obtained by Turner and Bissett, J. Chem. Soc. 103, 
1904 (1913). 


in each case below the transition temperature 
with the samples wholly liquid. Under this 
condition of temperature and concentration the 
solution is in each case that of the alcoholate. 


DISCUSSION 


A reasonable deduction can be made directly 
by comparison of Fig. 5 with Figs. 1, 2 and 3, 
without the use of any assumptions. Each of the 
curves in Fig. 5 is closely like the curve for the 
same solution at 24°C as shown in the other 
figures. The ratio of major and minor peak 
intensities for a curve in Fig. 5 was computed and 
then compared with the like ratio in the corre- 
sponding curve in Figs. 1, 2 or 3. In the solutions 
with each alcohol the ratios are the same within 8 
percent and the peak positions are alike within 
measurement error. There is a small difference in 
a noticeable lessening of the minor peak width at 
the lower temperature, which, other experiments 
show, may readily be only a direct temperature 
effect. But since the solutions at the lower 
temperatures in Fig. 5 are obviously alcoholate 
solutions, then because of the likeness of diffrac- 
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S Relative Ionization Currents 
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Fic. 4. Most prominent Mo Ka diffraction peaks with 
powdered lithium chloride alcoholate crystals. 


tion, or the same arrangement of scattering 
centers, those in Figs. 1, 2 and 3, at the higher 
temperatures, would appear to be atso the same 


alcoholate solutions. 

A second conclusion can be obtained by refer- 
ence to previous experiments with miscible 
solutions. It is that the ionic LiCl alcoholate 
forms a liquid structure with the solvent mole- 
cules and that the diffraction curves arise from 
this structure or grouping. This is fundamentally 
important, if true. For it is a conception of the 
role of the solvent molecules, of the ionic cloud 
and of association, all in one. It has been 
conclusively shown" that when two thoroughly 
miscible liquids are mixed, one does not get two 
corresponding sets of x-ray diffraction peaks for 
the solution, but only one, this corresponding to a 
cybotactic structure in which both constituents 
participate. Indeed, the evidence is that the case 
in liquids is analogous to that in solid solutions.” 
Thus one would expect that varied mixtures of a 
hypothetical liquid alcoholate and a n-alcohol 


* Meyer, Phys. Rev. 38, 1083 (1931). 

* Of course the reference here is to ihe simple case of like 
crystals in the solid solutions. Experiments with liquids 
have not proceeded far enough to obtain analogies with all 
the cases found in solid solutions. 
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Fic. 5. Relative diffraction intensities, Mo Ka radia- 
tion, n-ethyl, -propyl, and -butyl solutions of lithium 
chloride at temperatures below the respective alcoholate 
transition points. 


would give shifting peaks with the pattern 
altering from that of the n-alcohol to that of the 
alcoholate as the amount of alcoholate in the 
solution is varied. The position of the diffraction 
peaks for the hypothetical liquid alcoholate can 
be judged from the diffraction with the crystal 
form. For it is well known in liquid diffraction 
experiments that the liquid curve has its broad 
diffuse bands approximately in the regions of the 
greatest ‘‘density’’ of peak energies on the crystal 
diffraction curve. With increasing concentration 
of the alcoholate in the n-alcohols, the peaks 
shift as if the liquid solute alone would have 
peaks at the small angle of 2° or 3° and in the 
region greater than 10°. But this is just what the 
crystal alcoholate diffraction curves would 
suggest on the ground of the peak density rule 
already stated. Thus the change in diffraction 
curves with concentration is clearly consistent 
with the tentative conclusion that the ionic 
alcoholates of LiCl forms a liquid (cybotactic) 
structure with the solvent molecules. It would 
seem that the conclusion from miscible liquids is 
justified whether or not the alcoholate may exist 
as a separate liquid. It is upon the effect of the 
final mixture and the crystal itself that the 
argument is based. 
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One. justification of the selection of the alco- 
holate as the solute in reaching the second 
conclusion has been expressed in the first 
conclusion above. Another is that the chief 
peaks for LiCl crystal for Mo Kea radiation 
would occur at 13.9° and 16° and liquid peaks 
in this region could not, according to the above 
argument, give the observed effect of the 
concentration of LiCl. The solute must be the 
alcoholate. It thus seems that, from our present 
knowledge of x-ray diffraction in liquids, one is 
led to the definite but of course tentative 
conclusion that the liquid groupings in the 
solutions are those of the alcoholates in normal 
alcohol. The important feature is, however, that 
there is a liquid structure in which both Lit and 
Cl- participate with the alcohol molecules. 

A confirmation of this view is that there is no 
evidence whatever for gaseous scattering due to 
Cl- near 0°.!7 If there were random ions of Cl- 
and if they were in sufficient number this gaseous 
scattering would be found. But in the case of the 
saturated solution of n-ethyl alcohol the number 
of chlorine ions is approximately one-fourth the 
number of alcohol molecules, and hence the 
number of scattering electrons in the former is 
about one-fifth of the latter. If the Cl~ ions give 
incoherent scattering as in a gas, there would be a 
noticeable relative increase in scattering near 0°. 
(The relative intensities in the various curves 


6 According to the measurement of Davey, Phys. Rev. 
21, 143 (1923). 


17 Attention has already been called by Prins (reference 


5) to this same absence of scattering in aqueous solutions. 
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have been reduced to the same value at the 
higher peaks.) In Figs. 1, 2 and 3, one finds no 
evidence whatever for gaseous scattering near the 
zero, for in Fig. 1 the curve for normal ethyl 
alcohol, if extended, would meet 6=0° at a point 
about 20 percent of the intensity of the chief 
maximum. For 20 percent saiuration this is 
slightly increased but for 80, 75 and 100 percent 
saturation the corresponding value is actually 
decreased. In n-propyl and n-butyl alcohol the 
differences are negligible. Although these curves 
have not been corrected and presented so the 
diffraction with the same number of alcohol 
molecules is shown, yet the effect of the Cl ions, 
if at random, would clearly be in evidence. 

Judging from these consistent evidences, one 
may then regard the conclusion as at least 
strongly indicated that the LiCl occurs in liquid 
cybotactic groups, in which both the alcoholate 
and alcohol participate. Presumably the LiCl is 
ionized and hence these temporary cybotactic 
structures may be regarded as analogous to ionic 
crystals. Indeed, it is already known that 
conduction in ionic crystals has definite simi- 
larities to conduction in electrolytes. Just how far 
this conception of the structure of concentrated 
electrolytic solutions will be found general one 
can only surmise. Experiments by the author in 
water solutions are at present in entire harmony 
with the view here expressed, but the details will 
be presented in a subsequent paper. 

I wish to acknowledge with appreciation my 
indebtedness for the observational work recorded 
herein to Dr. R. M. Zabel, research assistant. 
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An X-Ray Study of the Structure of Liquid Benzene, Cyclohexane and Their Mixtures* 


HERBERT K. WarD, Pennsylvania State College 
(Received January 15, 1934) 


Carefully purified benzene and cyclohexane exhibit 
x-ray diffraction peaks corresponding to preferred spacings 
of 4.68A and 5.09A, respectively, which are interpreted as 
the effective thicknesses of the respective rings. Mixtures 


of these two hydrocarbons give rise to two independent 
diffraction peaks indicating emulsion-type solutions, in 
which the disperse phase is too small to show even a 
Tyndall effect. 


INTRODUCTION 


TUDIES of the x-ray diffraction effects ex- 

hibited by liquids have given much informa- 
tion concerning the nature of the liquid state. 
These studies, at least in the case of liquids 
composed of molecules not having spherical 
symmetry, appear to substantiate Stewart’s 
cybotactic theory.' 

The spacings corresponding to the diffraction 
peaks exhibited by liquids composed of disk 
shaped molecules, such as benzene and cyclo- 
hexane, have been shown? to correspond to the 
effective thicknesses of the respective rings or 
disks. The positions of the diffraction peaks given 
by these two liquids have been measured also 
by other workers* with fairly concordant results. 


*From a thesis submitted to the Graduate School of the 
Pennsylvania State College in partial satisfaction of the 
requirements for the degree of Ph.D. in Physical Chemistry 
and Chemical Physics. 

'G. W. Stewart, Chem. Rev. 6, 483 (1929); Rev. Mod. 
Phys. 2, 116 (1930). 

For other reviews see: C. Brucker, Phys. Zeits. 29, 373 
(1928); W. Good, Helv. Phys. Acta 3, 205 (1930). 

*G. W. Stewart, Phys. Rev. 33, 889 (1929). 

-  §&P,. Debye and P. Scherrer, Nachr. Kgl. Ges. Wiss. 

Géttingen (1916); p. 1. 

E. D. Eastman, J.A.C.S. 46, 917 (1924). 

C, W. Hewlett, Phys. Rev. 22, 688 (1922). 

J. R. Katz, Zeits. f. Physik 45, 97 (1928). 

W. H. Keesom and J. de Smedt, Proc. Amst. Acad. Sci. 
25, 118 (1922); 26, 112 (1923). 

P. Krishnamurti, Ind. J. Phys. 4, 104 (1928). 

, C. V. Raman and C. M. Sogani, Nature 119, 601, 
1927), 

S. S. Ramasubramanyam, Ind. J. Phys. 3, 127 (1928). 

W. C. Pierce, Phys. Rev. 38, 1413 (1931). 

C. M. Sogani, Ind. J. Phys. 1, 357 (1927); 2, 97 (1927). 

E, W. Skinner, Phys. Rev. 36, 1625 (1930). 


Several binary liquid systems have been 
studied.‘ The earlier work of Wyckoff® indi- 
cates a pure addition of the diffraction effects of 
the components. Krishnamurti‘ finds some cases 
belonging to this class and postulates the 
coexistence of the original groupings, as: AAAA, 
BBBB. For other solutions he finds one halo or 
diffraction peak for the solution located at a 
continuously variable position between those 
corresponding to the components; for this 
(second) class, groupings such as ABABAB, 
A» B, Am Bn, AB, are postulated. Meyer,° 
however, finds all of the systems studied by him, 
excepting phenol-water below the critical solution 
temperature, belonging to Krishnamurti’s second 
class. Other studies of liquid solutions include, in 
addition to binary liquid systems, also cases of 
solid solutes in various solvents, chiefly water.’ 

Krishnamurti* suggests that the action of a 
solution toward the salting out effect (he here 
refers to aqueous solutions) might serve to 
indicate to which class the solution belongs. 
Meyer,’ however, perhaps more logically, places 


4P. Krishnamurti, Ind. J. Phys. 3, 331 (1929); 3, 507 
(1929). 

5 A. W. Meyer, Phys. Rev. 38, 1083 (1931). 

®R. W. G. Wyckoff, Am. J. Sci. 5, 455 (1923). 

7K. Banerjee, Ind. J. Phys. 3, 399 (1929). 

W. Good, Helv. Phys. Acta 3, 305 (1930). 

P. Krishnamurti, Ind. J. Phys. 2, 501 (1928); 3, 209 
(1928). 

W. C. Pierce, Phys. Rev. 38, 1413 (1931). 

J. A. Prins, Nature 123, 84 (1929); Phys. Zeits. 30, 525 
(1929). 

H. Shiba and T. Watanabe, Sci. Papers. Inst. Phys. 
Chem. Research (Tokyo) 10, 187 (1929). 

M. Ishino, S. Tanaka and A. Tsuji, Mem. Coll. Sci. 
Kyoto Imp. Univ. 13, 27 (1930). 
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solutions into the two classes: emulsion type and 
true solution type. In the former, two distinct 
types of cybotactic groups are postulated to 
exist, while in the latter both components enter 
into the same cybotactic groups, supposedly the 
groups have about the same composition as the 
macro-composition of the solution. 

There is no reason to doubt the technique, 
results, or interpretation of any of the workers 
as there is little duplication in the systems 
studied by the different workers. However, only 
qualitative agreement exists between the work 
of Meyer and of Krishnamurti on the two 
systems: phenol-water and tetranitromethane- 
cyclohexane. 


APPARATUS 

X-radiation 

A molybdenum anode Coolidge type x-ray tube 
(G. E. X-ray Corp.) was operated at 20 milli- 
amperes and at about 30 kilovolts r.m.s. The 
usual zirconium oxide in celluloid filter (supplied 
by G. E. Company) transmits selectively the Ka 
group of molybdenum (7.10A), however, there is 
an appreciable transmission of the white radi- 
ation at about 0.4A. This was corrected for by 
the use of a balanced strontium oxide filter* made 
by soaking several pieces of blotting paper, free 
from filler containing heavy atoms, in strontium 
nitrate solution, drying, and assembling with the 
use of Duco Household Cement. The whole filter 
was then impregnated with celluloid. 


Spectrometer 


The spectrometer was constructed with an 
horizontal axis of rotation in order that it may be 
used for certain studies planned for the future. 
Essential features of the spectrometer are given 
by Fig. 1. The scale given in the figure applies 
only in directions parallel to the x-ray beam. 
Dashed lines have been used to indicate an 
alternative position of the spectrometer arm. 
The mirror, lamp and scale used for the measure- 
ment of the angles of diffraction have not been 
indicated. The Soller slits? used for both colli- 
mator and ionization chamber each consist of 


8 P. A. Ross, Phys. Rev. 28, 425A (1926). 
W. C. Pierce, ibid. 38, 1409 (1931). 
® W. Soller, Phys. Rev. 24, 159 (1924). 
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Fic. 1. Diagram of spectrometer. (Scale applies only to 
directions parallel to the x-ray beam.) A, Mo anode of 
Coolidge type x-ray tube. S/it, Soller type slit (see text). 
Lig., Liquid being measured. F, ZrO, or SrO filter (see 
text). Jon.Ch., ionization chamber. 


ten parallel slits 19.7 cm long by 0.0795 cm wide 
by 1.27 cm, with lacquered iron foil 0.0076 cm 
thick as separators. 


Sample holder 


The sample holder was made of nickel plated 
brass with cellophane windows, cemented on 
with Duco Household Cement. The sample 
holder could be filled, emptied and washed 
without removing it from its position on the 
spectrometer. The spacing (1.3 cm) between the 
two cellophane windows is definitely less than the 
optimum thickness (about 1.6 cm) in order to 
avoid undue filtration effect on the x-ray beam by 
the sample. 


Ionization chamber 


The ionization chamber used (see Fig. 2) was 
blown from Pyrex with a Pyrex window about 
0.1 to 0.2 mm thick and was filled with methy! 
bromide at 74 cm pressure. Electrodes made from 


(0) 5 


Scale of cm 


Fic. 2. Ionization chamber. A, A’, Electrodes made of 
10 mil (0.25 mm) sheet nickel (polished). B, B’, 40 mil (1.0 
mm) stub nickel leads electrically welded to Alectrodes 
(A, A’) and to 40 mil tungsten leads used to make the 
metal-to-Pyrex seals. C, C’, Tungsten-to-Pyrex se@ls. DD’, 
Copper disk-to-Pyrex seals (see text). EL, Thith Pyrex 
window. 
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polished sheet nickel each 20 cm long by 2.6 cm 
wide and about 2.0 cm apart were used. A copper 
disk sealed’® into the Pyrex side tube provides 
electrical guarding for the collector electrode and 
lead. The side tube used for preliminary evacu- 
ation, testing for leaks, filling of the chamber is 
not shown; this was later sealed off. 


Vacuum tube electrometer 


The vacuum tube electrometer used to measure 
the ionization currents is shown schematically in 
Fig. 3. This amplifier uses the General Electric 


Std Cell 


Fic. 3. Vacuum tube electrometer. R,; 710" ohms. 
Rs, 2.5X 10° ohms. R;, 10 megohms. Ry, 1.0 megohm. R;, 
L. and N. student potentiometer. Rs, L. and N. dial 
decade resistance box 9(100-+10+1-+0.1) ohms. R;, Gal- 
vanometer shunt—10+100+1000+10,000 ohm fixed re- 
sistors in series with an inductance switch connected to all 
junctions. S; to Ss, Sealing wax insulated switches with 
mercury-copper contacts. S; to Soo, Switches of suitable 
mechanical construction—‘‘jack"’ knife, and Pyrex in- 
pee mercury copper contacts (these in the filament 
circuit). 


electrometer tube—Pliotron FP54." A number of 
circuits using this tube have been described 
recently,” being used for various purposes. The 
overall sensitivity of the amplifier here described 
is about 6X10-'* amperes per division using a 
rugged galvanometer* having a sensitivity of 


‘© W. G. Housekeeper, Trans. Am. Inst. Elect. Eng. 42, 
870 (1923). 

"G,. F. Metcalf and B. J. Thompson, Phys. Rev. 36, 
1489 (1930). 

A. W. Hull, Physics 2, 409 (1932). 

” A. DuBridge, Phys. Rev. 37, 392 (1931). 

Heil and Edwards, Phys. Rev. 43, 1 (1933). 

S. E. Hill, Science 73, 529 (1931). 

F. J. Moles, Genl. Elec. Rev. 36, 156 (1933). 

W. Soller, Rev. Sci. Inst. 3, 416 (1932). 

* Leeds and Northrup—Catalog number 2400a. 
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1.3X10-* amperes per division. This has been 
accomplished without recourse to rate of deflec- 
tion technique, usually used for the higher 
current sensitivities. 

The several sections of the circuit have been 
separately shielded as indicated by the dashed 
lines. The ionization chamber battery (B;) is also 
separately shielded, similarly the several storage 
batteries (Bz, to Bs) are shielded as a unit. 
Shielded flexible wires are used between the 
several shield boxes. All shielding and electrical 
guarding are connected together and jointly 
grounded. 

Sealing wax insulators made from a tested 
supply of sealing wax have been used throughout 
the high insulation part of the circuit. For 
measuring currents of various magnitudes any 
one of several grid leaks (Ri, Re, R3, Rs) may be 
selected by suitable switches (S: to Sj). Re- 
sistors supplied commerciallyt have been found 
satisfactory for use as grid leaks. Good quality 
commerical wire wound resistors have been used 
elsewhere in the circuit, excepting those used for 
protection purposes only. Potential for the 
ionization chamber is supplied by a number of 
large radio B batteries in series (B;), while other 
voltages required are supplied by storage 
batteries of large capacity (Bz to Bg). 

Leads X and Z, normally connected to ground, 
have been provided for test purposes and the 
intercomparison of the several grid leaks when 
using the FP54 as the indicating instrument.” 
The filament current of the FP54 may be 
measured with the potentiometer (R;) by the JR 
drop through a 10 ohm resistor in the filament 
circuit. When desired the plate (anode) current 
may be measured by transferring these leads to 
the 100 ohm resistor in the plate circuit as 
indicated by the dotted lines. 


MATERIALS 
Benzene 


Benzene was purified by successive treatments 
with sulfuric acid, treating with mercury, washing 
with water, drying, distillation, and several 
fractional recrystallizations. This method is 


tS. S. White Dental Supply Co., 154 W. 42nd St., 
New York, N. Y. 
13G, M. Rose, Rev. Sci. Inst. 2, 810 (1931). 
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similar to that used by Richards and coworkers.!4 
It had a boiling point of 80°C (uncorrected), a 
freezing point of 5.49°+0.02°C, and an index of 
refraction (°°) of 1.5008 as determined by an 
Abbe refractometer. 


Cyclohexane 


Cyclohexane was thought to contain benzene 
as one of its impurities, so it was first nitrated 
(cyclohexane is not affected) with an equal 
volume of a mixture consisting of two parts of 
sulfuric acid and one part of nitric acid to one 
part of water. Then the cyclohexane was given 
several successive treatments with sulfuric acid, 
washed with water, dilute sodium hydroxide 
solution, and several more water washes. At this 
stage a test portion did not discolor clean 
mercury upon prolonged shaking. The cyclo- 
hexane was then dried with calcium chloride, 
distilled, and further purified by several frac- 
tional recrystallizations. It had a boiling point 
of 80°C (uncorrected), a freezing point of 
6.45°+0.05°C, and 1.4262. 


Solutions 


Three solutions (or mixtures) of benzene and 
cyclohexane were made up by volume from the 
purified components. The refractive index of 
each was determined at 20°C by means of an 
Abbe refractometer. Table I gives the compo- 
sition and refractive index of each. Weight and 
mol percentages are calculated from the volume 
percentages and the densities of the components. 


TABLE I. 


Refractive 
Index 


1.5008 
1.4790 
1.4595 
1.4413 
1.4262 


Percent Cyclohexane 
Liquid Volume Weight Mol 
Benzene 
Solution No. 1 
Solution No. 2 
Solution No. 3 
Cyclohexane 


EXPERIMENTAL PROCEDURE 


The amplifier was warmed up by operation for 
about one hour before use in the measurement of 
ionization currents. During this period the 


“4 T, W. Richards, J.A.C.S. 36, 1825 (1914); 41, 2019 
(1919). 


amplifier was balanced by the adjustment of the 
variable 10,000 ohm and 400 ohm resistors in the 
plate circuit, after adjusting the filament current 
of the FP54 Pliotron to 0.110 ampere. The 
ionization current for the zero beam was checked 
periodically using low sensitivity at this time; 
however, no appreciable change was found over 
the period of this work. Then for the rest of the 
warm-up period, the collector electrode was tied 
to the control grid of the FP54 and the rest of the 
grid circuit connected for highest sensitivity, i.e., 
grid connected to ground through R,; and R; in 
series. After the warm-up period the initial 
current was balanced out with the use of R;; this 
initial current is due in part to the grid current of 
the FP54 and in part to the natural ionization in 
the ionization chamber in the absence of x-rays. 

Ionization currents were then measured by 
direct deflection of the galvanometer, waiting a 
period of one minute between readings. This 
wait is required since the relaxation time RC, 
(where R is in ohms and C is in farads the time 
unit will be seconds) of the grid circuit is longer 
than the galvanometer period and hence de- 
termines the speed of response of the system. 
This method of operation is speedier than the null 
method under these conditions (namely, with a 
grid leak of over 10'° ohms), and introduces no 
errors since the galvanometer deflection-grid 
voltage curve is linear over the range involved. 
The drift after the first hour of operation (warm- 
up) is negligible over a period of fifteen minutes 
or more under the worst conditions encountered 
during the work. 

At each setting of the spectrometer the 
intensity of the diffracted x-rays passed by each 
of the filters was measured by the vacuum tube 
electrometer and ionization chamber. The ioniza- 
tion current obtained using the strontium oxide 
filter was then subtracted from that obtained 
using the zirconium oxide filter giving that due to 
the molybdenum Ka radiation alone. This follows 
since they have the same transmission factor for 
the white radiation and the ZrO: filter selectively 
transmits this spectral group. Zero readings in 
each were obtained by inserting a lead screen in 
the path of the x-ray beam. At least three check- 
ing deflections were obtained at each setting of 
the spectrometer arm with each filter, and the 
three sets of readings were averaged. 


Intensity 


~ 


I 
spe 
pos 
anc 
bas 
sys 
the 
dec 
wit 
The 
age 
( 
I 
Or 
corre 
and 
is ne 
at la 


X-RAY STUDY OF THE STRUCTURE OF LIQUIDS 


In measuring the diffraction curves, the 
spectrometer arm was moved to the several 
positions. The angles were measured by a lamp 
and scale, and determined by calculations on the 
basis of the geometry of the lamp and scale 
system. 

Each diffraction curve was measured both in 
the direction of increasing and in the direction of 
decreasing angles. Both of the curves checked 
within the precision of measurement in all cases. 
The curves given herewith (Fig. 4) are the aver- 
age curves obtained from these pairs. 
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Fic. 4. Diffraction curves. 


Over the range included in this study the 
correction for scattering by the collimator slit 
and the cellophane windows of the sample holder 
is negligible as determined by experiment. At 
smaller angles the correction is not negligible and 
at larger angles has not been determined. 


RESULTS AND DISCUSSION 


Diffraction curves for benzene, cyclohexane, 
and their solutions are given in Fig. 4. All curves 
have been plotted to the same scale, hence 
comparisons may be made between the several 
curves. The angles of diffraction as given are one- 
half the respective angles of deviations and are 
thus equivalent to @ in the Bragg formula 


nd = 2d sin 0. (1) 


The x-ray diffraction curves for the solutions 
indicate the presence of each peak, although they 
are close enough together so that the resolution is 
not perfect. The curves for the 21.5 and 45.2 mol 
percent cyclohexane solutions show clearly both 
peaks in both of the original experimental curves 
for each solution. In the case of the 71.3 mol 
percent cyclohexane solution, the top of the peak 
is identical in position with that of the pure 
cyclohexane while the central portion is asym- 
metrical with respect to the top. This asymmetry 
is present in both of the original experimental 
curves while none of the curves for pure benzene 
or cyclohexane is asymmetrical in this region. 
This shows that the diffraction curves for the 
solutions arise from the combined effects of two 
preferred spacings, which are identical with 
those of the two pure components. 

Positions of the diffraction peaks for benzene, 
cyclohexane, and their mixtures are given in 
Table Il. The positions of the secondary peaks 


TABLE II. 


Main Peak 
Liquid d 


Benzene 
Solution No. 1 
Solution No. 2 
Solution No. 3 
Cyclohexane 


are, of course, not so well defined as those of the 
main peaks but are given as those of the com- 
ponents from which they arise, with which the 
curves are consistent. Table III gives the peak 
limits and widths (half intensity used arbitrarily). 
The values given in parentheses for the mixtures 
are the corresponding values for curves calculated 
by assuming that the diffraction effects of the 
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TABLE III! 


Peak width 
(at half 
intensity) 


1.20° 


Peak limits 
(at half intensity) 


3.73° 4,93° 
3.53 (3.50)* 5.04 (5.00) 
3.50 (3.50) 4.78 (4.70) 
3.48 (3.50) 4.63 (4.65) 
3.52 4.48 


Liquid 


Benzene 
Solution No. 1 
Solution No, 2 
Solution No. 3 
Cyclohexane 


0.95 


* Values in parentheses are from curves calculated by 
assuming strict addition of the diffraction effects of the 
components. 


1 Quotation from laboratory notebook. 


Synthetic 
1-3 


Collie 


0 

1.0 
15 
2.0 
25 
3.0 
3.3 
3.6 
3.8 
3.9 
4.0 
4.1 
4.2 
4.3 
4.4 
4.5 
4.6 
4.8 
5.0 
5.5 
6.0 
6.5 
7.0 
7.5 
8.0 


Column 1 headed @ represents various selected angles of 
the spectrometer. Column 2 headed C¢H¢ gives values of the 
intensity of the diffracted beam for the angles given in 
column 1. These intensities are read from the smooth curve 
drawn through the experimental points found for benzene. 
Column 3 gives corresponding values of the intensity of the 
diffracted beam from cyclohexane. Columns 4, 5 and 6 
give the calculated values for 3 mixtures. In column 4 the 
values are obtained by adding 50 percent of the values in 
column 2 to 50 percent of the value in column 3. Column 5 
is obtained by adding } of the values in column 2 to } of the 
value in column 3. Column 6 is obtained by adding } of the 
values of column 2 to } of the value of colume 3. The values 
in parentheses in Table III come from plots of the values in 
the last three columns which were then treated as if they 
were real data. 
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two components are strictly additive. These 
differences are more apparent than real as may 
be shown by plotting the experimental points on 
top of the synthetic curves. Trial shows no 
difference in any portion of any of the curves 
large enough to be significant. 

It will be noted that the peak for cyclohexane 
is distinctly higher and its width is distinctly less 
than in the case of benzene. There are over 
seventeen percent less molecules per unit volume 
of cyclohexane than of benzene, and since the 
same thickness of sample was used in each case, it 
means that there were less cyclohexane molecules 
in the diffracting volume. While the comparative 
diffracting power of the benzene and cyclohexane 
molecules acting as diffracting centers is not 
known, it is probable that the diffracting power 
of the cyclohexane molecule is not more than 
about fifteen percent greater than that of the 
benzene molecule. The sharper peak means that 
the semi-orderly arrangement in cyclohexane is 
more orderly than in the case of benzene and that 
probably a greater fraction of the cyclohexane 
molecules (more intense diffraction) are included 
in the cybotactic groupings. It was noticed that 
during the recrystallization the crystal nuclei of 
cyclohexane were much more persistent above 
the melting point than those of benzene. This is 
consistent with theory. 

The spacings corresponding to the diffraction 
peaks exhibited by benzene and cyclohexane 
have been shown’? to be the effective thicknesses 
of the respective rings. This is an intermolecular 
spacing as contrasted with an intramolecular 
spacing.* Hence it may be seen that the cybo- 
tactic groups of both benzene and cyclohexane 
coexist in their mixtures. Thus we may picture 
benzene-cyclohexane solutions as submicroscopic 
emulsions. The minimum size of the groups which 
are able to give rise to regular x-ray diffraction is 
smaller than the limit of resolution of an ultra- 
microscope even for the case where the difference 
of refractive indices is larger tian for the liquids 
used here. The extreme smal!ness of the cybo- 
tactic groups in the present experiments is shown 


16 W. P. Davey, Trans. Am. Soc. Steel Treating (1933). 

*If, as it seems unlikely, the peak corresponds to an 
intramolecular spacing, the following interpretation of the 
data taken on the solutions is not demanded, but it is still 
not inconsistent with the data. 
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by the fact that the ordinary rough tests showed 
no Tyndall effect in any of the solutions. It is 
believed that such emulsions are described here 
for the first time. 

The independent existence of the cybotactic 
groups of benzene and cyclohexane molecules in 
their solutions may be possibly connected with 
their deviations from the laws of perfect solutions 
in very dilute solution. Professor M. R. Fenske 
and Mr. M. R. Cannon" have found that the 
cryoscopically determined molecular weight of 
benzene in cyclohexane is nearly twice normal. 
It is here suggested that some of the cybotactic 
groups may have a regular enough structure to 
act as single molecules in their effects on the 
thermodynamic (vapor pressure connected) prop- 
erties of solutions, thus making it appear that the 
solute molecules are associated. True association 
should give rise to secondary peaks corre- 
sponding to distances in the associated molecule 
such as has been found, for example by 
Stewart! |? in the case of some polar compounds 
such as certain of the aliphatic alcohols. 


SUMMARY 


1. Carefully purified benzene and cyclohexane 


give x-ray diffraction peaks at 4.68 and 5.09A, 
respectively. 


‘6 Private communication, July, 1933. 
7G, W. Stewart, Phys. Rev. 35, 296 (1930). 
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2. A possible connection between the shape 
and intensity of the x-ray diffraction peaks given 
by liquids and the persistence of crystal nuclei 
above the melting point has been suggested. 

3. Solutions of benzene and cyclohexane have 
an emulsion type structure, in which the disperse 
phase is considerably too small to be seen with a 
microscope, and too small to show a perceptible 
Tyndall effect. 

4. A possible connection between the devia- 
tions from the laws of perfect solutions and the 
x-ray diffraction patterns of solutions has been 
suggested. 

5. A vacuum tube electrometer has been 
designed and tested, which has a direct deflection 
(or null method if desired) sensitivity of 6 10~™ 
ampere per division. 
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The lattice energies of NH,F and of the CsCl-type 
lattices of the other ammonium halides are computed by 
the method of Born and Mayer, assuming ionic crystals 
with a van der Waals potential. The calculations appear to 
be quantitatively satisfactory for the chloride, bromide and 
iodide. The method fails, however, for the asymmetrical 
fluoride lattice, and another method is employed. The 
theoretically calculated lattice energies in kilocalories are, 


NH,Cl, 161.6; NH,Br, 154.0; NH,lI, 145.5. From these 
lattice energies by mears of a thermal cycle the energy 
change at 0°K of the reaction NH3,,, +H *gas>NH4"*,.., 
which may be called the proton affinity of ammonia, is 
found to be — 206.4 kcal. With this proton affinity and the 
same cyclic process the lattice energy of NH,4F is calculated 
to be 191.4 kcal. 


INTRODUCTION 


HE ammonium halides, with the exception 

of the fluoride, are dimorphic, crystallizing 
in the CsCl-type lattice at low temperatures and 
in the rocksalt-type lattice at elevated tem- 
peratures. The fluoride possesses a Wurzite-type 
lattice. The interest in the anomalous behavior 
of the ammonium halides has created during 
recent years a considerable literature describing 
the physical properties of these salts. These data 
now make it possible to evaluate the lattice 
energies with some degree of accuracy. This cal- 
culation is carried out here following the method 
of Born and Mayer.' Reference should be made 
to their paper since the detailed equations and 
notations will not be given here. The lattice 
energies, U, of these uniunivalent ionic crystal- 
line salts, namely the energy changes at 0°K of 
the reactions are 
given by 


a’ C D 
U=9(ro) = + B(ro) + (1) 


ro 


The term varying inversely with the lattice 
distance, ’o, is the electrostatic energy. The terms 
varying as ro ° and are, respectively, the 
dipole-dipole and the quadrupole-dipole con- 
tributions to the van der Waals potential. B(r,) 
is the repulsive potential and hv,/2 is the zero- 


point energy. 


1M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 


CALCULATION OF THE VAN DER WAALS 
POTENTIAL 


The evaluation of the constants C and D from 
the polarizabilities, a, and from certain energies, 
e, of the ions in the lattice has been discussed by 
Mayer.” ¢€ is an energy corresponding to a 
weighted average frequency of absorption of the 
ion in the crystal. The values of ¢ for the negative 
ions of the ammonium halides are taken to be 
the same as the values found by Mayer? for the 
negative ions of the rubidium halides. The values 
used for the matrix elements, Qoo”, of the negative 
ions are the average values found by Mayer’ for 
the negative ions of the alkali halides. The 
polarizabilities, a, of the negative ions are then 
calculated from the relation 


(3/2) (ae/e?). 


The energy, ¢, of the positive ammonium ion 
is given the value 90 percent of the ionization 
potential (the potential required to produce a 
doubly charged ion) of the gaseous ammonium 
ion. This ionization potential of NHy,* is es- 
timated to be 31.7 volts by proportion from the 
ionization potentials of its electronic isomers 
CH,, Nat and Ne, which are, respectively, in 
volts, 14.5,3 47.04 and 21.47.4 The polarizability, 


2 J. E. Mayer, J. Chem. Phys. 1, 270 (1933). 

3 T. R. Hogness and H. M. Kvalnes, Phys. Rev. 32, 942 
(1928). 

4R. F. Bacher and S. Goudsmit, Atomic Energy States. 
McGraw-Hill, New York, 1932. 
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LATTICE ENERGIES OF THE AMMONIUM HALIDES 


a, of the ammonium ion in the lattice is cal- 
culated from its molar refractivity given by 
Wulff and Cameron? as 4.30. The matrix element, 
Qo, of the ammonium ion is then computed by 
means of Eq. (2). 

The sums of the polarizabilities of the negative 
and positive ions, Ya(calc.), Table III, row 9, 
have been compared with the experimental 
values, La(yellow), row 10, calculated from the 
indices of refraction, mp”, for yellow light.® 
The experimental values of Ya(yellow) will be 
lower than those which would be calculated 
were the index of refraction extrapolated to 
\= x, because of the infrared vibrations of the 
ammonium ion. An estimated value of Za(A= ~), 
calculated from the dielectric constants’ meas- 
ured at radiofrequencies, is given in row 11. 

The electron numbers, p, of the ions, which are 
needed to compute D, are taken as the geometric 
mean of the number of outer electrons, assumed 
here to be 8 for all ions, and the value calculated 
from the energy, e, and from the polarizability, 
a, by the method of Herzfeld and Wolf.? The 
lattice constants S¢’, S¢’’, Ss’ and Ss’’, needed in 
the calculation of C and D, and the Madelung 
constant, a, needed in the calculation of the 
electrostatic energy, are given in Table I. These 


TABLE I. Lattice constants. 


Crystaltype Wurzite Crystaltype Wurzite CsCl 
a (Madelung Se” 0.253 2.148 

constant) 1.638 a 1.633 1.155 
Si 4.354 z M 4 8 
0.762 M’ 12 6 


S;' 4.104 


constants as given for the Wurzite-type lattice 
are in reality those for the very slightly different 
Zincblend-type lattice. This approximation in- 
troduces no error in the final results, however, 
since it will be shown later that the lattice energy 
of NH,F must be calculated by a different 
method. 


°P. Wulff and H. K. Cameron, Zeits. f. physik. Chemie 
B10, 347 (1930). 
* International Critical Tables, Vol. VI, p. 99. 
je K. F. Herzfeld and K. L. Wolf, Ann. d. Physik 78, 195 
25). 


THE REPULSIVE POTENTIAL AND THE 
LATTICE ENERGY 


To calculate the repulsive potential, B(7o), it 
is necessary to approximate the ionic radii. In 
Table II the lattice distances 79, the distance 


TABLE II. Rocksalt-type lattice distances. 


Salt 


RbCI 
RbBr 
RbI 


ro X 108 cm 


3.267 
3.427 
3.662 


between the centers of closest ions, of the NaCl- 
type lattices of the ammonium and rubidium 
halides are compared. These lattice distances, 
taken from the Strukturbericht,> agree closely, 
and hence the ionic radius of the ammonium ion 
is taken to be the same as that of the rubidium 
ion. The ionic radii given in Table III, rows 22 
and 23, are taken from a paper of Huggins and 
Mayer.? 

The various steps in the calculation of the 
repulsive potentials and of the lattice energies 
are given in Table III. The values in parenthesis 
are averages of the values for NH,Cl and 
NH,Br. The compressibility data are computed 
from the measurements of Bridgman.” The 
thermal coefficients of expansion are from the 
work of Simon and Bergman." Both data are 
for atmospheric pressure and 273°K. The lattice 
distances, ro, as determined from x-ray measure- 
ments’ and as calculated from densities,*® are 
given in Table III, rows 20 and 21, respectively. 
The x-ray values are used here to calculate the 
lattice energies, except in the case of NH,F 
where the density value is used since the x-ray 
value appears to be of doubtful accuracy. The 
values of ro given for NH,F represent weighted 
means since this crystal is so deformed that one 
of the four neighbors of any ion is closer than 
the other three. The constants a, JJ and J’, 
needed in the calculation of the repulsive poten- 


8’ P. P. Ewald and C. Hermann, Strukturbericht 1913- 
1928, Akad. Verlagsgesellschaft, Leipzig, 1931. 

®°M. L. Huggins and J. E. Mayer, J. Chem. Phys. 1, 643 
(1933). 

10 P, W. Bridgman, Phys. Rev. 38, 182 (1931). 

" F, Simon and R. Bergman, Zeits. f. physik. Chemie B8, 
255 (1930). 
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Salt roX 108 cm 

| NH.Cl 3.266 

NH,Br 3.45 

NHgl 3.622 

= 


WILLARD E. 


TABLE III. Calculation of the lattice energies. 


NH«Br NHal 
CsCl CsCl 


NH,Cl 
CsCl 


Salt 
Crystal type Wurzite 
Qoo? (neg. ion) X 10!6 cm? 1.44 
Qoo? (pos. ion) X 1016 cm? 5.08 
(neg. ion) X10!2ergs/molecule 18.5 
(pos. ion) X 10!2ergs/molecule 45.3 
(neg. ion) X 10-4 1.18 
(pos. ion) X 10% 
Ya X 10% (cale.) cm? 
Ya(yellow) X10™ (exp.) 
= ) X10" (exp.) cm* 
C++ X 10° ergs cmé 
X 10° ergs 
X 10° ergs 
d44X 1076 ergs cm$ 
X 1076 ergs cm$ 
d,—~X1076 ergs cm$ 
C X10 ergs cmé 
D X1076 ergs cm$ 
ro X108 (from x-rays) cm 
ro X 108 (from density) cm 
X108 cm 
r+ X108 cm 
8 X10!2 (compressibility) 
barye™ 
B~'(ap/aT) p degree“ 
X10" barye! 
V-\@V/aT) p degree! 
(3T/NB) (AV X10"2 
ergs/molecule 
ro*{d2@(ro) /dr?] X 
ergs/molecule 
ae?/ro X10!? ergs/molecule 
C/ro® X 10!2 ergs/molecule 
/ro’ X 10!? ergs/molecule 
X10!2ergs/molecule 
o X10!2 ergs/molecule 
o/r 
p cm 
ro/p calculated 
B(ro) X 10!? ergs/molecule 
X 10!2 ergs/molecule 
U X10" (calc.) ergs/molecule 
U (calc.) kcal. 
U (thermal cycle) kcal. 


wre 


145.5 


tial, are given in Table I. The values of c¢ and &, 
which give the dependence of the repulsion on 
the charge, are computed from a theoretical 
equation of Pauling’s'? assuming 8 electron shell 
positive and negative ions. These values are 1.25 
and 0.75, respectively. The small zero-point 
energy is easily estimated. 

The values of the lattice energies in ergs per 
molecule and in kcal. calculated by means of 
Eq. (1) are tabulated in Table III, rows 40 and 
41. The lattice energy thus found for NH,F is 
more than 16 kcal. too small. The error is due to 
the fact that Eq. (1) holds only for perfectly 


"L. Pauling, Zeits. f. Krist. 67, 377 (1928). 


BLEICK 


symmetrical crystals, whereas NH,F is deformed 
or asymmetrical. It is obvious for reasons of 
stability that the nature of the forces producing 
the asymmetry must be such as to increase the 
lattice energy above that of a symmetrical 
crystal of equal lattice distance. The correct 
lattice energy of NH,F, as given in Table III, 
row 42, is calculated in the next section by using 
a thermal cycle involving the average proton 
affinity of ammonia found from the lattice 
energies of the other ammonium halides. 


Tue Proton AFFINITY OF AMMONIA 


In Table IV are given the energies of the reac- 
tions used to calculate the energy change at 0°K 
of the reaction NH3,,,+Ht eas which 
may be called the proton affinity of ammonia. 
The lattice energies are those of the salts cal- 
culated in Table III, and are the energies of the 
indicated reaction at 0°K. The heats of formation 
of the crystals are taken from the International 
Critical Tables, as are also the heats of formation 
of ammonia and of the halogen gaseous atoms. 
The heat of dissociation of fluorine is from the 
work of von Wartenberg and Taylor.'® The 
energy changes of the crystals between 291°K 
and 0°K are calculated by integrating the heat 
capacities measured by Simon, von Simson and 
Ruhemann." The small energy changes of the 
gases for this temperature change are estimated. 
The electron affinities of the halogen atoms are 
taken from a paper by Mayer and Helmholz’ 
where they were calculated from the lattice 
energies of the alkali halides. The sum of the 
AH’s of all of the first eleven reactions gives the 
AH, which at 0°K is the proton affinity of am- 
monia. The value in parenthesis in the NH.F 
column is the average value of the proton affinity 
of ammonia found from the other three halides. 
From this average value of 206.4 kcal. the A// 
of the first reaction, or the lattice energy of 


NH,F, is found to be 191.4 kcal. 


13H. v. Wartenberg and J. Taylor, Goett. Naclr., 
Math.-phys. KI., p. 119 (1930). 

4 Simon, v. Simson and Ruhemann, Zeits. f. physik. 
Chemie A129, 339 (1927). 

16 J, E. Mayer and L. Helmholz, Zeits. f. Physik 75, 19 
(1932). 
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3.40 4.08 5.27 
1 
15.7 13.8 12.0 
3.29 «4.49 6.66 Ne. 
4.99 6.19 8.36 | 
11 9.09 108 — 
13 127 208 399 = 
14 98 121 161 sum 
16 266 510-1192 + 
17 177 240 363 ee 
18 1251 1597 2286 
19 1898 2681 4423 
20 3.342 3.514 3.785 
21 3.347. 3.508 —— oni 
22 1.503 1.624 1.810 
dist 
25 14 — 
26 — — 
27 1.42 1.61 = 
28 
1.16 1.44 (1.30) 
29 
99.3 $92 — of 
30 12.00 11.41 10.60 By 
31 0.90 O85 0.78 
32 0.12 O12 O11 and 
33 17.2 16.0 14.8 whi 
34 170 156 i 
35 9.87 oul 
36 0.347 0.371 (0.359) evel 
37 9.63 947 10.54 
38 1.75 1.65 1.36 ca 
39 0.12 0.10 0.08 less 
40 11.15 10.63 10.04 
41 161.6 1540 
42 
(192 
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TABLE IV. Calculation of the proton affinity of NH; and the NH,F lattice energy. AHI in kilocalories of the reactions.* 


Reaction 


NH,Br 


NHuXer. 291°K—N O°K 

Nog. +21 Toy, 
NH 2 Nog. +3 Hog. 


2 

NH3g. O° ye, 291°K 
Hg. 0°K—+Hy, 291°K 

Xe, 291°K 

Hy, 

+e7 

Sum = NHsg. +H*g. 


0°K 
0°K 


| | 
w& 
a 


* Subscripts, cr. indicate crystal; s.s., standard state; and g., gas. 


The values found here for the lattice energies 
and the proton affinity may be in error by as 
much as 3 kcal. due to uncertainties in the lattice 
distances and compressibilities. 


DISCUSSION OF THE RESULTS 


Two references concerning the proton affinity 
of ammonia have been found in the literature. 
By means of a similar cyclic process van ArkeL 
and de Boer" find a proton affinity of 209 kcal., 
which is in good agreement with the value 
found here. The agreement is fortuitous, how- 
ever, since they employ a lattice energy of 148 
kcal. for NH,Cl. This lattice energy is 14 kcal. 
less than the one found in this paper and must 


A, E. van Arkel and J. H. de Boer, Physica 7, 354 
(1927). 


surely be in considerable error. By the method of 
molecular orbitals Mulliken’? obtains an am- 
monia proton affinity of 9.0 volts or 207 kcal., 
again in good agreement with the value found 
here. 

The good agreement obtained between the 
ammonia proton affinities deduced from the 
lattice energies of the ammonium halides is 
regarded as support for the accuracy claimed in 
the lattice energy calculations. It also shows that 
the ammonium halides are definitely ionic 
crystals to which the equations of Born and 
Mayer quantitatively apply. 

The author is greatly indebted to Professor 
Joseph E. Mayer for suggesting the foregoing 
problem, and is very grateful to him for the 


advice received in working it out. 


17 R,S. Mulliken, J. Chem. Phys. 1, 492 (1933). 


0°K 161.6 15 145.5 191.4 
~3.6 —5.0 —2.7 
291°K -75.1 —48.5 =111.7 
291°K 11.0 I 11.0 11.0 
Hoe. 291°K —51.0 -5 —51.0 —51.0 
291°K —28.8 —25.5 —31.8 
2.3 2.3 2.3 
1.4 14 1.4 
1.4 1.4 1.4 
0°K —312.0 —312.0 —312.0 
a 86.5 74.2 95.3 
— 206.3 — 206.2 (— 206.4) 
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The Theory of the Electrolytic Separation of Isotopes 


Ronatp P. BELL, Balliol College and Trinity College, Oxford 
(Received January 26, 1934) 


(1) Gurney’s treatment of overvoltage has been applied 
to the relative rates of discharge of two isotopes. As a 
first approximation it is found that the separation efficiency 
should be nearly independent of the electrode material and 
the current density. This agrees with the data at present 
available for the separation of the isotopes of hydrogen. 


(2) A more rigid treatment predicts small variations of the 
efficiency with the nature of the electrode and the current 
density, but is not at present able to predict these vari- 
ations quantitatively. (3) The application of the theory to 
other isotopes is indicated, and a rough estimate given of 
the possibility of separating Li and Li’ by electrolysis. 


INCE the original discovery of the electro- 

lytic separation of the isotopes of hydrogen 
by Washburn and Urey,' various suggestions 
have been made as to the mechanism involved. 
A difference in the behavior of the two isotopes 
may reasonably be expected at one or more of 
the following stages in the process of electrolysis 


(a) The mobility of the hydrogen and diplogen 
ions. 


(b) The rate of neutralization of the two ions by 
electrons. 


(c) The rate of combination of hydrogen and 
diplogen atoms on the electrode. 


In order to decide between these alternatives it 
is necessary to have some knowledge of the effect 
of various factors upon the separation, and it is 
only very recently that any work has been done 
in this direction. 

Possibility (a) cannot be an important factor 
in the separation, since it can be carried out 
equally well in acid or in alkaline solution.’ It is 
more difficult to decide between (b) and (c); 
however, it has recently been found by Bell 


1 E. W. Washburn and H. C. Urey, Proc. Nat. Acad. Sci. 
18, 496 (1932). 

2See C. H. Collie, Nature 132, 568 (1933); K. F. Bon- 
hoeffer and G. W. Brown, Zeits. f. physik. Chemie B23, 171 
(1933), 


and Wolfenden* and by Topley and Eyring‘ that 
the efficiency of separation is substantially the 
same for a number of different cathode materials, 
which makes it very improbable that the ob- 
served phenomena are due to the rate of com- 
bination of atoms on the electrode surface.’* It is 
much more likely that the difference in the rates 
of discharge of hydrogen and diplogen is due to a 
difference in the rate at which their hydrated 
ions are neutralized by electrons. Gurney® has 
previously shown that it is this rate of neutral- 
ization which is chiefly responsible for hydrogen 
overvoltage, and has developed a quantitative 
theory of overvoltage which agrees well with 
experiment, notably with the fact that the 
variation of current density with applied poten- 
tial is independent of the nature of the cathode. 
The author has previously pointed out*® that 
Gurney’s treatment can also be applied to the 
separation of the isotopes of hydrogen, when it 
predicts that the efficiency of separation will be 
independent of the cathode material, and gives 
a value of the right order of magnitude for the 
efficiency itself. The present paper deals with 
this point in more detail. 

For the current density due to the discharge of 
a hydrated cation, Gurney derives the expression 


3R. P. Bell and J. H. Wolfenden, Nature 133, 25 (1934). 

4B. Topley and H. Eyring, J. Am. Chem. Soc. 55, 5058 
(1933). 

5 As suggested, e.g., by A. Sherman and H. Eyring, J. 
Chem. Phys. 1, 345 (1933). 

6 R. W. Gurney, Proc. Roy. Soc. A134, 137 (1932). 
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ELECTROLYTIC SEPARATION OF 
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E,—«V 
exp 
E 


0 


where Ey is the energy of neutralization of the hydrated ion 
in its lowest energy state by an electron, 
EF, is the work function of the metal, 
V is the applied cathodic potential, 
x’ is the mean distance from the electrode at which 
the ions are discharged, 
e and m are the charge and mass of the electron, 
y is a correction factor a little greater than unity, 
introduced to account for the repulsive potential 
energy of the discharged ion and the water 
molecule immediately after discharge. 


If we apply this expression to the discharge of 
two isotopes X and Y under comparable con- 
ditions, the only factor which will not be the 
same in the two cases is Ep which will differ on 
account of the effect of the mass on the lowest 
vibration-rotation state of the hydrated ions. 
Gurney has shown that Eq. (1) can be brought 
into the form 


+log 
y-1 


(2) 


log i= 
ykT 


where ® is a complicated function having the 
property that its variations with V and 7 are 
negligible compared with the variations of the 
rest of the expression. Similar considerations 
show that the variations of leg 7 with E» are also 
determined essentially by the first term of Eq. 
(2), so that the relative rates of discharge of the 
two isotopes are given by 


log a=log (ix/iy) =[(Eo)x —(Eo)y J/ykT. (3) 


a may be termed the efficien¢y of separation of 
the isotopes. 

This equation predicts at once that the value 
of a should be independent of the nature of the 
cathode, and is thus in agreement with the 
experimental results for the separation of the 
hydrogen isotopes. It also agrees at least quali- 
tatively with the actual magnitude of a. Taking 
the mean value as a= 0.27 and putting y= 2 (the 
value found by Gurney to give agreement with 


7See G. N. Lewis and Macdonald, J. Chem. Phys. 
1, 341 (1933); Bell and Wolfenden, reference 3, Topley and 
Eyring, reference 4. 


ISOTOPES 


the observed values of d log 1/d V), Eq. (3) leads 
to (Eo)n—(Eo)p= 1900 cal. per mole. This is a 
reasonable figure for the difference between the 
zero-point energies of H and D in the ions H;O* 
and DH,0*, and is not very different from the 
value of 1400 cal. given by Eyring and Sherman® 
for the links O—H and O—D. It might be 
possible to obtain an accurate value for this 
difference from the Raman spectrum of the ion 
H;O0+, but the necessary data are not yet avail- 
able.* It may be pointed out that any theory in 
which the zero point energy is regarded as con- 
tributing directly to a heat of activation will 
lead to the expression 


log a=[(Eo)x — (Eo)y J/RT, 


(2mE)} 


E 1 
kT ¥ h 


which will only give the experimental value of 
a=0.2 by giving (Eo)——(Eo)p the improbably 
low value of 950 cal. per mole. 

The above treatment assumes that y has 
exactly the same value for the two isotopes, 
which would probably not be justifiable in a 
strict examination of the problem. The signif- 
icance of y appears from Fig. 1. The upper 
curve represents the potential energy of the ion 
H;0* plotted against the distance of one of the 
hydrogen atoms from the center of gravity of the 
rest of the molecule, while the lower curve repre- 


Distance 


Fic. 1. 


8 The fundamental frequency in the Raman spectrum of 
leads to — (Eo)p = 1600 cal. per mole. 


sents the potential energy of a discharged hy- 
drogen atom at different distances from the 
center of gravity of an H2O molecule. AB is the 
lowest vibrational level in the ion, and CD is the 
highest level to which an electron in the metal 
can make a transition. Gurney supposes that the 
curved portions EF and BD can be replaced by 
straight lines: if the ratio of the slopes of these 
straight lines is x, then y=1+ x. The potential 


(Eo)y—EiteV 
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energy curves in the figure will be unaltered for 
the systems DH,O+ and D+H.,0, but the lowest 
energy state will now be somewhat lower, say at 
A’B’. The curved portions to be replaced by 
straight lines are now E’F and B’D, and it is 
obvious from the general shape of the curves that 
the value of y will be slightly greater than in the 
first case. Eq. (3) should therefore be replaced 
in general by the equation 


yx (vy —1) 
(4) 


log a= 


yxkT 


where yx and yy are slightly different. Without 
a knowledge of the exact form of the potential 
energy curves it is impossible to predict the 
magnitude of the difference between yy and yp, 
but it may be pointed out that since Ey)— Ei +eV 
has a value of the order of 20,000 cal. per mole, 
a very slight change in y will produce a con- 
siderable difference in the value of a. Eq. (4) 
predicts that the separation factor will vary 
slightly with the work function E, (and hence 
with the nature of the cathode), and also with 
the potential V, so that it should be possible to 
decide by experiment how far the simple Eq. (3) 
is applicable. Meanwhile it is probably unwise to 
attach much importance to the exact value 
obtained for (Eo)——(Eo)p from Eq. (3). 

Both Eq. (3) and Eq. (4) give for the variation 
of a with temperature 


d log a/dT = —a/T. (5) 
For the separation of hydrogen and diplogen, 
Eq. (5) predicts a decrease in efficiency of about 
35 percent over the temperature range 0°—100°. 
There are at present no accurate data on this 
point, but this figure is not incompatible with 
the results of Bell and Wolfenden.* 

The above considerations will of course apply 
to the separation of any isotopes which can be 
discharged from a solvated ion. For a rough 
estimate of the efficiency of such a process for 
elements other than hydrogen, Eq. (3) will be 
sufficiently accurate and variations in y may be 
neglected, so that the factor which essentially 
governs the separation is (Ey)x—(£o)y, the 
difference in the zero-point energies of the two 


’ 
yykT vy (yx —1) 


isotopes in the hydrated ions. This quantity may 
be written 


(Eo)x (Eo) shvy = 
= (6) 


where f is the force constant of the link between 
the ion and the water molecule and Mx, My are 
the reduced masses of the two isotopes. It is 
immediately obvious that the separation will in 
no case be as efficient as it is for the hydrogen 
isotopes, since both the force constant and the 
value of the expression in {..... } will be smaller 
for all other elements. There are no data available 
for any reliable quantitative predictions, but the 
point may be illustrated by a very rough inves- 
tigation of the possibility of separating the 
isotopes Li® and Li? by electrolysis. The heat of 
hydration of the H+ ion is about 200,000 cal. 
per mole. The corresponding quantity (per water 
molecule) for the Lit ion is not known directly, 
but is probably about 15,000 cal.? Assuming that 
the force constant is roughly proportional to the 
heat of hydration, Eq. (6) gives for the relative 
separation efficiencies 


log 15, (- ~)/(i 
log ay “(x 000 
which leads to a,i;=0.96. Thus in this case 
separation by electrolysis will probably be less 


efficient than other methods, and conditions will 
be still more unfavorable for heavier elements. 


* Estimated from the heat of hydration of the solid 
lithium halides. 
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The Limiting Gas Density Method for Molecular Weights 


R. T. BrrGe ann F. A. Jenkins, University of California, Berkeley 
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In this, the first of a series of papers, we give a general 
discussion of the limiting gas density method for the 
determination of molecular weight. The various methods 
that have been used, or might well be used for the reduction 
of the data, are compared in a systematic way, and the 


most desirable methods for such work are pointed out. 
These methods are then applied to selected data given in 
the literature, and are compared with the actual published 
calculations and results. The paper closes with certain 
conclusions drawn from the preceding analysis. 


INTRODUCTION 


NE of us! recently made a photometric 
measurement of the relative intensity of 
the and Swan bands. This result, 
combined with Aston’s determination of the mass 
of C and with the known abundance of the 
oxygen isotopes, leads directly to an atomic 
weight of carbon of approximately 12.01, as 
contrasted with the international atomic weight 
12.00. 

In searching for the cause of this discrepancy, 
we have made a critical study of certain previous 
determinations of the atomic weight of carbon. 
Most of the recent work has been carried out by 
the “limiting gas density’? method. We have 
found, in this study, that a number of different 
methods have been employed in the reduction of 
the data, and in many cases the logical basis for 
the method is not clear to us. We have, accord- 
ingly, made a certain number of recalculations, 
and have in some cases obtained results appreci- 
ably different from those published. It has not 
been possible to make a complete investigation 
of all previous data. In many cases the published 
information is too incomplete to allow a recalcu- 
lation. It does seem desirable, however, to give 
at this time a general outline of the limiting gas 
density method, with illustrations of the various 
types of calculation that can be used and have 
been used in work in this field. From the calcu- 
lations we have already made, it would appear 
that the best value of the atomic weight of 


'F, A. Jenkins and L. S. Ornstein, Proc. Amsterdam 
Acad. 35, 1212 (1932). 


carbon is appreciably higher than 12.00 and that 
the discrepancy with the band spectrum value, 
12.01, may therefore not be serious. 

The results that we have obtained will be 
presented in a series of papers. In this, the first 
paper of the series, we discuss the general 
principles of the limiting gas density method. 
The second paper of the series will be devoted to 
the limiting gas density of oxygen, and the 
resulting value of R, the normal mole volume. 
In the third paper we shall discuss the atomic 
weight of carbon, as determined especially by the 
gas density method. 


THE LimitiInG GAs DENsity METHOD 


We define the density of a gas, at any pressure 
p, by 
(1) 


where w is the mass of gas used, and 7 is the 
measured volume. Now let 


Ly= pp/p=w/ po. (2) 


Then if p=1, L,=p,, the density at unit pres- 
sure. This is regularly denoted as merely L. 
Moreover, if we had a perfect gas (pv = constant), 
L, would be constant regardless of the value of 
p. Actually the value of L, is found, in the case 
of all real gases, to vary with p. Hence L, is the 
density the real gas would have, at unit pressure, 
if it acted like a perfect gas between the pressure 
at which the density p, is actually measured and 
unit pressure. This is a useful concept, for at 
p=0 any gas becomes asymptotically a perfect 
gas. Hence the limiting value Zi; of L, at p=0 
is the density at unit pressure that the actual 


Pp=w/r, 
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gas would have if it were a perfect gas. This 
quantity Ziim is of course a function of tempera- 
ture, but if all work is done at 0°C, then the 
product of Ziim and the normal mole volume R 
is, by definition, the molecular weight M, i.e., 


M=Lhiim-R. (3) 


As previously discussed by one of us,? R is the 
normal mole volume in liter-mole~! units, with 
245 = 980.616 cm-sec.~*, the units commonly em- 
ployed in gas density work. 

Eq. (3) is the fundamental equation for the 
determination of molecular weights by this 
method. The value of R is derived from the 
limiting density of oxygen, as will be discussed 
by one of us in the next article of this series.’ 
Then for any other gas one needs to know only 
Liim. This quantity, in turn, can be determined 
only by an extrapolation to p=0 of values of L, 
taken at various different pressures. The problem 
is then one of extrapolation, with all the uncer- 
tainty attending such a process. The uncertainty 
is of two kinds,—first that connected with the 
assumed functional form of L,, and second that 
connected with the purely experimental errors. 

Fortunately in the case of the permanent 
gases, the first cause of uncertainty is virtually 
absent. From work over a long range of pressures 
it is known that the curvature of L, when 
plotted against p, between p=0 and one atmos- 
phere, is so small as to be quite negligible. We 
are therefore entirely justified in assuming a 
linear relation and this has been done in all 
recent work on these gases. 

In the case of easily condensible gases, how- 
ever, the plot of L, against p sometimes shows a 
definite curvature. This makes the problem of 
extrapolation more difficult, and incidentally 
greatly increases the number of possible methods. 
Another addition to the list of possible methods 
comes from the fact that 1/L, is often plotted 
against p. This custom arose in the following 
way. We see from Eq. (2) that pv is inversely 
proportional to L,. Hence, except for the con- 
stant factor w, a plot of 1/L, is also a plot of pv. 
If w=1, the two become identical. We have 


2 R. T. Birge, Rev. Mod. Phys. 1, 1 (1929). See p. 12. 
(This article will be denoted G.C. 1929.) 
Thus R=(M)o,/(Liim)o, = 32/(Liim)o,. 
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ourselves occasionally been confused by the use 
in the literature of pv for 1/L,, and in the present 
article we shall therefore write pv only when p 
and v are alone known (i.e., w unknown). We 
shall always write 1/L, for the value of pv with 
w=1. 

Because of the reciprocal relation, Eq. (2), of 
pv and L,,, one gets the following useful relations, 
which also define \ and Xo! 


L/Liim = (p2) 0/ (pv) 1 =1+A=1/(1—Ao), (4) 


where and refer to zero pressure, and 
L and (pv); to unit pressure. In this article unit 
pressure is always one atmosphere. It would be 
much more logical to use LZ; and Ly in place of 
L and Lyin, but the nomenclature of Eq. (4) is 
so well established (except that P and V, or p 
and V, are sometimes used for our p and 2), 
that in order to facilitate comparison with the 
many articles quoted, we here adopt it. From 
Eq. (4) it follows immediately that 


ho =A/(1+A) and X=o/(1—Ao). (5) 


Returning to the possible methods for deter- 
mining Lj;,,, we then have two general classes. 
In the first class we plot L, against p and 
extrapolate to get the intercept Zii, at p=0. 
In the second class we plot 1/L, against p and 
extrapolate to get the intercept 1/Lii, at p=0. 
Each of these classes applies also to pv data. 
That is, if we have made observations only on 
p and v, then we may plot the products pv 


4 The quantity \ was introduced originally in connection 
with a linear function of p, and is particularly appropriate 
in this case. It is, however, now used for the more general 
case, and this has led to certain misconceptions that make 
it desirable to emphasize its correct definition, as given by 
Eq. (4). 

The particular symbol ) is in general use, except among 
the Germans, who write 1—a in place of 1+ (see Geiger 
und Scheel, Handb. d. Physik 2, 493 (1926)). In G.C. 
1929 the German nomenclature was adopted and this 
choice has since been criticized by Batuecas, J. chim. 
phys. 29, 26 (1932). The symbol Apo is being introduced for 
the first time, in place of the standard symbol Ao!, for 
reasons that will appear in connection with Eq. (57) 
ahead. Meanwhile we wish to emphasize the fact that in 
this article \ and Xo are defined always and only by 
Eq. (4), although the actual numerical values obtained for 
and Xo, from a given set of data, will naturally depend 
on how the data are treated. 
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against p, and extrapolate to get the intercept 
(pv)o. This is, except for the constant factor w, 
precisely the same as plotting 1/L, against . 
Similarly we may plot 1/pv against p, which is 
equivalent to plotting L,. 

If values of L, had been determined in all 
investigations, there would be no need to discuss 
pv values at all. Actually, however, much density 
work has been done at atmospheric pressure 
only. This yields the quantity Z. Then in order 
to get Liim it is necessary only to get values of 
pv at or near one atmosphere, and at various 
lower pressures. This is often done in a separate 
investigation. From these values of pv one 
obtains, by any one of a variety of methods, 
(pv)o and (pv);. Their ratio by Eq. (4) gives 
1+. Then, also by Eq. (4) 


Liim =L/(1+”) 
and by Eq. (3) 
M=LR/(i1+4). 


(4’) 


(6) 


This last equation is the standard form used in 
the literature, and often, as we have said, the 
value of L comes from one investigation, or set 
of investigations, and the value of (1+ A) from 
another set. We do wish, however, to emphasize 
the fact that in Eq. (6) L/(1+ A) is only a long 
way of writing the Zi; of Eq. (3) and that the 
problem of determining M, aside from the value 
of R, is solely that of determining Liim. If values 
of L, have been obtained, then Liim follows 
directly by extrapolation, and it only obscures 
the issue to introduce separate values of ZL and 
1+). However, in order to compare our results 
with those given in the literature, we shall 
indicate how L and 1+, as well as Liin, may 
be obtained from the data. 

The importance of the above facts is so great 
that we here summarize the four general ways 
in which the data may be handled. In all four 
cases what is desired is a f(p) and this is obtained 
most reliably by the method of least squares. 
By f(p) we mean, however, the adopted analytic 
representation of the data, regardless of how 
obtained. 

(1) One derives L, =f(p) from observed values 
of L,. The value of f(p) at p=0 is Liim, and this 
completes the problem. 
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(2) One derives 1/pv=f(p) from observed 
values of pu. The value of f(p) at p=0 is 1/(pv)o 
and at p=1 is 1/(pv);. The ratio gives, by Eq. 
(4), the value of 1+. This is all that can be 
derived from pv data. 

(3) One derives 1/L,=f(p). The value of f(p) 
at p=O0is 1/Ziim and this completes the problem. 

(4) One derives pu=f(p), and as in case (2) 
the values of f(p) at p=1 and p=0 give 1+, 
by Eq. (4). 

Since L, and 1/pv differ by only a constant 
factor, whatever is the best functional form of 
f(p) in case (1) is necessarily the best in case (2), 
for any given gas. The same thing is true of (3) 
and (4). On the other hand, the best function 
for (1) and (2) may be different from that for 
(3) and (4) and regardless of whether it is 
different or not, the numerical value of Liin 
derived from (1) will never, except by accident, 
agree with that derived from (3). If, however, 
one has set up a really satisfactory f(p) in both 
(1) and (3), the difference in the two resulting 
values of Liin will be found entirely negligible. 
Exactly similar remarks apply to (2) and (4). 

We now come to the much more involved and 
uncertain problem of the derivation of a suitable 
f(p), in any one of our four cases. With the 
exception of an arbitrary graphical extrapolation 
of the data, by means of some curved line, all 
methods are equivalent to the assumption of 
some particular functional form. Of all such 
possible forms, only two seem to have been used, 
viz., the first and second degree rational integral 
functions (r.i.f.) commonly called linear and 
parabolic solutions. We shall consider specific 
forms for each of these two functions, and noting 
again that cases (1) and (2) necessarily require 
the same function, and similarly (3) and (4), 
we have the following new four-fold division. 

(a) L, or 1/pv parabolic 
(b) L, or 1/pv linear 
(c) 1/L, or pv parabolic 
(d) 1/L, or pv linear 


Case (a) 


We begin with (a), the parabolic equation for 


L,=at+bp+cp* (7) 


where, by definition a=Z);,,. If now a value of 
a is obtained by least squares, or by some other 
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reliable method, this completes the problem. 
The following discussion, which is unfortunately 
somewhat involved, is necessary because L, data 
have quite regularly been used, in various differ- 
ent ways, to calculate also values of \ or Xo, and 
it seems desirable to consider the relative relia- 
bility of the values so obtained. 

To obtain \ or Xo it is necessary to evaluate 
a, b, and c, and as usual we recommend for this 
a least squares solution. Then \ and Xo are given 
by Eq. (4), in which both L and Lyin represent 
calculated values. Very often in the literature, 
however, the observed value of L has been used, 
and in order to distinguish the two resulting 
values of \ or Xo thus obtained, we shall for the 
remainder of this section write L, and L, for the 
adopted (observed) and calculated values of L, 
respectively.° Then, by definition 


L.=a+b+e (8) 


and as noted Liim =a. Hence by Eq. (4) 
=Liin/L-=a/(at+b+c) (9), 


or 


ho=(b+c)/(atb+c) and rA=(b+c)/a. (10) 


Eq. (10) gives the only correct values of \» and 
dX corresponding to Eq. (7). There is also a 
simple geometrical relation between A» and Eq. 
(7) which can be shown as follows. Let us define 


dy’ (1/L.)(dL,/dp) (b+2cp)/(a+b+c). (11) 


,’ can be thought of as the variable slope of the 
curve obtained by plotting L,/L. against p. The 
limiting value of this slope, at p=0, is 
(12) 
and this is not simply related to either \ or Xo, 
as is seen by comparing Eqs. (12) and (10). 


But the value of 4,’ at p=1/2 is exactly Xo, by 
Eqs. (11) and (10). This and other relations are 


5 In the general Eq. (4) any value of Z can be used, but 
it is obvious that the best value of A or Ao will be obtained 
only when one uses the best calculated value of L(=L-,). 
When, as is often the case, the value of Z is known more 
accurately than the values at lower pressures, a proper 
weighting of the various points will obviously lead to a 
calculated value of L(=L,) close to the observed value 
(=I). The same remarks apply to (pv);. 


R. T. BIRGE AND F. A. JENKINS 


—10 
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Fic. 1. Schematic diagram, applying to case (a), see Eqs. 
(11), (12), (24) and footnote 6. 


shown in Fig. 1.6 In place of the ,’ of Eq. (11) 
we might have defined another variable slope, 
obtained by replacing L. by Liim. In that case 
the value of this new variable, at p=1/2, is 
exactly . (See Fig. 2.) This interchange of 
and when and are interchanged, is 
very general. 


Fic. 2. Similar to Fig. 1, except that Zim replaces L. (or 
(pv)o replaces (pv)i-), and \ replaces do. In both figures the 
intercept values at p=0 and at p=1 atmosphere are 
indicated. 


Our chief purpose in writing Eq. (11) is to 
make a comparison with the following Eq. (13). 
which has been used in the literature. In Eq. (13 
the derivative of Eq. (11) is replaced by a 
special type of finite difference, and for con- 
venience L, is used in place of L,. We then have. 
by definition, 


Na’ = (13 


6 It is easy to show from Eq. (4) that Xo is the slope o! 
the secant of the L,/L.: p curve, drawn between p=! 
and 0, and the slope of the secant equals that of a parabol2 
at the middle point. 
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If now the values of \q’ are considered as a f(p) 
and are extrapolated by any method to p=0, 
one gets for this limiting value, remembering 
that p;=one atmosphere. 


(La—Liim)/La 
(14) 


This method then gives a value of Xo or of X, 
but it is not the best value, even if the best 
possible extrapolation is made, since it uses an 
arbitrarily adopted value of L, ordinarily the 
observed value, in place of the best calculated 
value. 

When a nonlinear equation is assumed in gas 
density work, it is quite customary to reduce it 
to a dimensionless form, by dividing the f(p) 
by some common factor. Such a method, how- 
ever, merely complicates the matter, so far as an 
evaluation 6f Lii,, is concerned, and leads to a 
value of X» or \ suffering from the same defect as 
the value given by Eq. (14). This may be shown 
as follows. The customary procedure is to divide 
each observed value of L, by the observed value 
of L. We then have, in place of Eq. (7) 


From Eggs. (7) and (15) we get 
L.=a/a" =b/b" =c/c”. 


(15) 


(16) 


The limiting value Liin/La at p=0, or a’, is 
commonly used to get \ or Xo as follows: 


Liim/La=a" =1—Xo=1/(1+d). (17) 


The desired value of Liim is evidently La’, and 
it is quite immaterial what numerical value of 
L, is used in Eq. (15), so long as the same value 
is used in getting Liim=L.a’’. Moreover, the 
value of Liim so obtained is identical with the a 
of Eq. (7). This is shown directly by Eq. (16). 
The division and subsequent multiplication by 
L.is thus entirely useless, so far as the calculation 
of Liim is concerned. On the other hand, the 
value of \y and \ given by Eq. (17) suffers from 
the same defect as does that given by Eq. (14), 
ie., the observed value of L has been assumed 
to be the best value. The correct expression for 
the best value of Xo or A resulting from Eq. (15) 
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is obtained, as usual, by first putting p=0 and 
then p=1, and taking the ratio. This gives 


Liim/L-=a" +b" +c") 
=1—do=1/(1+A). (18) 


This is identical with Eq. (9), because of the 
relations of Eq. (16). Jf now, by chance, L,=L, 
we have, by Eq. (15), a’+6’’+c’’=1, and Eq. 
(18) becomes identical also with Eq. (17).? 

Let us now consider the parabolic equation 
in 1/pv, corresponding to Eq. (7), viz. 


(19) 


To get the best value of \ or Ao we treat this 
equation just as we did Eq. (7), i.e., we evaluate 


7 There are a number of different possible values of \ 
and Xo, derivable from the same set of data, that should 
be distinguished. Taking Eq. (4) as our basis, we note 
that the value of \ or A» depends on, (1) the value adopted 
for Liim, and (2) the value adopted for LZ. In Eq. (7) the 
value of Liim(=a), which as just noted is identical with 
that obtained from Eq. (17), corresponds to a parabola 
passed freely through the experimental values L», prefer- 
ably by least squares. On the other hand, the value of Liim 
in Eq. (14), which comes in turn from Eq. (13), corresponds 
to some curve (it may be a parabola) that is required to 
pass through the adopted value L,. Hence the Liim of 
Eq. (7) or (17) is preferable to that of Eq. (14). The A or 
Xo derived from either Eq. (17) or (14) is, however, definitely 
inferior to that derived from Eq. (9) or (18), since in both 
Eqs. (17) and (14) we assume the adopted L,, whereas in 
Eqs. (9) and (18) we assume the calculated L, value given 
by Eqs. (7) and (15). 

We note also that the non-dimensional L,/Z, data may 
be treated by finite differences, just as in the case of Eq. 
(13). In fact, by merely putting L, inside the parentheses, 
Eq. (13) becomes the appropriate equation. We see that 
the A» of Eq. (14), expressed in this way, is merely the 
limiting value of the slope of a set of secants, drawn from 
the adopted L,/L,.=1 value at p=1 to the successive 
experimental L,/L, points. See Fig. 3. 

Certain further remarks about Eq. (13) may appropri- 
ately be made here. It follows immediately from the form 
of Eq. (13) that if one assumes for A,’ a Jinear equation 
in p, and gets a least squares solution, then the calculated 
values of L», or of L,/La, that may be derived by Eq. (13) 
from the calculated values of \,4’, satisfy a parabola. This 
calculated parabola, as already noted, necessarily passes 
through the observed L, point, but it is mot the least 
squares parabolic solution of the L, data under the 
condition that the parabola pass through the LZ, point. 
In general if Ay’ is assumed as a r.i.f. of the rth degree, 
this is equivalent to assuming L, to be a r.i.f. of the 
(r+1)th degree which is required to pass through a 
definite point. 


|| 

| 
qs. 
1) 
ise 

is 

Xo 

is 
. (or 
the 

are 
, to 
13). 
13 
a 
‘on- 
ive. 
(13 
pe of | 
rbola 
> 


172 R. T. BIRGE AND F. A. JENKINS 


a’, b’, and c’ by least squares or otherwise. We 
again differentiate between an adopted and a 
calculated value of (pv);, by writing (fv)1. and 
(pv) respectively. Then, by Eq. (4), 
1—do=1/(1 +A) = (p2)1-/ 

=a'/(a’'+b’+c’). (20) 


This corresponds to Eq. (9), and is obviously in 


- agreement with Eq. (9) since, from Eqs. (2), (7) 


and (19) 
L, (pv) =w=a/a' =b/b' =c/c’. (21) 


Corresponding to Eq. (10) we write, from Eq. 
(20) 
ho = (b’ +c’) /(a’ +b'+c’) 
and (22) 
(b’+c’)/a’. 


Corresponding to Eq. (11) we again get a simple 
geometrical relation by writing 


d(i/pv) 0b’+2c’p b+2cp 
dp a’+b’+c’ at+b+c 


The value of \,’ at p=1/2 is, by Eq. (22), 
exactly Xo. In this case \,’ is the variable slope 
of the (pv),./puv curve, considered as a function 
of p, but this is identical with the L,/L. curve 
mentioned in connection with Eq. (11) and 
shown in Fig. 1, since, from the inverse pro- 
portionality of L, and (pv), Eq. (2), we have 


L,/L-. (pv) pv. (24) 


Accordingly both sides of Eq. (24) are given as 
the ordinate in Fig. 1. We may now define a 
new variable slope, by replacing (pv):- by (v)o 
in Eq. (23). This slope, at p=1/2, will give \ 
(see Fig. 2). We may also treat 1/pv data by 
finite differences, but as noted near the end of 
footnote 7, we are thus effectively treating the 
data in a non-dimensional form, and we ac- 
cordingly take this up first. 
Dividing each (1/pv) by an adopted (1/(pv) 10) 


value we get 
(pv) 1a/(pv) =L,p/La=a" +b" pt+e''p. (25) 


In other words, due to the inverse proportion- 
ality of L, and pv, the non-dimensional Egs. (15) 
and (25) are identical. In the case of pu data we 
cannot of course evaluate Ljim, but the remarks 


(23) 


regarding the evaluation of \ or \» from Eq. (15), 
by means of Eq. (17), apply also to the present 
case. To be specific, 


(pv) 1a/(pv)o = Din /L 
= 1—d»=1/(1+A) (17’) 


gives a poor value of X» or A. If Eq. (25) is 
properly handled, the value of (fv)» corresponds 
to a parabola passed freely through the pv data, 
but the value of \» or \ in Eq. (17’) suffers from 
the use of the arbitrarily adopted (pv)1a value. 
The correct equation is 


(pv) Liim/Le=a"/(a"’ +b" +c"") 
=1—dy=1/(1+d). 


It is quite common practice to handle the non- 
dimensional (pv):./pv data by means of finite 
differences.* In analogy to Eq. (13), remembering 


8See A. Guye, J. chim. phys. 17, 141 (1919), for an 
extended discussion of this and other special methods for 
handling the data. For a criticism of some of his conclusions 
see G. P. Baxter, J. Am. Chem. Soc. 44, 595 (1922). It is 
well known that it is convenient to handle finite differences 
only when the abscissas are equally spaced, and the use 
of Eq. (26) in the literature is confined to this special case. 
We should like at this point to make some general remarks 
on the handling of finite differences. Let us suppose that 
a difference table shows that the successive second differ- 
ences are constant, within the limits of experimental error. 
Then one is justified in adopting a parabolic formula. 
Common practice is now to take the average second 
difference and apply it to the final first difference and the 
final observation to get the first extrapolated value, and 
so on for further extrapolated values. In such a solution 
only the first two and last two observations are used. No 
other observations have any effect on the result. The 
parabola passing through each of the extrapolated points 
has the following properties. It passes exactly through the 
last two observations and lies on one side of the first two 
observations by some small but equal amount. This 
particular solution cannot be obtained by least squares 
with any set of weights but it does effectively give infinite 
weight to the last two observations, equal but finite weight 
to the first two, and zero weight to all the others. A special 
case is that of three equally spaced observations. In such 
a case the parabola necessarily passes through each of the 
three points and the question of weighting does not enter. 
Also the calculated and observed values of pv coincide and 
there is necessarily only one possible value of Xo or >. 
Guye used the method only in this special case. The real 
question in such a situation is whether, with only three 
observations, there is any reason to assume a parabolic 
formula. The deviation from linearity may be due entirely 
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also the remark in footnote 7, we write 


Aa’ = [1 — 1a/pv p). (26) 


That this expression is identical with Eq. (13) is 
seen from the relation between pu and L, given 
in Eq. (25). Hence the limiting value again gives 
do. In this case Xo is the limiting value of a set 
of secants drawn on the (pv) i./(pv) curve from 
unity at p=1 to the successive experimental 
points (compare final paragraph of footnote 7, 
and see Fig. 3). 


0 ve 1 


Fic. 3. When L,/L, is plotted against p, the limiting 
value of the slope of a set of secants, drawn from unity at 
p=1 to the experimental points is one (but not the best) 
estimate of Xo (see footnote 7). This is also true if (pv)ia/pv 
is used as ordinate (see Eq. (26)). 


This completes the discussion of case a, which 
may be summarized as follows. One fits data on 
L, to a parabolic function, Eq. (7), and evaluates 
the coefficients a, b and c, preferably by least 
squares. Then the best value of the desired Lyin 
is given by the coefficient a, and the best value 
of \ or Xo, if desired, by Eq. (10). Xo is geometri- 
cally the slope, at p=1/2, of the L,/L. plot, 
where L, is the calculated value at p=1 given 
by Eq. (7). \ is similarly the slope, at p=1/2, 
of the L,/Liim plot. A non-dimensional form of 
the parabola, Eq. (15), leads by a longer process 
to the same value of Ziim as that given by Eq. 
(7), but as usually employed, to a less trust- 


to experimental error. An example of this situation is 
given in the last part of the section on Illustrations. 

In a similar way a set of first differences showing no 
trend justifies the use of a linear equation, but if one uses, 
for the extrapolation, the average first difference, this is 
equivalent to passing a straight line exactly through the 
first and last observations and ignoring all the others. 


worthy value of or A», Eq. (17), since an 
arbitrary value L, is used in place of the calcu- 
lated L.. The best value is given by Eq. (18). 
A commonly employed finite difference method 
for handling the data, Eq. (13), not only uses L, 
in calculating \ or Xo, but forces the adopted 
curve, which may be a parabola, to pass through 
the adopted L, point. 

In the case of 1/pv data one can calculate only 
dor Ao, the best values resulting from the adopted 
parabola, Eq. (19), being given by Eq. (20) or 
Eq. (22). The non-dimensional form, Eq. (25), 
is identical with Eq. (15), and leads to the 
correct value of \ or Ao given in Eq. (18’) or to 
the poorer value of Eq. (17’). The finite difference 
method for handling such non-dimensional data, 
Eq. (26), although often employed, suffers from 
the same defects as Eq. (13). Some of the 
relations in the case of the non-dimensional 
forms are given in Figs. (1) to (3). 


Case (b) 


The linear case (b) is merely a special form of 
case (a), in which the coefficient of the second 
degree term in p has become negligibly small. 
We therefore assume this coefficient to be zero, 
and thus get a linear equation. Every equation 
in case (a) can be used for case (b) by setting 
c=c’=c'’=0. Thus one takes Eq. (7) in the form 


L,=a+bp. (27) 


The best value of Liim is as before merely a, 
while the best value of Xo or \ is given by 


1—Ap=1/(1 +A) =Liim/L-=a/(a+b) (28) 
or 
ho=b/(a+b) and rA=b/a. (29) 


In place of a variable slope for the L,/L, plot, 
we now have a constant slope 


ho= (30) 


We can use Fig. 1 also for this case, if we assume 
that the secant shown in that figure is now the 
best representation of the data. 

In the case of the linear equation, the more 
general coefficients of Eq. (27) may conveniently 
be replaced by \» and L,. The result is 
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If in Eq. (30) we replace L, by Liim we get 
A= (1/Liim)(dL,/dp) =b/a, (32) 


i.e., \ is the constant slope of the L,/Liim plot, 
and the secant of Fig. 2 can be taken to represent 
such a plot. Introducing \ and Liim in Eq. (27) 
we get 


Lyp=Liim(i+ Ap). (33) 


There is no particular advantage in using, for 
case (b), the finite difference Eq. (13), although 
one might adopt for Xo the average value of the 
various secant slopes (A,’) thus obtained. Such 
a value of > has naturally the same defect 
discussed in connection with Eqs. (13) and (14). 

We may also use the data in the non-dimen- 


sional form 
L,/La=a"+b''p, (34) 


but as in the case of Eqs. (15) and (17), the 
resulting value Liim=a’’La is identical with that 
given by Eq. (27), while the use of a’’ as a value 
of 1—Xpo is inadvisable since it assumes the 
correctness of the point L,. The best value of Xo 
or ) is given by 


+b") (35) 


and Eqs. (35) and (28) are really identical, just 
as were Eqs. (18) and (9). 

The linear equation in 1/pv may now be given 
a similar treatment. We begin with Eq. (19) in 
the form 


1/pu=a'+b'p (36) 


and note that a trustworthy value of \» or A can 
be obtained only by evaluating @’ and b’ by 
some good method. The result is 


1—Xo=1/(1+A) = 
=a'/(a’+b’), (37) 
which is really identical with Eq. (28), on 


account of the general relations given in Eq. 


(21). Hence 
ho=0'/(a’' +b’) and dA=)'/a’. (38) 


By analogy with Eq. (30), we have in place of 
the variable slope of Eq. (23), the constant slope 


(pr) +b’) (39) 


and the secant of Fig. 1 can be used to represent 
this constant slope of the (pv):./pv plot. Intro- 
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ducing A» and (pv);,. into Eq. (36) we get 
(40) 


If in Eq. (39) we replace (pv):. by (pv)o, the 
result is 


d= (pv) 


as shown by the secant of Fig. 2. Introducing \ 
and (pv)o into Eq. (36) one has 


As already noted, it is rather common to use 
pv data in non-dimensional form, and in this 
linear case Eq. (25) becomes 


(pv) 1a/ pv L,/La =a’'’+b"p, 


which is identical with Eq. (34). The remarks in 
connection with Eqs. (34) and (35) therefore 
apply here. 

We note finally that Eqs. (30) and (39) are 
special definitions of \o, holding only for a linear 
equation, and similarly Eqs. (32) and (41) are 
special definitions of \. These special definitions 
appear implicitly in the linear Eqs. (31), (33), 
(40) and (42). 

Case (b) may be summarized as follows. One 
fits the data on L, to a linear function, Eq. (27), 
and evaluates the coefficients a and b. The best 
value of the desired Liim is given by a, and the 
best value of \ or Xo, if desired, by Eq. (29). 
Xo is geometrically the constant slope of the 
L,/L, plot, and \ that of the L,/Liim plot, as 
shown explicitly by Eqs. (30) and (32) and 
implicitly by Eqs. (31) and (33). A non-dimen- 
sional form of the data, Eq. (34), is of no practical 
advantage, and is often used to give a poor 
value of 1—X» (=a’) in place of the correct 
Eq. (35). In the case of 1/pv data, the linear 
Eq. (36) leads to values of Xo and X given 
properly by Eq. (37) or Eq. (38). Xo is geometri- 
cally the constant slope of the (pv),./pv plot, 
and dX that of the (pv)o/(pv) plot, as shown 
explicitly by Eqs. (39) and (41), and implicitly 
by Eqs. (40) and (42). The secants of Figs. 1 
and 2 also show these relations. 


Case (c) 

All of the necessary mathematical forms are 
included in cases (a) and (b). To get the analo- 
gous formulas for cases (c) and (d) one merely 


(41) 


(42) 


(43) 
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substitutes 1/L, for L, and pv for 1/pv. This 
change often results in the replacement of Ao by 
or of 1—X» by Aside from changes of 
this type the equations of cases (a) and (c) are 
strictly analogous. The analogy holds also for the 
assumptions made and for the reliability of the 
results. Hence it is unnecessary to repeat the 
greater part of the discussion of case (a), and in 
general we merely list each equation of case (c), 
with the analogous equation of case (a) indicated 
in square brackets. 
The basic parabolic equation is 


where, by definition, a=1/Zi;,, and 
1+A=1/(1—Yo) = (1/Liim)/(1/L-) 
=a/(a+B+y7) (45) [9] 


(44) [7] 


or 
A= —(8+7)/(at+8+7) 
(46) [10] 


and 


— (B+7)/a. 
For pv data the basic equation is 


po=al +B'p+y'p’. (47) [19] 


We define 
d(1/Ly)_ dm) 
dp (pric dp 
B+2yp +2y'p 
a’ 


(48) [11] [23] 


Hence i,” is the variable slope of the L./L, or 
of the pu/(pv);. plot, and by Eq. (46) equals 
—-\ at p=1/2. This and other relations are 
shown in Fig. 4, which is analogous to Fig. 1 of 


Fic. 4, Illustrative of case (c), Eq. (48), and strictly 
analogous to Fig. 1, of case (a). 


1 


Fic. 5. Illustrative of case (c), and related to Fig. 4 in 
exactly the same way that Fig. 2 is related to Fig. 1, so 
far as the ordinate is concerned. 


case (a). Replacing ZL, by Liim, or (pv)1- by (pr)o, 
in Eq. (48) gives a variable slope with a value 
— Xp at p=1/2. This is shown in Fig. 5, which 
is analogous to Fig. 2. 

In non-dimensional form we have 


L,/Ly= pr/ (pr) la 
=a" +B"p+y"p (49) [15] [25] 


and L,/a’’ gives a value of Liim identical with 
that obtained from Eq. (44). On the other hand, 
a’ is a poor value of 1+, or of 1/(1—Ab), 
(compare discussion of Eqs. (17) and (17’)). 


The best value is given by 


L./Liim = (pv) o/ (pr) le 
(50) [18] 


A common use of the non-dimensional data is in 
the finite difference form 


1 —(L./L») — pr/ (pv) 1a. 
pi-p 


ha” is the slope of the secant drawn from unity, 
at p=1, to any other experimental point, in 
terms of the L,/L, or pu/(pv);, data, and its 
limiting value at p=0 gives a poor value of —\, 
just as the limiting value of \,4’ of Eqs. (13) and 
(26) gives a poor value of Xo. This is illustrated 
in Fig. 6. 

Case (c) has been given in such condensed 
form that a summary seems unnecessary. 


Case (d) 
In a similar manner we write down, for case 
(d), equations analogous to those of case (8), 


(51) [13] [26] 


A 
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0 
Fic. 6. Illustrative of case (c), Eq. (51), and strictly analo- 
gous to Fig. 3 of case (a). 


starting with the assumed special form of Eq. 
(34), in which y=0. 


1/L,=a+Bp. (52) [27] 


The best value of 1/Zii is a, and the best value 
of X or Xo is given by 


=a/(a+B) (53) [28] 


A=—B/(a+B) and Ay=—B/a. (54) [29] 


For pu data we use the special form of Eq. (47), 


pv=a'+p'p, (55) [36] 
giving 

=a’/(a’ +8’). (53’) [37] 
r,”” of Eq. (48) now equals —X, the constant 
slope of the L./L, or pu/(pv):- plot, as given by 
the special definition 


d(i/Ly)_ 
dp (pv) le dp 
=B/(a+B) =B'/(a' +p’). (56) [30] [39] 


The secant of Fig. 4 can be taken to represent 
such a plot. If in Eq. (56) we replace L. by Liin, 
and (pv)i. by (pv)o we get 


d(1/Ly) 1 d(pv) 
dp (pv)o dp 
=B/a=p'/a’. (57) [32] [41] 


Hence —Xp is the constant slope of the Liin/L, 
or pv/(pv)» plot, which may be represented by 
the secant of Fig. 5. Eqs. (56) and (57) give 


—Ag= 


lim 


A. JENKINS 


the best values of \ and do, as obtained from the 
coefficients a’ and ~’ of Eq. (55). We note also 
that Eq. (57) is the origin of our symbol Xj, 
which bears the same relation to (pv)) in Eq. 
(57) that \ does to (pv);. in Eq. (56).° 

Putting L., Liim, and Xo into Eq. (52) we 
get the special equations 


(58) [31] 


1/L,=(1/Liim) (1 — op). (59) [33] 
Putting (pv);., (pv)o, and into Eq. (55) we 
get the special equations 


(60) [40] 
(61) [42] 


We may use the 1/L, or pv data in the non- 
dimensional form 


La/Ly= (pr) (62) [34] [43] 


As in Eq. (49), La/a” gives a value of Liin 
identical with that obtained from Eq. (52), and 
as in the case of Eq. (49), a’ is a poor value of 
1+ or 1/1—Xp. The best value is given, as in 
Eq. (50), by 


Le] = (pv) o/(P?) +B") 
=1+=1/(1—X). 


and 


and 


(63) [35] 


There is no particular advantage in using, for 
case (d), the finite difference Eq. (51), although 
one might adopt for —\ the average value of the 
various secant slopes (A4’’) thus obtained. This 
value of —\ corresponds to a linear graph forced 
to pass through an accepted point at p=1. 

This completes the discussion of the four cases, 
and we now turn to the literature for illustrations 
of the use of the more important of the preceding 
equations. 


ILLUSTRATIONS 


As our first illustration we consider the 
possible methods of handling a set of L, data, 
and for this purpose we choose five values of L, 


°On this same ground }, is a more logical symbol than 
\, but we do not wish to depart from such a well established 
nomenclature. 
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for (CH;)Cl, published by Batuecas.’° These 
are listed in Table I. 


TABLE I. 


p (at.) 1 3 2 3 i 
2.2895 2.2802 2.2713 2.2664 


L, (g. liter) 2.3075 


Each of these L, values is itself the average of 
many determinations, all of which appear in the 
original paper. Batuecas gives a plot of L, as a 
f(p). When replotted on a large scale, both L, 
and 1/L,, as f(p), show a small but perceptible 
deviation from linearity. This deviation may of 
course be due to some systematic experimental 
error, a question to which we revert later. For 
the moment we accept the data at their face 
value, and proceed to the analytic derivation of 
a value of Liim. In this work we use the method 
of least squares, giving equal weight to each 
point. In order to indicate certain relations, we 
have carried out the calculation for all four types 
of suggested solution, i.e., L, and 1/L, linear, 
and L, and 1/L, parabolic in p. In the following 
sections, unless otherwise stated, L always means 
the calculated value L.. 

(1) L, linear. Eq. (27) 


L, = (2.25291 +0.0001 2) + (0.05467 + 0.00019) p, 
s=1.33X10~ or one part in 17,300 of ZL." 


Hence 


Liim =a = 2.25291+0.00012, L=a+b=2.30758, 
\, from Eq. (4) or (29), =0.024266+0.000084. 


(2) 1/L, linear. Eq. (52) 


1/L,, = (0.443793 +0.000031) 
— (0.010454 +0.000051)p, 


s=0.353 or one part in 12,300 of L. 


T, Batuecas, Anales soc. espan. fis. quim. 24, 528 
(1926), 

" The standard deviation (S.D.) of a set of observations 
is s=(Zwe?/Sw)? where v=residual =Obs.-Calc., and w is 
the weight of each observation. When all observations are 
given equal weight (unity) as here, s becomes (2v?/n)}, 
where » is the number of observations. The S.D. differs 
from the average deviation of the observations from their 
calculated values only in that the root mean square 
average is used. The probable error of each coefficient 
evaluated and of the function itself is proportional to s. 


Hence 


Liim = 1/a= 2.25330+0.00016, 
L=1/(a+8) =2.30766, 
d, from Eq. (4) or (54), =0.02412+0.00012. 


(3) LZ, parabolic. Eq. (7)” 


L,, = (2.25249 -+0.00030) +0.056316p 
—0.0013042p?, 
s=1.087 X10~4 or one part in 21,230 of L. 


Hence 


Liim =a = 2.25249 +0.00030, 
L=a+b+c=2.30751, 
X, from Eq. (4) or (10), =0.024424. 


(4) 1/L, parabolic. Eq. (44) - 


1/L,, = (0.443948 +0.000056) —0.011068p 
+0.0004891 p2, 


s=0.201 X10~ or one part in 21,560 of L. 
Hence 


Liim = 1/a = 2.25252+0.00028, 
L=1/(a+6+y) =2.30750, 
X, from Eq. (4) or (46), =0.024408. 


The above data show the following significant 
relationships : 

(a) The proportional value of s for the linear 
and parabolic solutions of L, is not very different, 
and the two resulting values of Zii,, agree as well 
as is to be expected.“ On the other hand the 
proportional value of s for the linear solution of 
1/L, is distinctly larger, and the corresponding 
value of Liim differs from the parabolic solution 
by more than is to be expected from the probable 
errors. Hence we conclude that the linear solution 
of the 1/L, data is unsatisfactory. 

(b) Both parabolic solutions lead to almost 
the same value of Liim and to the same probable 
error. The latter results from the near equality 
of the two proportional values of s. We conclude 


In the case of parabolic solutions it has seemed 
unnecessary to calculate probable errors except for the 
important first coefficient giving Lim. 

'8 Their difference just equals the sum of their probable 
errors. 


177 
f 


178 


therefore that either parabolic solution is satis- 
factory, and slightly better than the linear 
solution of L,. 

(c) In spite of the smaller proportional value 
of s for the parabolic solutions, as compared 
with the linear, the corresponding probable 
errors are much larger. This is an illustration of 
the fact that as the number of undetermined 
coefficients in the adopted function is increased, 
the reliability of any extrapolation very rapidly 
decreases.'* We may elaborate this very im- 
portant point as follows. In the region of interpo- 
lation, the probable error of any adopted function 
is a correct indication of the reliability. But in 
the region of extrapolation we are assuming the 
correctness of the adopted function, and any 
calculated probable error is based on such an 
assumption. Hence if the assumption is incorrect, 
the probable error will be quite misleading. So 
in the present case, the relatively small probable 
error of Liim for the two linear solutions is a 
true criterion of the reliability of Ziim only if 
the function is known to be linear. In the case 
of (CH;)Cl, however, it is possible that the 
function is not linear, and therefore the probable 
error Of Liim, as calculated by either parabolic 


solution, is a much more conservative indication 
of the reliability. 

We accordingly adopt, as the best value of Liin 
obtainable from the data of Table I, the average 
of the two parabolic solutions, i.e., 


DLiim = 2.2525 +0.0003. 


It will be shown in the second article of this series 
that the best value of R is 22.4146+0.0003 
liter-mole™. Hence by Eq. (3) M=50.4889 
+0.0067 for (CH;)Cl. The international atomic 
weights for H (=1.0078) and Cl (=35.457) seem 
well established and should have negligibly small 
errors compared to that in M itself. Hence 
C= 12.0086+0.0067 from the data of Batuecas.'° 
In order to compare our results with the calcu- 
lations of Batuecas, we note that the best value 
of \, as obtained from the average of the two 


14 This fact has been well shown by H. Schultz, J. Am. 
Statistical Assoc. 25, 139 (1930) (see plots on p. 159). 
This important article was unfortunately overlooked in a 
later discussion of the same subject by one of us (R. T. 
Birge, Phys. Rev. 40, 207 (1932)). 


R. T. BIRGE AND F. A. JENKINS 


parabolic solutions is 0.02442. All four types of 
solution give a value of L close to the observed 
2.3075, and the parabolic solutions give this 
exactly, to five digits. 

Batuecas makes the calculation in a quite 
different manner. He first plots L, against p and 
makes a graphical linear extrapolation, obtaining 
Liim = 2.2527, compared to our 2.2529 for the 
same function. Batuecas’ article also includes a 
number of sets of measurements of pv at various 
pressures, the highest pressure being usually 
close to one atmosphere. These data have not 
been considered in our calculations. From such 
data he calculates 4” of our Eq. (51).!° The 
values so obtained are found to decrease rapidly 
with decreasing pressure. This corresponds to a 
pv curve with a negative (and large) coefficient 
of the second degree term, whereas his 1/L, data 
lead, as we have just seen, to a positive (and 
small) second degree coefficient. The two sets of 
data are therefore definitely incompatible, and 
since the average deviation of the pv data is 
much larger than for the L, data, we have 
discarded the former data entirely. Batuecas, 
on the other hand, considers that the decrease 
in the \4” values is due to experimental diffi- 
culties (surface effects, etc.) and adopts for \ 
the highest ,’’ value in each set. This happens 
to agree fairly well with his previous determina- 
tion of \. The average of all such highest values 
is 0.0245. 

Batuecas now uses his observed value of L 
and his graphically extrapolated value of Ljj,,, to 
get, by Eq. (4), \=0.0243. This he averages 
with his other value, 0.0245, and combines the 
new average \ with his observed L to get finally 
Liim = 2.25254. By chance this is close to our 
adopted Liin value, but we fail to see the 
justification for his handling of the pv data, or 
for his averaging of the results.'® 


‘6 Batuecas calculates the ratio of the change of px to 
that of p, using always the highest pressure (which is no! 
one atmosphere) for one of the points. This ratio is then 
divided by (pv):, a value he evidently gets by some sort 
of extrapolation. His resulting 4’ values are therefore 
not exactly those defined by Eq. (51), and the limiting 
value at p=0 is not exactly \, but the theoretical! dis- 
crepancy is negligible. 

16Tn particular one is never justified in combining an 
observed value of L with a value of \ obtained from ?% 
data in order to calculate Lym, when L is one of a set of 
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In order to leave no false impression regarding 
the properties of (CH;)Cl, we may note that 
later work on this gas has been performed bv 
Crespi.!? He obtains for L, at p=1, 1/2 and 1/4 
atmospheres the values 2.30734, 2.27998, and 
2.26617 respectively. These lie quite accurately 
on a linear graph, giving 


L, = (2.252489 +0.000062) 
+ (0.054868 +0.000093) p, 


s=0.43 X10~ or one part in 53,600 of L. 


Hence Liim = 2.252489 +0.000062, 
4 =0.024359+0.000041, 
L=2.30736, 

M =50.4886+0.002, 
C= 12.0082 +0.002. 


Because of the small number of points, the 
stated probable errors are rather unreliable, but 
it would appear that Crespi’s work is somewhat 
more accurate than that of Batuecas, and hence 
that the curvature in the L, plot of Batuecas’ 
data is due to some systematic experimental 
error. Fortunately both agree very well in the 
resulting values of Ziin and the atomic weight of 
carbon. 

Moles,'* in summarizing the work on (CH;)Cl, 
handles the data of Batuecas!® and Crespi!’ in 
an entirely different way. Although each of these 
investigations yields directly a value of Liim, 
and as just noted, practically the same value, 
Moles does not attempt to get a value of Liim 
from either. Instead he adopts Batuecas’ own 
average value of \ (0.0245), and uses our finite 
difference Eq. (51) to get from Crespi’s data a 


L, values. Such a procedure forces a slope on the L, plot 
that may be quite incompatible with that obtained directly 
from the L, data. In such a situation the only reliable Lj im 
value is that obtained directly from the L, data. In the 
language of statistics, the value of \ obtained from pv 
data may be considered as a priori information regarding 
the value of \ to be obtained from the L, data. A method 
for incorporating such information into the solution of the 
L, data for Liim could doubtless be devised, but the 
method would probably not be simple and so far as we 
know has never been investigated. 

" Crespi, Anales soc. espan. fis. quim. 25, 25 (1927). 

'SE. Moles, Rec. trav. chim. 48, 864 (1929). The Ly 
data of Crespi!? are listed here. 
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value }\=0.0244.'"° The average of these two 
(=0.02445) is then combined with Crespi’s 
observed value of L to obtain Liim. The resulting 
value (2.25227) differs appreciably from the 
value 2.2525 given directly by each investigation. 
This well illustrates the danger of working with 
an observed value of Z and a calculated value 
of \, rather than with a directly calculated 
value of Liim. Moles’ resulting value C= 12.003 
differs from ours solely because of the difference 
in the adopted values of Liim. 

It seems desirable to interpolate, at this point, 
a brief description of the method for calculating 
from (pv) data, ordinarily used by Batuecas 
and others who have been associated with 
Professor A. Guye. The method is described in 
detail in a paper by Batuecas, Maverick and 
Schlatter?® on the value of 1+ for oxygen. 

One obtains just three values of pv at different 
pressures, the first pressure being close to one 
atmosphere. These are fitted to a linear equation 
for po—our Eq. (55)—using a and 6 for our 
a’ and #’. The value of 1+. then follows from 
our Eq. (53’). The method used is really a 
variation of the finite difference method, Eq. 
(51), and is as follows. From points 1 and 2, 
using Eq. (55), 8’ is evaluated, i.e., 8’ = A(pv) /Ap. 
From points 1 and 3 a second value of #’ is 
obtained. The resulting average value of 6’ is 
substituted in Eq. (55) to give, for each experi- 
mental point, a value of a’. The three resulting 
values are averaged. Then our Eq. (53’) is used, 
with (a@’ aver +’ aver) in the denominator, but with 
each of the three values of a’ in the numerator, 
to get three values of 1+, which are finally 
averaged. The last two steps are obviously 
unnecessary. The average 1+ is given directly 
by awer/ aver +8’ aver): 

By making a graphical representation of this 
method one can see easily that it is equivalent 
to a least squares solution in which a high weight 
has been given to point 1 (at the highest pres- 
sure), a much smaller weight to point 2 and a 
still smaller weight to point 3. This is probably 
a more reasonable method of weighting than 


19 The values of are assumed by Moles to vary 
linearly, and the correct extrapolated value at p=0 is 
then actually 0.02444. 

20 T, Batuecas, G. Maverick and C. Schlatter, J. chim, 
phys. 22, 131 (1925). See p. 137. 
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that of equal weights. Hence we have no real 
objection to the method, although we feel that 
it is preferable to make a deliberate weighting 
of each point rather than to use this concealed 
and wholly arbitrary system of weighting.” 

We conclude this section with a discussion of 
a paper by Dietrichson, Orleman and Rubin,” 
entitled The Density of Ammonia at Reduced 
Pressures and its Relation to the Atomic Weight of 
Nitrogen, the Gas Constant R, and the Limiting 
Molal Volume, vo, which illustrates the confusion 
on this subject that still remains in the literature 
and the consequent need of an extended treat- 
ment of the theory such as we have just given. 
These investigators determine the density, p,, of 
ammonia at p=1, 2/3 and 1/3 atmospheres.”* 
They list and recalculate also previous data** on 
pp at p=1, 2/3, 1/2, and 1/3 atmospheres. 
D.O.R. give values of pv for w=1. Hence by 
our Eq. (2) this is merely 1/L, and in general 
we shall write 1/L, where they write pv.”> They 
also use d, for our L, vo for our R, and R for the 
gas constant which we shall denote Ro. (R= RoTo, 
or = RT> in their notation.) 
D.O.R. first write our Eq. (6) in the form 


L= M(1+ )/RoTo 


and remark that the accuracy with which M 
can be determined depends on the accuracy with 
which \, Ro and J) are known. We have, however, 
shown in detail that the value of \ is of only 
secondary importance. It is Ziim that is needed. 
But what is more important is that 7») has 
nothing to do with the problem. It is only the 
product Ro7>=R that is required. This product 


21 The question of weighting, and the actual experimental 
results of reference 20 will be discussed more fully in the 
second article of this series. 

2G. Dietrichson, C. W. Orleman and C. Rubin, J. 
Am. Chem. Soc. 55, 14 (1933). To be denoted D.O.R. 

23 Their actual results, which will often be used in this 
discussion, are pp=0.77126, 0.51161, and 0.25458 respec- 
tively. These lead by our Eq. (2) to Lp=0.77126, 0.767415 
and 0.76374 respectively, or 1/L,=1.296580, 1.303075 
and 1.309346 respectively. 

24E. Moles and T. Batuecas, Anales soc. espan. fis. 
quim 28, 871 (1930). To be denoted M.B. 

25 They also use 1 cc as the unit of volume and therefore 


all their values of pu and coefficients in equations for pu 


as a f(p) must be divided by 10* to obtain our 1/L, in 
liter-atmos, g.~' units, We make such a division in quoting 
their results, 
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R is known far more accurately than either 
factor,2® and as is shown in the next article of 
this series, has a probable error of only about 
one part in 64,000 (R=22.4146+0.0003). The 
final error in M is therefore due almost solely to 
the purely experimental error in Liim. 

D.O.R. next describe three different ways in 
which \ may be calculated, from the same set 
of data, and label the three different numerical 
results thus obtained \”’ and \’” respectively. 
Their \’ is gotten from our Eq. (62), in which 
as we have remarked, a” is a poor value of 1+. 
To get their a’’ they make a linear extrapolation 
of the L,/L, data to p=0. Their \” is gotten 
from our Eq. (34), in which again a’”’ is a poor 
value of 1/(1+2). Here also they make a linear 
extrapolation of the L,/L, data.?” Their first 
two methods then correspond to our cases (d) 
and (b) respectively, i.e., an assumption of 1/L, 
and L, linear in p, respectively. 

In the third method, D.O.R. say ‘‘the authors 
feel that the coefficient of deviation from Boyle's 
law should be defined as (1/pv) ps0L (pv) /0(p) Jr 
=)’’’.”” Since they use our Eq. (6) to get J, 
there is only one possible definition of 1+), 
namely the one given in our Eq. (4). Their 
definition of \’”’, as written, is ambiguous, but 
we find, by actual trial, that their published 
result (’’’=0.015109) can be obtained omly if 
one writes = —[1/(pv)o ][d(pv)/dp The 
true definition of \, by the discussion following 
our Eq. (48), is however 


—[1/(p2): (pv) /dp 0.014516. 


These results are based on their published 
equation 


pv =1/L, = 1.315403 —0.017767p— 0.001054". 


26 Actually Ry is determined from R/T» and has therefore 
practically the same probable error as has 7», i.e., one part 
in 9100. (T>=273.18+0.03 by G.C. 1929.) 

27 We have remarked in footnote 8 on the inadvisability 
of using any average first difference for such an extrapola- 
tion. Quite aside from this question, we are unable to 
check the results given by D.O.R. In the case of 1+)’ 
the two available first differences are 0.00501 and 0.004836, 
but to get their published 1.014901 one must assume 
0.005055 in the extrapolation. Similarly for 1/(1+’’) the 
available first differences are 0.004985 and 0.004765, but 
0.005031 is needed to get their published 0.985219 (=1/ 
1.015003). 
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In getting this equation the authors state that 
each point was weighted in proportion to the 
pressure and a least squares solution obtained. 
Since, however, there are three undetermined 
coefficients and only three points, the correct 
equation necessarily satisfies each point, regard- 
less of the weighting, and could hardly be called 
a least squares solution. Their published equation 
gives a small positive residual for every point. 
The correct solution, which can be obtained 
very quickly from a difference table, is 


1/L,=1.315390~0.017793p —0.001017p°. 


Having obtained their three different values 
of 4, which lead in turn to three different values 
of M, the molecular weight of NH:, the authors 
remark “It will be noted that the at.wt. of N 

. was found to show considerable variation. 
It therefore appears that the gas density method 
for the determination of atomic weights cannot 
at present be considered as reliable as that based 
on gravimetric analyses, not because of insuffi- 
cient accuracy as regards the gas densities, but 
rather because of uncertainties in connection 
with A, Ro and J» (also R=RoT>).”” But the 
variation in the value of \, as calculated in 
different ways, has nothing whatever to do with 
the value of R, and again we note that it is only 
the values of R and Lim that concern us in this 
problem. 

Having reached the above conclusion, D.O.R. 
now proceed to use our Eq. (3) to calculate R, 
from the accepted best value of M for NHs 
(17.0312). In this calculation they use, for the 
first time, the value of Liim directly calculated 
from the data. This value they calculate in two 
ways, (1) 1/Z, parabolic, giving their 1/Liim 
= 1.315403 in the equation already quoted, and 
(2) 1/L, linear, for which they get 1/Liim 
= 1.316484. From these two results they get 
R=22.4032 and 22.4213, respectively. They 
make similar calculations for the M.B. data, 
getting R=22.4136 and 22.4120. As we shall 
show immediately, the most probable value of 
Liim Obtainable from this work has a probable 
error of about one part in 12,000 and obviously 
cannot be used to improve the evaluation of R, 
which is already known to about one part in 
64,000, from the work on oxygen to be discussed 
in the second paper of this series. 
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Let us now turn to a consideration of the most 
probable value of Liim, and of M, resulting from 
the work of D.O.R. and of M.B. As just noted, 
D.O.R. have given solutions based on our Eqs. 
(44) and (52), i., 1/L, parabolic and 1/L, 
linear, respectively. We have made similar calcu- 
lations, with results somewhat different from 
those that they publish. We have already given 
our parabolic solution for the D.O.R. data. 
The corresponding linear solution with the 
weighting (proportional to pressure) used by 
D.O.R., is 


1/L,= (1.31582 +0.00011) 
— (0.01922+0.00011)p, 


s=0.50X10~‘ or one part in 26,000 of L. 

Liim = 0.759983 +0.000064 or one part in 12,000. 
(from Eq. (53)) = 0.01482 +0.00010, 

M =17.0347+0.0014, N= 14.0113+0.0014. 
For the reasons to be stated presently, we believe 


this to be the best result obtainable from the 
D.O.R. data. Their own equation is 


1/L, = 1.316484 —0.020066p 


and gives s=2.17X10-*, or 4.3 times that of 
the necessarily minimum value of the least 
squares solution. 

For the M.B. data we obtain as the linear 
solution 


1/L, = (1.315776 +0.000088) 
— (0.01991 +0.00012)p. 


s=0.60 X10~‘ or one part in 21,600 of L. 
Liim = 0.760007 +0.000051, 
4=0.015371+0.000089, 


M =17.0352+0.0012, N=14.0118+0.0012. 


D.O.R. give for the same treatment of the M.B. 
data 


1/L,=1.315940—0.020177p. 


We find from this equation, s=0.9610-‘, or 
1.6 times the true least squares result. 

Now a plot of the M.B. data shows no trend 
whatsoever from linearity. In fact an unweighted 
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parabolic solution of their data’ gives Liim 
= 1.315846, which differs from the linear result 
by less than the probable error of the latter. 
On the other hand, the value of s for the D.O.R. 
linear solution is even smaller than that for the 
M.B. data, so that unless the former data are 
by far the more accurate, the deviation from 
linearity should be ascribed to experimental 
error. Moreover, the value 1/Liim = 1.315390, as 
already given for the parabolic solution of the 
D.O.R. data, leads to M=17.0403, a result 
greater than the accepted value, 17.0312, by one 
part in 1870, whereas the linear solution gives a 
result greater by only one part in 4860. Finally, 
the linear solutions of the D.O.R. and the M.B. 
data give almost the same value of Liim. We 
conclude, therefore, that the mean deviation of 
the D.O.R. data from a linear curve, i.e., one 
part in 26,000, is to be ascribed to experimental 
error. We are not in a position to discuss critically 
the cause of the high resulting value of M. We 
may note, however, that if the D.O.R. point at 
p=1/3 is discarded, and a straight line passed 
through the other two 1/L, values (at p=1 and 
2/3), the result is Zjim = 0.759841, or M=17.0315 
and N=14.0081. It is therefore possible that 
systematic errors enter into the determination 
at p=1/3.29 

The value of Liim, as published by M.B.** is 
0.75990, as contrasted with our 0.760007 ob- 
tained from a linear 1/L, solution. We have 
made solutions of their data corresponding to all 
four of our cases (a) to (d), but none of the 


results is as low as theirs.2° From it, with their 


usual R= 22.414, M.B. obtain N =14.009. It is 
possible that M.B. obtained their result in a 
manner similar to that just mentioned in con- 


28 We find that the weighted and unweighted linear 
solutions give values of L)jm differing by only one part in 
200,000, so that the weighting is here immaterial. All our 
parabolic solutions have accordingly, to save time, been 
made with unweighted data. 

29 The observed p,=1/3 is given by D.O.R. as 0.25458 
+0.00001 g. liter! and the calculated value, on this new 
assumption, is 0.254536. The discrepancy is 4.4 times the 
stated probable error, but the data quoted on p. 16 of 
D.O.R. indicate the possibility of systematic errors of 
this magnitude. 

8° This value, 0.75990, is quoted by the Inter. At. Wt. 
Committee, J. Am. Chem. Soc. 53, 1627 (1931). We do 
not have available the original M.B. paper. 


nection with the D.O.R. work, for if one uses 
only the two 1/L, points at p=1 and 2/3, the 
resulting value of Liim is 0.759886, giving, with 
our value of R, N= 14.0091. So again there may 
be systematic errors at their two lower pressures 
(1/2 and 1/3 atmospheres). In any case it would 
appear that ammonia is not a suitable material 
to use in precision measurements of the at.wt. 
of nitrogen, by the limiting gas density method. 


CONCLUSIONS 


From the preceding discussion we may con- 
clude that, starting with the fundamental Eq. 
(3), the best value of Liim is obtained from the 
intercept of the L, curve, at p=0, after fitting 
the L, data to an analytic function by means of 
least squares. In the case of the permanent 
gases, a linear function in p is known to be 
satisfactory, in the region p= 1 to 0 atmospheres, 
but in the case of the easily condensible gases, 
it may be necessary to use a parabolic function. 
Further, if Z, is given by a linear function in 9, 
1/L, is likely to be an equally satisfactory linear 
function in ~, and this form has often been used 
in the literature. A similar remark applies to 
L, and 1/L, as parabolic functions in p. These 
four possible functions have been treated in 
detail in cases (a) to (d). In each of these cases 
there are several! ways in which the function 
may be written and handled, thus accounting for 
the wide variety of methods used in the litera- 
ture. In each of the cases (a) to (d) however, 
we have shown that there is only one method 
that leads directly to the most reliable value of 
Liim. Unfortunately some of the more commonly 
used methods of calculating Liim are by no 
means the most reliable, and numerical examples 
of this have been given in the concluding section 
on Illustrations. 

If one has only values of pv, at various pres- 
sures, it is possible to calculate only \ or Xo, as 
defined by Eq. (4). To get Liim we then need a 
value for L, the density at one atmosphere, as 
shown by Eq. (4’). It is, however, questionable 
procedure to use a value of L obtained from one 
investigation in combination with a value of }, 
or Xo, obtained in some other investigation. 
Furthermore, if experimental values of L, are 
available, these lead directly to a value of Liim: 
and the calculation of L and 1+. from such 
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data, as is often done, only complicates the 
matter, and has sometimes resulted finally in an 
untrustworthy value of Liim. These points have 
also been discussed in connection with cases (a) 
to (d), and have been illustrated by examples 
from the literature. It is, in fact, easily possible 
to obtain a value of Zin, and hence a molecular 
weight, that differs by several times its probable 


error from the most trustworthy value given by 
the data. Hence for data of the accuracy now 
obtainable in precision work on gas densities, 
the best analytic methods are not only suitable, 
but demanded. The actual application of such 
methods to the data on oxygen and on carbon 
will be given in detail in the next two papers of 
this series. 
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On the Statistical “Interaction Between Ions and Molecules” in Media of Small 
Dielectric Constant 


Otto HALPERN, New York University, AND Puitipp Gross, University of Vienna 
(Received November 24, 1933) 


In this paper the general theory of the interaction between ions and molecules in solution is 
discussed with particular reference to the theory of Debye on the salting-out effect. Formulae are 
developed for the free energy of the interaction between ions and dipoles. The formulae are 
applied to the dissociation equilibrium of an electrolyte in a medium of low dielectric constant. 


I. INTRODUCTION. COMPARISON WITH A THEORY 
OF DEBYE 


S is well known, the solubility of neutral 
molecules is very frequently changed by 
the presence of ions (salting-out effect). Debye! 
was the first to ascribe this phenomenon to the 
‘local changes in the average dielectric constant 
of the medium produced by the different at- 
tractions exerted upon the molecules of different 
types in the nonhomogeneous electric fields of 
the ions. From macroscopic electrostatics we 
learn that substances of high dielectric constant 
(high dipole moment) accumulate in regions of a 
higher field of force, while substances of lower 
dielectric constant are expelled therefrom. In 
perfect analogy, the molecules of higher dipole 
moment will accumulate in the neighborhood of 
the ions, while those of smaller dipole moment 
will mostly be found in more distant parts. Since 
the dipole moments of most solute substances are 
smaller than that of water these molecules will 
be repelled by the ions, or in other words the 
volume accessible to them will become dimin- 
ished by the presence of charges. This explains 
qualitatively the reduction of solubility of most 
solutes in water. Inversely, solutes having a 
dipole moment larger than that of water show 
an increase in solubility in the presence of 
charges, as we can predict by the simple picture, 
and as has been verified experimentally.? (Salt- 
ing-in effect.) 


1 P, Debye and I. McAulay, Phys. Zeits. 26, 23 (1925); 
P. Debye, Zeits. f. physik. Chemie 130, 56 (1927). 
2 Gross and Schwarz, Monatshefte f. Chemie 55, 287 


Debye has succeeded in expressing this picture 
in the form of a quantitative scheme. For this 
purpose he considered the influence of the various 
ions as being independent, and determined the 
free energy of a mixture the dielectric constant 
of which varies with the distance from the ion 
which is assumed as being at the center. In 
statistical equilibrium, the gradient of the 
osmotic pressure which is produced by the 
variable concentrations must be just balanced 
by the electric forces produced by the ion. He 
thus derives as a function of the distance from 
the ion the law of distribution into which the 
ionic diameter, the dielectric constants of the 
components, the charge of the ion and the 
absolute temperature enter as parameters. For 
practical evaluation some assumptions have to 
be made as to the variation of the dielectric 
constant with the relative concentrations of the 
molecules. 

There can be no doubt that the principle of 
Debye’s views is valid and leads to a correct 
understanding of the phenomena observed. On 
the other hand, some of the assumptions made 
for the purpose of calculation do not seem to fit 
all cases of importance. In aqueous solutions the 
effect generally is so small that we have to take 
ionic concentrations which are comparatively 
high (0.1 N and higher). In such solutions the 
average distance between the ions becomes s0 
small that we can no longer assume the effect of 
each ion to be independent of the presence of the 


(1930); Sitzber. Ak. Wien. IIb 139, 179 (1930); Gross 
and Iser, ibid. 55, 329 (1930); IIb 139, 221 (1930). 
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others.* Furthermore, it is an essential assump- 
tion on which Debye’s theory is based that we 
shall have at every point in the mixture a com- 
position corresponding to the equilibrium dis- 
tribution. This assumption is the more dangerous 
since the volume elements with which we are 
concerned are of a molecular order of magnitude, 
so that fluctuations will invariably play an 
important part and will influence the free energy 
of the composition decisively. The following cal- 
culations will show that this point is of para- 
mount importance, and it could have been 
suspected to play a decisive part on a priori 
grounds. In fact, it changes completely the 
quantitative form of the result for the case of a 
“salting-in’’ effect to which the formulae given 
below can be directly applied. 

Without any calculations it is easy to under- 
stand the qualitative difference between Debye’s 
treatment and that which we are going to 
propose, by considering the following example: 
The energy density of an ion carrying the charge 
ein a medium of a dielectric constant D in the 
distance r from the ion is given by the expression 


é/(82Dr'). 


The maximum change which can be reached by 
the presence of molecules of higher dipole 
moment is precisely equal to the value given 
above. This maximum change would correspond 
to an infinite dielectric constant. On the other 
hand, let us take a molecule with a dipole » and 
put it at a distance 7 from the ion. The work 
gained by carrying the molecule from infinity to 
r will increase with the dipole moment. The 
difference in the results is due to the fact that, 
in statistical equilibrium for the case of an 
ideally homogeneous complete mixture, the 
various dipoles influence each other, and in this 
way reduce the work which could be gained by 
the attraction of one single dipole. 

To sum up: Debye’s treatment would be rigor- 
ous if we were dealing with infinitely small 
molecules carrying infinitely small dipoles so 
that there is always a complete mixture even in 


= f 


a 


R 


*Gross, Monatshefte f. Chemie 53 and 54, 445 (1929). 
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molecular dimensions according to the laws of 
thermodynamics. We shall on the other hand 
neglect the ‘competition for space around an 
ion” which the different molecules are making 
each other. The formulas derived in the following 
paragraph therefore apply directly only to the 
case of a solute of very low concentration in a 
medium, the dielectric constant of which is 
negligibly small compared to the dielectric 
constant of the solute molecule. This permits us 
to consider the ‘‘fluctuation-terms” which are of 
decisive importance for dissociation—equilibria. 


II. StraTIsTICAL CALCULATION OF THE FREE 
ENERGY 


Let us consider the dilute solution of an elec- 
trolyte in a solvent of the small dielectric con- 
stant D. The electrolyte shall be partially dis- 
sociated, the concentrations expressed in mols 
per liter of the ions, and the undissociated mole- 
cules shall be denoted respectively by c; and 
cu. The charge of the ion, which shall be char- 
acterized by the index 7, shall be ¢;. The total 
volume of the solution shall be V, so that the 
volume which can be ascribed on the average to 
each ion will be equal to 


vi= 4rR°/3= 1000/N 


where N=Loschmidts number. The minimum 
distance between the molecules and the ion 7 
shall be aKR. We shall assume for the purpose 
of calculation that the value of the dielectric 
constant D shall be constant up to the surface 
of the ion, an assumption to be discussed in 
detail in due course. We now can divide the 
total volume into N>¢c; cells each containing one 
ion; and it shall furthermore be assumed that 
on the boundary of this volume the force exerted 
by the ion shall have become vanishingly small. 
Each volume v will contain Nceyv/1000 undis- 
sociated molecules. The excess free energy Wr of 
each molecule due to the interaction between ion 
and molecules is then given by the following 
equation 
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To evaluate this integral we expand the exponential function in a power series, and first carry out 
the integration over 3. We thus arrive at the following expression: 


"cos 
(n+1)! Jo 


a 


=4r 2dr} 1 


This series can be integrated term by term, with the result 


4rr 
M/kT — 


oo 1 1 eu 
( 
(2n+1)!(4n—1)\ DRT 


2n R 
a 


We will now discuss the value of the integrated 
power series separately for the two boundaries. 
For r=R the first term gives (47/3)R*=v while 
all following are smaller by a factor of the order 
of (eu/DkT R?)*. It is therefore sufficient to limit 
ourselves to the first term. 

The lower boundary necessitates a more elab- 
orate treatment. For practical reasons we are 
mostly interested in such values of the dipole 
moment and ionic diameter for which eu/DkT a’ 
becomes large compared to 1, since only then 
can an appreciable effect be expected. If we now 
rewrite the expression for our power series at the 
lower boundary in the following form 


(2n+2) '(1—7/2(2n+2)) 
(o = eu/DkTa) 


we notice that it differs from the power series 
for the exponential function e’ for three reasons. 
The first four terms are missing, the factorials 
in the denominator have the additional factor 
1 —7/2(2n+2) and every second term has dropped 
out. If now o is a large number compared to 1, 
as assufned, the value of the exponential function 
is almost exclusively determined by those terms 
in its power series the index of which is almost 
equal to a. In this case the first few missing terms 
do not contribute anything appreciable, and the 


(Qn +2) '(1—7/2(2n+2)) eu 


additional factor in the denominator can be 
replaced by a constant factor 1—7/2(20+2). 
(It will appear later on that it is without influence 
in any case.) For the missing terms we can 
account by introducing the factor 3. Keeping 
this in mind we can sum the power series for the 
lower boundary also, and thus arrive at the final 


expression 


rare’ 
(4a) 


mare’ 


(4b) 
vo*(1 —7/2(20+2)) 


Yu=—hkT log {1+ 


Since 
v=1000/NSc; 


we arrive at the following expression for the 
additional free energy of each molecule 


} 
100002(1—7/2(20+2))} 


Yu= log| 1+ 


For sufficiently low ionic concentrations we can 
expand the logarithm, and we thus obtain the 
following expression for the additional free energy 
of the total solution: 


100062(1—7/2(20-+2)) 
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IlI. REMARKS AND APPLICATIONS 


The most objectionable point in our derivation 
seems to be the use of the microscopic dielectric 
constant of the solvent up to the surface of the 
ion. For all practical applications where ionic 
diameter and the dipole of the molecule will be 
constants which must be determined empirically 
we have to keep in mind that the values thus 
obtained are not strictly correct, but might be 
considerably influenced by the probably wrong 
assumption as to the value of the dielectric 
constant. At the present stage we do not think 
that matters can be improved in this respect. 

The application to the change of solubility is 
quite obvious. The equilibrium condition for the 
molecule in solution can be simply obtained by 
differentiating VY, with respect to Ncy and adding 
the resulting expression ~y to the free energy 
of the solute. The dissolved molecules have an 
activity coefficient e¥”/*” the logarithm of 
which is proportional to the concentration of the 
ions, and which depends in a very sensitive 
way on the dielectric constant of the solvent. 

A more important application, for which our 
treatment would seem to be particularly appro- 
priate, is the dissociation equilibrium, and its 
dependence on concentration. In media of low 
dielectric constant the effect calculated above 
will be of great importance, and might account, 
as we have anticipated, for a number of observed 
anomalies. The equilibrium condition between 
ions and non-dissociated molecules can be written 
in the form 

f2c?/fucm=K (7) 


fu stands for the activity coefficient of the 
molecule, f; for that of the ion. f; is due to the 
interionic effects as well as to the interaction 
between ions and molecules. We shall here omit 
the interionic contribution to the activity coef- 
ficient, since it can always be inserted easily, 
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provided that the theory of ionic interaction has 
been developed for the solvent under consider- 
ation. The contribution to the activity coefficient 
resulting from the ion-molecule interaction is 


obtained by differentiating YW, with respect to 


' New and Nce;. Inserting these expressions into 


(7), we arrive at the following equilibrium con- 
dition 
c? 


log —=log K 
Cm 


kT Mee 
10000°(1 —7/2(20+2)) 


(8) 


In the special case of slight dissociation (8) sim- 
plifies to 


log A+ Be (9) 


(a=degree of dissociation and c= analytic con- 
centration of solute). 

This equilibrium condition shows interesting 
features, insofar as the dissociation sometimes 
increases with increasing concentration, because 
of the presence of the term Be on the right-hand 
side. The concentration for which the dissociation 
does not change with concentration is given by 
the relation 


da/dc=0, 
which can be written in the form 


1/Cmin.= B. (10) 


If we substitute this expression for B in (9) we 
notice that the change of dissociation with con- 
centration is very unsymmetrical on both sides 
of the point Cnin., and that more particularly 
the dissociation increases very rapidly if we 
Pass Cmin. 

We shall make use of the formulae here derived 
in the following paper, to explain conductivity 
phenomena in media of low dielectric constant. 
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The authors give a general summary of some recent 
theoretical views concerning the characteristic features 
of electrolytes in media of small dielectric constant. 
After a discussion of conductivity observations in water- 
dioxane mixtures, previous attempts for their theoretical 


interpretation are criticized. It is shown that in agreement 
with previously proposed views of the authors the experi- 
mental data can be qualitatively and partly quantitatively 
understood by considering the activity-coefficient arising 
from ion-dipole interaction. (See preceding paper.) 


I. RECENT CONDUCTIVITY OBSERVATIONS AND 
THEIR INTERPRETATIONS 


HE behavior of electrolytes in solvents of 

low dielectric constant is known to show a 
number of anomalies when compared with 
aqueous solutions. The electrolytes are appar- 
ently of high molecular weight when judged by 
osmotic phenomena. At the same time these 
seemingly highly associated electrolytes show 
comparatively high electric conductivity. This 
conductivity does not depend on concentration 
in the way to which we are accustomed in the 
case of water, or media of considerable dielectric 
constant. Here the equivalent conductance 
decreases with increasing concentration, in 
accordance with Ostwald’s law of dilution or 
with the theory of the conductivity coefficient 
(interionic forces). The authors! have some time 
ago (1925) given a general discussion in which 
these various contradictory phenomena are 
brought into harmony. By assuming interaction 
forces between the ions and between the ions and 
molecules they provided a suitable explanation 
for the apparent contradiction between “high 
association’’ on the one hand and anomalous 
conductivity on the other. The most powerful 


1 Phys. Zeits. 26, 636 (1925). 

The authors wish to make use of this opportunity to 
correct a point in their previous presentation. In the paper 
quoted above they made use of an expression for the 
activity coefficient of the ions as given by Debye’s theory 
of electrolytes. This application is certainly untenable, 
since the behavior of ions in media of small dielectric 
constant cannot even approximately be described by 
Debye’s formula. 


tool for treating these problems was offered by 
the conditions for thermodynamic stability.’ 

Quite recently, several papers have appeared 
which are of importance for the problem under 
consideration. An investigation by V. K. La Mer 
and Harold C. Downes? has again drawn atten- 
tion to the appearance of anomalous conductivity 
of acids in benzene. These authors explain their 
results by making use of the views which have 
been given by the writers. Still later in a paper by 
Kraus and Fuoss‘ the conductivity of tetraiso- 
amylammonium nitrate (and some other salts) 
has been studied with a solvent composed of 
dioxane and water in shifting proportions. 
Furthermore in several theoretical investigations 
the activity coefficients which are due to inter- 
action between dipoles and ions have been 
studied by Fuoss and Kraus.® In this paper we 
shall mostly be concerned with the theoretical 
interpretation of the experimental results ob- 
tained by Fuoss and Kraus. 


2 Rodebush (Rodebush and Ewart, J.A.C.S. 54, 419 
(1932)) claims that the ideas presented above have been 
anticipated in a paper by himself and Latimer (J.A.C.S. 42, 
1419 (1920)). We readily admit that Rodebush and 
Latimer have in a general way mentioned the influence of 
electric forces the existence of which was in fact obvious. 
The determining influence of the activity coefficient on 
anomalous conductivity, the connection between anomalous 
conductivity and small osmotic pressure, and the general 
thermodynamic relations were not mentioned. 

3 La Mer and Downes, J.A.C.S. 53, 891 (1931). 

4 Kraus and Fuoss, J.A.C.S. 55, 21 (1933). 

5 Fuoss and Kraus, J.A.C.S. 55, 476 (1933). 

Fuoss and Kraus, J.A.C.S. 55, 1019 (1933). 
Fuoss and Kraus, J.A.C.S. 55, 2387 (1933). 
Fuoss, Zeits. f. Elektrochemie 39, 513 (1933). 
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The experimental results obtained by these 
two authors can perhaps be simply summarized 
as follows: The conductivity of tetraisoamylam- 
monium nitrate is rather small at most concen- 
trations, if pure dioxane or mixtures with a small 
percentage of water are used as solvents. If the 
solvent contains a large percentage of water the 
conductivity (i.e., the dissociation) increases very 
rapidly. Keeping the solvent constant the equiv- 
alent conductivity shows a minimum with 
increasing concentration, and increases fairly 
rapidly after passing this minimum. The location 
of this minimum shifts markedly to a higher 
concentration with increasing dielectric constant 
(i.e., higher water percentage of the solvent). 
For sufficiently high dielectric constant the 
equivalent conductivity shows no more minima, 
and changes only slightly with concentration. 

Fuoss and Kraus have attempted to explain 
these observed facts by making use of the con- 
cept of “ionic association’? which has been 
introduced into the theory of electrolytes by 
Bjerrum.® According to this author ions which 
are at a distance R from each other smaller than 
€€./2DkRT (e:, charge of ion, D, dielectric 
constant, k, Boltzmann’s constant, 7, abs. 
temp) shall be considered as forming a neutral 
molecule which contributes neither to the inter- 
ionic forces nor according to Fuoss and Kraus to 
the electrical conductance. As Bjerrum already 
mentioned, and as actually worked out by Fuoss 
and Kraus, in their paper IV, this formation of 
complexes was not considered to cease with 
molecular formation, but should proceed to 
complexes of higher type (triple ions, quadruple 
ions, and so on). Thus, for example, a dipole 
carrying molecule and an ion were considered to 
form a “triple ion,” provided the distance R’ 
between them is less than a quantity R’ defined 
by the relation 


1/R’=2DkT/2+R’/(R’' +5)’, 
where 5=dipole moment u divided by e. 


With the help of Bjerrum’s formulae Fuoss and 
Kraus calculated the ‘‘dissociation constant’’ for 


‘Bjerrum, Copenhagen, Akad. fys. Medd. 7, 9 (1926). 


the various equilibria 
crfs/eu= Ki, 
f= Keo. 


(cy, Cu, Cy denote respectively the concentrations 
of free ions, ‘“‘molecules’’ and triple-ions; f; and 
fx are the ionic activity coefficients due to 
interionic forces.) Into these formulae, param- 
eters enter which denote the minimum distance 
between two ions, between an ion and a molecule, 
and so on. These were so determined as to give 
the best fit with their observation. In this way 
the authors were able to account for the increase 
of the dissociation with increasing dielectric 
constant, for the minimum in the conductivity- 
concentration curve, and for the shift of this 
minimum to higher values of the concentration 
with increasing dielectric constant. The authors 
furthermore made use of the Debye formula for 
the activity coefficient of the ions arising from 
the interionic forces; since the ionic concentration 
is very small it would follow (assuming correct- 
ness and applicability of this formula) that the 
interionic effect is of secondary importance 
only.” 

The main points which Fuoss and Kraus 
consider as experimental verifications of Bjer- 
rum’s ionic association may be enumerated as 
follows: (a) With a suitably chosen value a; for 
the minimum distance of approach between the 
ions it is possible to calculate the variation of the 
logarithm of the molecular dissociation constant 
for changes in the dielectric constant of the 
medium; the following relation is found to hold: 


log K,~1/D. 


(b) By introducing another molecular diameter 


7 Into the formula for the activity coefficient there enters 
the actual ionic concentration to be determined from the 
analytic concentration with the aid of the law of dilution 
which itself contains the activity coefficients. The authors 
Fuoss and Kraus quote, as their sources, papers by 
Davis, J. Phys. Chem. 29, 977 (1925). 

Sherill and Noyes, J.A.C.S. 48, 1861 (1926). 

MaclInnes, J.A.C.S. 48, 2068 (1926). 

It may be mentioned that a discussion of the influence of 
the activity coefficient and the law of dilution for weak 
electrolytes has been given in a paper by Gross and 
Halpern, Phys. Zeits. 25, 393 (1924). 
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dz they calculated a dissociation constant for the 
equilibrium 
Ko. 


This dissociation constant K» permits the deter- 
mination of the shift of the minimum in the 
conductance-concentration curve with changing 
dielectric constant of the solvent. The authors 
furthermore claim that there exists a range in 
which the product D?/cnin. is constant according 
to their theory of interaction between molecules 
and ions (triple ion formation). This they con- 
sider to be a theoretical explanation of an 
empirical law previously found by Walden. 

For the reasons given in the following dis- 
cussion we do not consider an explanation given 
on the basis of Bjerrum’s theory of ionic associ- 
ation to be tenable. 


II. BJERRUM’s CONCEPT OF IONIC ASSOCIATION 
AND Its APPLICATION TO CONDUC- 
TIVITY PHENOMENA 


In a previous paper® by one of the writers a 
discussion of Bjerrum’s ‘‘ionic association”’ has 
been given to which the reader should refer for 
details. It has there been shown that the Helm- 
holtz-relation between energy and free energy 
of the ions is not satisfied in Bjerrum’s theory. 
This result is due to the fact that the maximum 
distance at which two ions are considered as 
forming a molecule depends on the temperature 
according to the relation R= e;e;,/2DkT. Notice 
that R depends on T and D only in the form 
R=f(D-T); since this is a general consequence 
of statistical mechanics elimination of the tem- 
perature-dependence by some modification of the 
expressions for R would likewise eliminate the 
influence of a variation of D and thereby com- 
pletely invalidate the application of the theory 
to media of low dielectric constants. The same 
objections must of course be raised against the 
calculations of higher complexes (triple-ion 
equilibria etc.). 

A correction which is of slight importance in 
the case of aqueous solutions (or in the mathe- 
matically equivalent case of van der Waals 
forces in non-ideal gases®) becomes of decisive 
importance for media of small dielectric constant. 


8Q. Halpern, J. Chem. Phys. 2, 85 (1934). 
9 Cf. e.g., R. H. Fowler, Statistical Mechanics, § 9.6, 


Here we are referring to the determination of 
the volume which is accessible to the free ions. 
According to the general concept of ionic associ- 
ation an ion is considered to be free if its distance 
from the closest ion of opposite sign is larger 
than R. For determining the free energy of the 
ions one therefore has to consider them as con- 
tained in a volume which is equal to the total 
volume minus the ‘‘volume occupied by the 
molecules.”’ (This corresponds to the correction 
for the proper-volume in van der Waals theory). 
The ‘molecules’ in Bjerrum’s theory of ionic 
association have a “‘radius’”’ R which for media 
of low dielectric constant (e.g., dioxane D= 2.2) 
takes a value 1.3X10~® cm. Since almost all the 
salt dissolved is forming ‘‘molecules” the volume 
accessible to the ‘‘free ions’” becomes very much 
reduced (sometimes even negative); the free 
energy of the ions therefore is very much in- 
creased as compared with the value it would 
take if the total volume were accessible to the 
ions. 

To illustrate this fact let us consider an 
actual case of practical interest: a normal 
solution of tetraisoamylammonium nitrate in 
dioxane (D= 2.2). With the value R we obtain 
for the volume per liter occupied by the “mole- 
cules’’ 5000 cc. There is no rational way of cor- 
recting for these quite enormous errors. They are 
fundamentally associated with the unsound basis 
of the theory which treats long range (Coulomb) 
forces like short range (van der Waals) forces. 

Furthermore Fuoss and Kraus have to 
assume that these associated ions do not con- 
tribute anything to the conductivity of the solu- 
tion. They have not given any reasons to support 
this contention, nor is it by any means intuitively 
obvious why two ions with some forty molecular 
layers of liquid between them and_ negative 
potential energy (2k7)) only slightly in excess of 
their average kinetic energy should not par- 
ticipate in electric conduction. On the other 
hand, they have to assume for example that a 
molecule and an ion are bound to a triple ion 
even though the binding energy in some cases is 
smaller than kT, that is to say to less than two- 
thirds of the average kinetic energy. 

The idea in the papers of Fuoss and Kraus 
of explaining conductivity curves by partial 
dissociation is not new. What is new is only the 
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ELECTROLYTES IN MEDIAL OF 
interpretation of partial dissociation as ionic 
association in Bjerrum’s sense. This method of 
computation we have already criticized from a 
theoretical point of view. The fact that the 
logarithm of the dissociation constant varies 
proportional to 1/D, and that this variation is 
implied in Bjerrum’s picture cannot be thought 
to corroborate the latter. Every model will show 
the same dependence, because the work necessary 
to remove two ions from each other varies as 
1/D. Deviation from this simple relation could 
only be expected if the ions had a different 
solvation energy in the various media. Since 
some traces of water appear to be always present, 
the ions would be mostly surrounded by water 
dipoles (salting-out effect), so that the energy 
of solvation probably varies but slightly. This 
view is supported by experimental evidence 
gained through observations of the index of 
absorption of ions. The index of absorption of 
the anion of picric acid dissolved in pure ethyl 
alcohol (D=26) shifts to the value which it 
shows for aqueous solutions if only a small 
amount of water (0.7 m/lit.) is added to the 
alcohol. 

Remembering the difficulties introduced by 
the contradiction with thermodynamics, by the 
impossibility of a correct statistics (volume- 
correction) and by the crude assumptions of 
widely distant but nevertheless nonconducting 
ion-pairs, we feel justified in rejecting the quan- 
titative application of the hypothesis of ionic 
association. 


Il. Tue INFLUENCE or THE Activity-CoEr- 
FICIENT ARISING FROM THE ION-DIPOLE 
INTERACTION 


As mentioned in the introduction, it has been 
pointed out by the authors that anomalous con- 
ductivity can be explained by the influence of 
the activity coefficient due to the interaction of 
the ions with each other and with neutral 
molecules. The ion-dipole interaction we have 
treated in the preceding note; the formulae for 
the dissociation-equilibrium and for the concen- 
tration of minimum dissociation apply directly 
to the present case. 

The evaluation of the activity-coefficient due 
to interionic forces is at present impossible. It 


“Ph. Gross, Zeits. f. Elektrochemie 36, 789 (1930). 
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6 5 44345 5 25 


0.2 0.3 0.4 

Fic. 1. Test of Eqs. (8, 9, 10) of the preceding paper. 
The connection between minimum-conductance concen- 
tration (Cmin.) and dielectric constant of the medium. 
Tetraisoamylammonium nitrate dissolved in various 
dioxane-water mixtures. 


has been shown® that in media of small dielectric 
constant the Debye-formula for the activity- 
coefficient does not give the right order of magnitude 
even at the lowest concentrations which are accessible 
to observation. From the curves given below we 
learn that it is possible to reproduce quan- 
titatively several features of the conductivity 
phenomena by tentatively neglecting the inter- 
ionic effects and merely considering the ion-dipole 
interaction. This very crude procedure cannot of 
course claim any accuracy; but it is significant 
that the agreement is best for very low dielectric 
constants and near the minimum of dissociation, 
in which cases the ionic concentration becomes 
very small. 

Fig. 1 has been constructed with the aid of 
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2 4 6 10 
c (105) 


Fic. 2. log \°c as a function of the concentration c of tet- 
raisoamylammonium nitrate dissolved in pure dioxane. 
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formulae (8), (9) and (10) of the preceding paper. 
The degree of dissociation is replaced by the 
equivalent conductivity \ to which it is here 
assumed to be (nearly) proportional. For values 
of the constants we have taken a= 4.2 cm 
and w=5.8X10-'8 cm e. s. u. In Fig. 2 the ob- 


served values of log Mc are compared with the 


theoretical values given by (8) and (9) of the 
preceding paper. 

This agreement is much less satisfactory than 
it appears, and is partly misleading, since we 
are not able to reproduce theoretically by means 
of (8) of the preceding paper the whole conduc- 
tivity-concentration curve for all media. Our 
formula gives satisfactory results principally for 
the region close to the minimum-conductivity 
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and for the case of very small dielectric con- 
stants. Otherwise the activity coefficient due to 
interionic forces becomes important. We have 
furthermore already pointed out in the preceding 
paper that the meaning of the dielectric constant 
at the surface of the ion becomes somewhat am- 
biguous. 

The formulae of Kraus and Fuoss can only 
be considered to be empirical laws. Their claim 
to have provided a theoretical basis for Walden’s 


D*/Cmin.= const. 
does not seem tenable, since this law relates to 


ionic concentrations which are so high that the 
interionic forces become dominant. 
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Acoustical Studies. III. The Rates of Excitation of Vibrational Energy in Carbon 
Dioxide, Carbon Disulfide and Sulfur Dioxide’ 


T. RicHARDS AND JAMEs A, Reip, Princeton University 
(Received January 6, 1934) 


Measurements on these three gases at various pressures, 
temperatures and frequencies are reported. In all cases 
double collisions provide the predominant mechanism of 
excitation. The van der. Waals or capillary forces follow 
the acoustical cycle without lag under the conditions of 
the experiments. The velocity of sound at 9 ke in carbon 
dioxide shows that a part of the heat capacity has already 
disappeared from the adiabatic elasticity at this frequency, 
and comparison with the measurements of other investi- 
gators leads to the belief that this term is missing even at 
considerably lower frequencies at ordinary temperatures. 
The remaining vibrational energy terms disappear together 
from the acoustical cycle, their rates of excitation being 
sensibly equal. There is a considerable activation energy 
necessary to excite the deformation vibration. The velocity 
of sound at 9 ke in carbon disulfide agrees well with that 
obtained by calculation from the results of band spectral 
analysis-and hence no important energy terms appear to 


be lacking at this frequency. The symmetrical linear 
and the deformation (S|CfS]) vibrations 
leave the acoustical cycle together at high frequencies or 
at low pressures. A study of the dispersive region caused by 
this disappearance shows that the rates of excitation of the 
various energy states are identical within limit of error, 
and consequently that the activation energy of collision is 
uniform throughout the dispersive region. The velocity of 
sound at 9 kc in sulfur dioxide conforms reasonably well 
to the somewhat conflicting band-spectral data for the 
equilibrium heat capacity. There is a serious discrepancy 
between the results of the velocity measurements in the 
dispersive region here reported and the current interpreta- 
tion of the few existing absorption measurements. A 
tentative argument is given whereby this discrepancy may 
be removed. Several theoretical questions which require 
further elucidation are discussed. 


I. METHOD OF EXPERIMENT AND 
OF PRESENTATION 


A. Apparatus 

The measuring instruments were modifications 
of the sonic interferometer devised by Pierce.? 
Since they have already been described in a 
communication from this laboratory,* only the 
briefest sketch is here necessary. Magneto- 
striction oscillators, cut from nichrome steel, 
were mounted in cylindrical tubes of Pyrex glass 
and driven by valve-tube circuits. Points of 
maximum reaction in the vibrating column of 
gas were explored with plane-ground Pyrex 
reflectors, and detected by means of a micro- 
ammeter in the plate circuit of the driving 
oscillator. The displacements of the reflecting 
pistons were measured to 0.05 mm from 9 to 92 


‘An advance notice of this work appeared in Nature 
130, 739 (1932). A few of the earlier conclusions have here 
been modified owing to the introduction of a more quanti- 
tative method for the analysis of experiment. 

* Pierce, Proc. Am. Acad. Sci. 61, 1 (1928). 

’ Richards and Reid, J. Chem. Phys. 1, 114 (1933). 


ke and to 0.005 mm from 94 to 451 kc. The 
frequency was measured to 0.05 percent at all 
temperatures; it is not necessary to quote here its 
exact value in each case. Determinations at 9 and 
92 kc were made exclusively in Pyrex tubes of 2.5 
cm internal diameter and of 1 mm wall-thickness. 
Determinations at 94 and 451 kc were made 
exclusively in Pyrex tubes of 1.0 cm internal 
diameter and 1 mm wall-thickness. The pressure 
and temperature control were suitable to repro- 
duce the velocity of sound in a non-absorbing gas 
to +0.05 m 

The use of tubes of small diameter greatly 
simplifies the process of measurement, but 
necessitates a tube-correction to compensate for 
the retarding influence of the walls on the com- 
pressional wave. The tube corrections for this 
apparatus have already been defined by measure- 
ments in pentane and air at 9 and 92 ke and in 
pentane, propane and air at 94 and 451 ke. A 
number of measurements were made on other 
gases, but since they do not alter the previous 
values they will not be printed. In the meantime 
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a paper has appeared by Kaye and Sherratt‘ and 
an investigation has been undertaken in the 
laboratory of Professor Pierce by Mr. G. A. 
Norton® which very strongly support the Helm- 
holtz-Kirchhoff reduction formula, and_ place 
measurements of the velocity of sound in small 
tubes in an extremely strong position. 

As before the measurements at 9 and at 92 ke 
have been made compatible by a tube correction, 
and the measurements at 94 and 451 kc have been 
similarly adjusted to give mutual concordance. 
No attempt has been made to correlate the 92 
and 94 kc measurements since a difference of a 
few tenths of a meter in the absolute velocities is 
considered of no importance. The tube cor- 
rections have already been applied to all the 
measurements quoted below. 


B. “Idealization” of the measured velocity of 
sound 
The velocity of a plane sinusoidal wave of 
infinitesimal amplitude propagated in an ideal, 
non-absorbing gas of molecular weight \/ and 
heat capacity at constant volume C is given by 


Viaeat — ((RT/M)(1 +R/C) } i, (1) 


This will be spoken of below as the ideal velocity of 
sound. 

In interpreting measurements on the free-space 
velocity of sound in real gases it has been 
customary to make use of this expression, and 
consequently to neglect two sets of corrections. 
The first set deals with the acoustical conditions, 
and includes corrections for the finite amplitude 
of the wave, viscosity and heat conduction of the 
gas, and the like. No attempt to apply corrections 
of this character will be made, since it is believed 
that they are of secondary importance for the 
interpretation of the type of dispersive effect 
which is here discussed. The second set of 
corrections result from non-ideal behavior of the 
gas. Since it is expedient, in order to economize 
space, to plot measurements at various pressures 
and frequencies, or at various temperatures and 
frequencies, on a single Cartesian surface, cor- 
rections of the second type become imperative. 


4 Kaye and Sherratt, Proc. Roy. Soc. A141, 123 (1933). 
°G, A. Norton, private communication. 
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They will be introduced by means of a factor w 
which is to be applied in the sense w= V?jdeu1/ 
V? eae It has been obtained by the following 
argument. 

If the van der Waals equation 


(P+a/v?)(v—b) = RT (2) 
is used for the description of a real gas of density 
p the velocity of sound is found to be 


RT+P(b-—2a/RT) R RT 
-| + --. (3) 
M-—pb C M—2pb 


if appropriate derivatives for substitution into 
the equation V=(0P/dp)*saia, are evaluated, 
neglecting terms in p*. Eq. (3) is suitable for the 
description of the velocity of sound at low 
frequency in a real gas, and may easily be 
adapted for use with a dispersive gas in which a 
part of the heat capacity Cy—C., fails to follow 
the adiabatic acoustical cycle at high frequencies, 
if the assumption is made that the van der Waals 
forces remain in equilibrium. Since according to 
(2) (0C,/dv)r=0, the heat capacity and_ the 
density remain independent variables, and the 
usual derivation according to the method of 
Einstein is not seriously disturbed. This yields 


RT +P(b—2a/RT) 
M-—pb 
RT 
(M—2pb) 


= 


| (4) 


and is suitable for the description of a real 
dispersive gas which has a heat capacity of Cy 
at equilibrium, of C,, at frequencies so high that 
wi? =e, and in which the energy states con- 
tributing Cy—C,, to the heat capacity have 
identical approximate relaxation times J. The 
factor uw is determined by the absorption.® 

It is, therefore, the ratio of the squared 
velocity according to (1) and that calculated 
from (3) which determines the van der Waals 
correction factor w for the low frequency 
velocity of sound. Corresponding values of w in 


6 The significance of » and of 8 have been discussed by 
Richards in J. Phys. Chem. 1 (1933).on pages 865 (Eq. 
1.8) and 876 (Eqs. 5.5 and 5.6), respectively. 
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the dispersive region may readily be obtained 
from the customary ideal dispersion equation 
developed from (1) and from its corrected 
equivalent (4). Measured velocities of sound 
which have been multiplied by a calculated 
quantity w’ will be spoken of as idealized 
measured values. 

A considerable advantage is gained by ideal- 
izing measured velocities if the energy levels of 
the molecule are known from band-spectral 
analysis. In this case (1) contains no unknown 
quantities, and the ideal value of the velocity of 
sound at very low frequency is defined with great 
accuracy. Any lack of correspondence with 
experiment at low frequencies becomes, there- 
fore, of immediate interest. 

In an ideal gas, the energy adiustments of 
which are effected by double collisions only, an 
increase of the acoustical frequency has the same 
effect on the velocity as a numerically equal 
decrease of pressure. This is, of course, the result 
of the fact that the number of double collisions 
per molecule is directly proportional to the 
pressure. Hence a given value of v/P, where r is 
the acoustical frequency, defines the velocity at 
any given temperature. A plot of V? against 
logio v/P is a symmetrical S-curve if the gas has 
the dispersive properties described by (4). This is 
easily constructed and any experimental devia- 
tions from the regular function will be readily 
apparent. 


C. The effect of traces of gaseous impurities 


If there is a steric or energetic hinderance in 
the adjustment of vibrational energy by collisions 
between like molecules, very small traces of 
impurities may greatly affect the transition 
probability which is inferred from experiment. 
Since irregularities are reported below which 
appear to have this cause, it is perhaps worth 
while to cite a simplified specific example. 

Let us consider a pure gas A at such a tempera- 
ture and pressure that each molecule suffers on an 
average 10" double collisions per second. Let us 
further suppose that on an average only one 
collision in 10° is effective in bringing about a 
given energy transition. This means that, for the 
transition in question, faa is 10* sec.—'. Let us 
now introduce 0.1 percent of a gas B which has 
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approximately the same molecular weight and 
heat capacity as A and hence does not appreci- 
ably alter the measured value of the velocities at 
very low and at very high frequencies. The 
probability that the energy transition may be 
effected is now faa The number of collisions 
with B molecules suffered by an A molecule is 
about 10’ per second. If each one of these is 
effective fin is 107 sec.~' and is so far greater 
than faa which remains 10* sec.~! that the 
dispersive region in the gas will shift to fre- 
quencies which are higher by almost three orders 
of magnitude. If, on the other hand, the efficiency 
of an AB collision is equal to or less than the 
efficiency of an AA collision fig becomes nu- 
merically far inferior to f44 and the presence of 
the trace of gas B will cause no noticeable effect. 

Traces of accidentally included impurities are, 
therefore, not likely to cause an apparent 
diminution of fi, but may cause an enormous 
apparent increase. The lowest value observed for 
faa is, therefore, most likely to be correct if the 
gas is known to be reasonably pure. 

Since an ordinary glass surface retains water in 
the adsorbed condition which will evaporate 
slowly into any pure gas in contact with the 
surface, irregularities are to be expected if water 
is an effective agent in bringing about energy 
transitions in the pure gas. Evidence is already 
at hand which indicates that collisions between 
water molecules and carbon dioxide molecules 
or oxygen molecules give rise to this effect.7? In 
such cases wholly trustworthy values of the AA 
probabilities may evidently be obtained only by 
the use of high vacuum technique. 


II. MEASUREMENTS ON CARBON DIOXIDE 


The carbon dioxide was prepared from calcium 
carbonate and a guaranteed brand of hydrochloric 
acid. It was passed through a drying train of 
calcium oxide and phosphorous pentoxide, and 
was, therefore, reasonably pure when admitted 
to the evacuated apparatus. An oil pump of 
conventional design was used for evacuation, and 
the apparatus was washed several times with 
pure gas before a series of measurements was 


7 Carbon dioxide, Rogers, Phys. Rev. 41, 369 (1932); 
Eucken and Becker, Zeits. f. physik. Chemie B20, 467 
(1933); oxygen, Kneser, J. Acous. Soc. Am. 5, 122 (1933). 
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made. It is significant that, even with these 
precautions, it was difficult to obtain repro- 
ducible results. Evidently traces of some sub- 
stance which is extremely effective in adjusting 
the transverse vibration were not eliminated by 
this preliminary treatment. 


A. The velocity of sound at 9 kc 

The measured velocity of sound at 9 kc, 30°C, 
and 774 mm was 271.8 m sec.~', and at 403 mm 
was 272.2 m sec.-! at the same temperature. 
When idealized these yield an average of 272.5 m 
sec.—' at 9 kc, and show no trace of dispersion in 
this frequency and pressure region. The ideal 
velocity of sound at very low frequency from (1) 
may in this case be calculated with some 
precision, since the energy levels of the COs 
molecule have been extensively studied. Strictly 
speaking the unperturbed Planck-Einstein ¢- 
function for the vibrational energy 


Co/R=5/2+2¢01/T)+ (5) 


must be modified, since the second energy 
state corresponding to the _ characteristic 
temperature 0,(0 | CTO) is so close to the 
first energy state of 02(O0<C-—O) that a per- 
turbation occurs in one of the three degenerate 
levels of the second state of the 0; vibration. The 
heat capacity is somewhat larger than that given 
by (5) in consequence. The error is, however, 
small, and there is reason to believe it inap- 
preciable at 30°C from the recent analysis of 
Adel and Dennison,’ and the approximate 
agreement at low temperatures of the measured 
heat capacity of carbon dioxide with that 
calculated from the Planck-Einstein function 
found by Eucken and Miicke.* The value of 
C,/R will, therefore, be considered to be 3.5116 
at 30°C, which is the figure resulting from (5) if 
0, = 956°, 02.=1841°, and 0;=3331°. This corre- 
sponds to an ideal low frequency velocity of 
271.2 m sec.~', and is significantly lower than the 
corresponding idealized experimental values at 
9 kc. The discrepancy strongly suggests that a 
part of the heat capacity has already disappeared 
from the adiabatic cycle at acoustical frequencies 
considerably less and at frequencies of molecular 


8 Adel and Dennison, Phys. Rev. 43, 716 (1933). 
®Eucken and Miicke, Zeits. f. physik. Chemie B18, 
186 (1932). 


impact considerably greater than those of the 9 
kc measurements here reported. 

A certain amount of corroborative evidence for 
this conclusion is already at hand. The Inter- 
national Critical Tables give two empirical 
formulae for the velocity of sound in carbon 
dioxide over a range of temperatures, which were 
set up from the measurements of Buckendahl!” 
and of Dixon, Campbell and Parker" at at- 
mospheric pressure. Both these agree that at 
30°C and 760 mm the velocity of sound is 271.8 m 
sec.!, which corresponds to an idealized velocity 
of 272.5 m sec.-! under these conditions, and 
agrees satisfactorily with the 9 kc measurements 
reported above. The frequencies at which these 
measurements were made were evidently well 
within the audible range, and the conclusion 
follows that the missing energy terms have 
already disappeared at frequencies even below 9 
ke at atmospheric pressure. Furthermore, Eucken 
and Miicke,? on the basis of heat capacity 
measurements at equilibrium, were led to the 
conclusion that the term corresponding to the 
symmetrical linear (O@C-0O) vibration was 
absent from the measurements on the velocity of 
sound at 3-4 kc made by Partington and 
Schilling and by Dixon. The missing term was 
believed to make a gradual reappearance only at 
high temperatures. On this assumption the ideal 
velocity at 9 kc should be 272.1 m sec.~! at 30°C, 
which agrees fairly well with the idealized 
measured value of 272.5 m sec.~'. In view of 
considerations later to be discussed the con- 
clusion of Eucken and Miicke will not be 
unconditionally supported, it being merely 
agreed that a term corresponding to about 1300 
wave numbers is lacking at 9 ke and ordinary 
temperatures. 


B. A revision of the customary notation 


In view of the fact that several dispersive 
regions, which occur in quite separate ranges of 
frequency, must be distinguished, a point of 
notation arises. In writing Vo for the velocity of 
sound corresponding to the equilibrium heat 


10 Buckendahl, Diss. Heidelberg (1906). This paper was 
unfortunately not available, and the frequency at which 
the measurements were made is, therefore, unknown to us. 

" Dixon, Campbell and Parker, Proc. Roy. Soc. A100, 
1 (1921). 
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capacity Co, and V., corresponding to a high 
frequency heat capacity C,, it has already been 
stressed that the subscripts do not indicate zero 
and infinite acoustical frequencies, but are 
intended only to specify whether wi? =0 or 
whether w3? = from the point of view of 
experiment. Since the values of the relaxation 
times of the excited states of the various modes of 
vibration may apparently be separated by many 
orders of magnitude, additional subscripts un- 
fortunately become necessary in order to 
avoid ambiguity. For example the transverse 
(OFC vibration may maintain sub- 
stantially its equilibrium distribution in the 
acoustical cycle at frequencies so high that the 
symmetrical linear (OeC-0O) vibration is 
wholly absent from the effective heat capacity. 
In order to avoid multiple subscripts the follow- 
ing convention has been adopted, the terms 

and “out” signifying equilibrium distri- 
bution and complete disappearance from the 
acoustical cycle respectively. 


Mode of vibration 
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These velocities designate a number of the 
possible terminal velocities of a colinear triatomic 
molecule, and may also be employed, by the 
usual conventions, to specify the corresponding 
quantities for any triatomic molecule. In this 
language the 9 ke velocity reported above 
corresponds approximately to both V; and to 
V; since the unsymmetrical linear (O3-—C-0) 
vibration is not appreciably excited at room 
temperatures. 


C. The variation of the velocity of sound with 
pressure and temperature in the dispersive 
region 

A detailed study was made of the variation of 
the velocity of sound with frequency and with 
pressure in the upper limit of the main vibrational 
dispersive region. This was undertaken because 
there is reason to believe that the 0-1 transition 
is somewhat less probable than the 1—2 trans- 
ition,’ and the transverse vibration, which 
predominates in this dispersive region, has an 
appreciable fraction of its molecules in the second 
excited state at 30°C. The results of this study 
are presented in Table I, together with the 
appropriate van der Waals factors cbhtained from 
the constants a= 3.61 X 10° cm® atm. mole and 
b=42.6 cc mole. The idealized measurements 
have been graphed in Fig. 1 by plotting V? in 

cm? sec.-?X10-* as ordinates against logio v/P, 

where v is the acoustical frequency in cycles per 


TABLE I. Measurements of the velocity of sound in carbon dioxide at 30°C and various pressures and frequencies. 


At 450.7 ke 


Pin mm Hg 70.5 
Vin m sec.~! 283.7 
w 1.0005 


776 
283.1 
1.0050 


1336 
282.2 
1.0087 


116 
283.6 
1.0007 


777 
282.9 
1.0050 


1640 
281.5 
1.0108 


237 
283.6 
1.0015 


778 
283.1 
1.0050 


70.5 
283.6 
1.0005 


800 
280.3 
1.0054 


116 
283.4 
1.0007 


880 
280.0 
1.0060 


237 408 
282.9 282.2 
1.0016 1.0027 


935 1012 
279.7 279.3 
1.0063 1.0068 


1.0026 


1.0052 


408 
283.4 


431 
283.5 
1.0028 


880 
282.8 
1.0057 


692 
283.2 
1.0039 


604 
283.4 
1.0039 


800 
282.8 


935 
282.9 
1.0060 


1012 
282.7 
1.0065 


"See for example, Richards, J. Chem. Phys. 1, 872 (1933). 
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At 94.0 ke 

if P 602 604 774 776 778 
V 281.3 281.5 280.4 280.2 280.0 

t w 1.0040 1.0040 1.0052 1.0052 1.0052 

s P 1180 1336 1474 1640 

h V 278.4 277.0 276.7 275.6 

w 1.0080 1.0091 1.0102 1.0114 
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Fic. 1. The variation of the velocity of sound in carbon 
dioxide with frequency and pressure at 30°C. Experimental 
points are indicated by plus signs at 451 kc, and by circles 
with lines through them at 94 kc, the extension of the 
vertical line being equivalent to +0.1 m sec.-!. The 
continuous line is drawn from theory, and the dotted line 
indicates the course of the experimental points at 451 kc if 
they are uncorrected for van der Waals forces. 


second and P the pressure in mm of mercury, as 
abscissae. The solid line in Fig. 1 is the theoretical 
V?, log v/P curve for an ideal gas in which 
V.,/V1=1.042, fio=fei etc. = 2.23 X10° per mole- 
cule per second at 760 mm, and in which the 
excitation is brought about wholly by double 
collisions. The dotted line shows the course of 
the 45 kc measurements if they are uncorrected 
for departures of carbon dioxide from ideal 
behavior. 

The idealized measured velocities average to 
give 283.7 m sec.~! as the terminal velocity. The 
ideal terminal velocity using an apparent heat 
capacity ratio of 1.400 is 283.2 m_ sec.—. 


Thus the ideal value of V.,/V; is 1.041 and the 
experimental value for the same quantity is 
1.043, which is considered a satisfactory agree- 
ment. It is to be noted that the ideal velocities 
are consistently lower than the idealized experi- 
mental values by about 0.4 m sec.~!. The reason 
for this is not understood, but it may be due, at 
least in part, to a discrepancy between the exact 
absolute units of length and of temperature and 
those used to obtain the experimental figures. 
The tube corrections permit intercomparison 
of measurements but do not, of course, eliminate 
absolute errors. 


TABLE II. The variation of the velocity of sound in carbon 
dioxide with temperature at 774 mm pressure, 
and at 9 and 92 ke. 


Yi 10 13.75 21.9 30.0 59.0 70.0 90.9 
Vmsec.'9ke 262.9 — — 271.8 279.2 286.9 295.0 
V msec.“! 92 ke 268.4 270.5 273.4 274.6 279.9 — — 

w 1.005 1.005 1.005 1.005 1.004 1.004 1.004 


Measurements were made on the velocity of 
sound at 770-776 mm between 10°C and 90°C at 
9 kc and between 10°C and 50°C at 92 ke. These 
are summarized in Table II, the van der Waals 
factors having been obtained from constants 
already given. The idealized measured velocities 
are plotted against temperature in Fig. 2. It 
is to be noted that the velocity at 30°C and 92 ke 
is 274.6 m sec. in Table II and at 94 kc and the 
same pressure is 280.2 according to Table I. 
Only a small part of this discrepancy is accounted 
for by the difference in frequency. This is 
unfortunate, for it indicates that the purity of 
the carbon dioxide used for the study of the 
temperature coefficient was not as great as that 
used for the study of the pressure coefficient. 
This is confirmed by the unsystematic character 
of the points in the dispersive region of the 
former study, as illustrated in Fig. 2. Calculation 
shows that ‘the probability of the 1-0 transi- 
tion has been increased by a factor of about 3. 
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Fic. 2. The variation of the velocity of sound in carbon 
dioxide with temperature and frequency at 773 mm 
pressure. The various continuous and dotted lines conform 
to the convention explained in Section IIB with the 
exception of the dotted line connecting the 92 ke points, 
which have been indicated by circles. 9 ke points are 
indicated by black dots. 


Whatever the impurity present, however, its 
effect on the 9 kc velocity was too slight to be 
experimentally manifest, for the 9 kc measure- 
ments could be reproduced quantitatively with- 


out difficulty. 


D. Conclusions from the study of the velocity of 
sound in carbon dioxide 

(1) The van der Waals forces follow the 
acoustical cycle quantitatively throughout the 
regions of pressure, temperature and frequency 
which were studied. This is based on the course 
of the dotted line in Fig. 1 and similar uncor- 
rected constructions which have not been given 
to avoid confusing the diagrams. 

(2) The mechanism by which excitation and 
de-excitation of the vibrational energy is effected 
is wholly that of double-collisions within the 
limit of error of experiment. This conclusion is 
based on the agreement of the 94 and 451 kc 
Measurements at low pressures, and on the 
conformity of both to the fragment of a sym- 
metrical S-curve drawn in Fig. 1. As a corollary it 
may be stated that the radiation life-time of the 
first excited state of the transverse vibration is 
greater than seconds.'* 


“An argument communicated to us privately by 
Professor E. U. Condon indicates that this lifetime is 
between 1 and 100 seconds. 
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(3) The velocity of sound at 9 kc and 774 mm 
pressure between 10°C and 90°C conforms 
closely to that calculated on the assumption that 
the symmetrical linear (O@C-O) vibration is 
absent from the heat capacity under these 
conditions. The reason for this conclusion has 
been discussed at some length above in con- 
nection with the 30°C velocity at 9 kc, at other 
temperatures it follows from a consideration of 
Fig. 2. 

(4) The velocity of sound at the highest 
frequencies and lowest pressures employed in this 
study corresponds to a heat capacity in which the 
rotational terms are fully active and are classi- 
cally excited. It is maintained constant for at least 
one order of magnitude of frequency or of 
pressure change. This statement is justified by 
measurements at large values of log v/P in Fig. 1. 
The magnitude of this terminal velocity corre- 
sponds to that measured by Kneser,'* Piele- 
and Eucken and Becker," although none 
of these authors have demonstrated its constancy 
over so wide a range of conditions. 

(5) The probability of the 1-0 transition of 
the transverse vibration is 2.23 x 10° per molecule 
per second at 760 mm and 30°C. This is reason- 
ably in accord with the somewhat less carefully 
defined values given by Pierce,'7 Kneser!* and 
Pielemeier.'* The chemical purity of the gas is so 
acutely important in determining the absolute 
value of this constant that exact agreement 
between different investigators is not to be 
expected at present. 

(6) The probability of the 2-1 transition is 
equal to this within a factor of about 5. This 
follows’ from the correspondence of the experi- 
mental points and the continuous line in Fig. 1. 
A slight tendency for the points to fall under the 
continuous line may perhaps be imagined at 
large values of log v/P, but it is certainly 
insufficient to justify quantitative conclusions. 

(7) The temperature coefficient of the dis- 
persion is adequately accounted for if it be 
assumed that the 1-0 transition of the trans- 
verse mode is effected only by molecules able to 


14 Kneser, Ann. d. Physik 11, 777 (1932); 12, 1015 (1932). 

' Pielemeier, Phys. Rev. 41, 833 (1932). 

‘6 Eucken and Becker, Zeits. f. physik. Chemie B20, 
467 (1933). 

17 Pierce, Proc. Am. Acad. Sci. 60, 271 (1925). 
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supply an activation energy of 6000 cal. mole. 
On this assumption the relative translational 
energy for the 0-1 transition must evidently 
exceed 8000 cal. mole~!. The probable error, of 
this determination is great, and an accuracy 
better than +2000 cal. mole“! may not be 
claimed. The slope of the Vs: curve against 
temperature is difficult to determine from the 
experimental points, and the presence of a trace 
of impurity may evidently also somewhat have 
obscured the result. 

(8) The activation energy of the 2-1 trans- 
ition can evidently not be far inferior to this 
figure, since no flattening of the S-curve in Fig. 1 
is observed at large values of log v/P. 

(9) The effective collision diameter is found to 
have about one-fifth of its value for ordinary 
kinetic processes from the data given above. 
This conclusion is highly qualitative since the 
experimental determination of the activation 
energy is not exact. 


III. MEASUREMENTS ON CARBON DISULFIDE 


The carbon disulfide was obtained from the 
best certified commercial source. It was frac- 
tionated through phosphorous pentoxide, only 
the middle third being retained for study. The 
apparatus was pumped out and washed several 
times with pure vapor before each determination. 
It is apparent from the scattering of points in the 
dispersive region that this method of preparation 
was not wholly adequate. Additional difficulty in 
obtaining accurate measurements in the dis- 
persive region was occasioned by the enormous 
acoustical absorption of carbon disulfide vapor. 
Absorption reduced the intensity of the electrical 
reaction to about 1 percent of its value after 5 
wave-lengths in the center of the dispersive 
region. This is a somewhat greater effect than is 
to be expected from the magnitude of the 
acoustical absorption coefficient, although so 
large a part of the heat capacity disappears from 
the acoustical cycle with a hundred-fold change 
in frequency that outstanding absorptive prop- 
erties must result. 

Table III summarizes the measurements made 
at 30°C and various frequencies and pressures; 
they are presented graphically in Fig. 3, after 
being corrected for deviations from perfect 
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Fic. 3. The variation of the velocity of sound in carbon 
disulfide with frequency and pressure at 30°C. 451 kc 
measurements are indicated by plus signs, 94 ke by 
circles with lines through them, 92 kc by circles and 9 kc 
by black dots. The continuous line has been drawn from 
theoretical considerations. The measurements are not suf- 
ficiently accurate, however, to define unequivocally the 
high frequency terminal velocity, and the upper part of 
the S-curve must therefore be considered tentative. 


behavior. The correction constants employed 
were a= 11.61 10° cc? atm. and b=76.8 cc per 
mole. The high frequency terminal velocity is in 
some doubt, since no oscillator above 451 ke was 
available, and it was found impossible to make 
reproducible measurements at pressures below 80 
mm. That indicated in the figure, which corre- 
sponds to an apparent heat capacity including 
rotational terms and terms due to the unsym- 
metrical linear (SCS) vibration in their 
equilibrium distributions, appears to accord well 


TABLE III. The velocity of sound in carbon disulfide at 30°C 
and various frequencies and pressures. 


451 ke 
200 200 128 128 128 = 8% 
212.0 212.7 213.9 214.0 213.9 214.1 
1.005 1.005 1.064 1.004 1.004 1.002 


94 ke 
P mm 321 200 128 # £128 
V msec. 201.0 205.1 203.6 206.4 
w 1.010 1.006 1.004 1.004 


92 ke 
P 309 =128 
V 201.8 208.5 
w 1.010 1.004 


P mm 215 
V m sec.! 213.6 
w 1.006 


9 ke 
314 309 128 
200.3 199.5 201.0 
1.010 1.010 1.004 
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TABLE IV. The velocity of sound in carbon disulfide at 


various temperatures. 


At 200 mm and 451 ke 


ie 20.0 20.0 30.0 30.0 40.0 40.0 
V m sec.*! 210.0 209.8 212.0 212.7 214.7 214.7 
w 1.006 1.006 1.005 1.005 1.005 1.005 

At 309-314 mm and 92 ke 
i. 25.0 26.0 27.5 30.0 34.0 37.6 50.0 
V m sec.~! 200.6 200.7 201.8 201.8 202.9 203.9 206.4 
w 1.012 1.012 1.011 1.010 1.010 1.009 1.008 


At 309-314 mm and 9 kc 
25.0 30.0 30.0 50.0 60.0 70.9 80.0 
V msec. 198.7 200.3 199.5 206.5 209.6 212.5 215.6 
w 1.012 1.010 1.010 1.009 1.008 1.008 1.007 


with the measurements. The velocities at 9 kc 
correspond well with those calculated from band- 
spectral analysis for equilibrium conditions.'® 

Table IV summarizes in a similar manner 
measurements made on the temperature coeffi- 
cient of the velocity under various conditions. 
Fig. 4 represents graphically the result of the 
measurements at 9 and 92 kc and 309-314 mm 
only. 

These data are the justification of the following 
conclusions. 

(1) The van der Waals forces follow the 
acoustical cycle without lag throughout the 
region experimentally investigated. This follows 
from the mutual discordance of measurements 
before the w correction has been applied. 

(2) The transverse (STCJST) and sym- 
metrical linear (S@C-S) vibration leave the 
acoustical cycle at the same frequencies of sound, 
and their characteristic adjustment rates are 
identical within limit of error. This follows from 
the conformity of the experimental points to the 
symmetrical S-shaped V? against log v/P curve 
which is illustrated by Fig. 3. 

(3) The mechanism by which the excitation of 
these states is effected is one of double collisions 
only over the range and within the error of 


'§ The fundamertal molecular frequencies were obtained 
from Villars, Chem. Rev. 9, 369 (1932). The empirical 
formula given in the Int. Crit. Tables for the velocity of 
sound does not agree either with our measurements or 
with the band-spectral results. Since the values are 
consistently too high, it may be assumed that a trace of 
some unknown impurity has been unintentionally included 
in the construction of the formula. Air, for example, owing 
to its smaller density and heat capacity could produce 
such an error if present even in very small quantities. 
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Fic. 4. The variation of the velocity of sound in carbon 
disulfide with temperature and frequency at 310 mm 
pressure. The conventions are the same as in Figs, 2 and 3. 


experiment. Here, again, conformity to the 
S-curve in Fig. 3 is sufficient proof. 

(4) The activation energy of collision for that 
process of excitation which dominates in the low 
frequency part of the dispersive region is about 
5400 cal. mole. The effective cross section 
appears to be of normal magnitude for this 
excitation process, although the error of the 
experiments is great enough to make this 
statement somewhat questionable. The activa- 
tion energy and cross section have been obtained, 
as usual, from the slope of the line connecting 92 
kc points in Fig. 4. 

(5) The activation energy of collision and 
effective cross section for the excitation process 
which determines the upper part of the dispersive 
region are identical with those dominating the 
lower part within limit of experimental error. 

There are two methods for justifying this 
statement. The measurements at 200 mm and 
451 kc quoted in Table IV are in accord with the 
92 kc measurements graphed in Fig. 4, but 
uncertainty concerning the high frequency termi- 
nal velocity somewhat obscures their significance. 
Also, the conformity of the 451 kc measurements 
plotted in Fig. 3 with the symmetrical S-curve 
may be used as a powerful supporting argument. 
A difference of 2000 cal. mole“ in the activation 
energy of collision between the processes de- 
termining the lower and higher parts of the 
dispersive region would result, with normal 
cross sections, in a difference of nearly thirty-fold 
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in rate. This difference would be manifest by a 
marked tendency of the S-curve to flatten in the 
high frequency region, and the measurements 
although inexact, are ample to demonstrate that 
no such effect is present. 

(6) The probabilities of the 1-0, 2-1, etc., 
transitions of the deformation vibration are 
evidently about 7.8 X 10° per molecule per second 
at 30°C and 310 mm pressure. This is signifi- 
cantly greater than the corresponding proba- 
bilities for carbon dioxide which are 9.1 X 10* per 
molecule per second under these conditions. 


IV. MEASUREMENTS ON SULFUR DIOXIDE 


Relatively little attention has here been de- 
voted to sulfur dioxide because the band spectral 
analysis appears to be controversial. The 
measurements reported below were, indeed, 
obtained merely during the process of defining 
the tube-correction. They are printed to permit a 
rough comparison with the absorption coefficients 
measured by Grossmann.'® 

The sulfur dioxide was obtained from tanks of 
the usual commercial product. It was dried by 
phosphorus pentoxide, but not further purified. 
The absolute values of the velocities are, there- 
fore, not unassailable, and a scattering of 
experimental points is to be expected. 


TABLE V. The velocity of sound in sulfur dioxide at 30°C 
and various frequencies and pressures. 


At 451 ke 
770 
222.6 
1.013 


At 94 ke 
770 
221.1 
1.013 


200 
224.9 
1.003 


Pinmm 770 
V in m sec.~! 222.4 
w 1.013 


200 
223.7 
1.003 


770 
V 221.9 
w 1.013 


Table V summarizes the measurements at 
30°C at 94 and at 451 kc. The van der Waals 
factor was obtained from the constants a=6.7 
cm® atm. mole, cm* mole. A 
measurement at 9 kc, 30°C, and 480 mm gave a 
velocity of sound of 222.3 m sec.~!, and requires 
multiplication by a factor of (1.008)? in order to 
give the corresponding idealized velocity. 


19 Grossmann, Ann. d. Physik 5, 681 (1932). 


T. RICHARDS AND J. A. REID 


The ideal velocity at very high frequencies of 
sound will be (1.3333 RT/M)! or 229.0 m sec.~! at 
30°C if no vibrations contribute to the heat 
capacity, and if the rotations are classically 
excited and follow the wave without lag. At very 
low frequencies the heat capacity ratio from band 
spectra is 1.269 if the analysis of Bailey, Cassie 
and Angus,” which gives the fundamental 
frequencies of 610, 1164 and 1369 cm™' is 
employed, and 1.263 using the more recent 
values of 510, 1150, and 1370 cem™ due to Chow?! 
and Smyth.” Thus the ideal velocity of sound at 
very low frequencies should lie between 223.5 
and 222.9 m sec.~! at 30°C. This accords with the 
idealized measured 9 kc velocity. 

In Fig. 5 the experimental data are sum- 
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Fic. 5. The variation of the velocity of sound in sulfur 
dioxide with pressure and frequency at 30°C. The two 
dotted lines at low values of V? represent the low frequency 
velocity of sound according to two different band-spectro- 
scopic analyses. The curved dotted line connecting 
experimental points is based on an unproved hypothesis 
which is tentatively advanced in the text. 


marized graphically in the usual way. Several 
conclusions are clearly indicated by this diagram. 

(1) The van der Waals forces follow the 
adiabatic cycle quantitatively at 30°C to fre- 
quencies of 450 kc and 400 mm pressure. This 
follows from a plot, which has not been repro- 


20 Bailey, Cassie and Angus, Proc. Roy. Soc. London 
A130, 142 (1930). 

21 Chow, Phys. Rev. 44, 638 (1933). 

22 Smyth, Phys. Rev. 44, 690 (1933). 
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duced, of log v/P against the non-idealized 
measured velocities. 

(2) The dispersive region must be due, at least 
in part, to sluggishness in the adjustment of the 
energy states characterized by the fundamental 
frequency 510 cm~'. Supposing that the two 
high frequency molecular vibrations disappear 
from the adiabatic cycle, leaving only the lowest 
in equilibrium, the ideal velocity of sound would 
be about 224.3 m sec.~! at very high frequencies, 
and this is significantly lower than the highest 
idealized measured value of 225.3 m sec.. 

(3) The position of the dispersive region in 
sulfur dioxide leads to the expectation that the 
absorption maximum will be in the neighborhood 
of log v/P =3.5 or about 2400 ke at atmospheric 
pressure. 

The last of these conclusions permits a 
comparison with Kneser’s prediction,”* from the 
absorption measurements of Grossmann, that the 
absorption maximum must occur at about 20°C 
(the temperature is not exactly specified) and 
atmospheric pressure at about 1000 kc. A 
combination of these two figures yields an 
activation energy of 15,000 cal. mole by the 
approximate formula” 


So large an activation energy of collision cannot 
seriously be considered, since it implies a 
fantastically large effective molecular cross 
section in view of the rapid rate at which 
transitions occur. It is, therefore, necessary to 
supply an explanation for the discrepancy. 

The argument given in the following paragraph 
is presented for this purpose, but it must be 
considered to be highly tentative in character. 

If the dispersive or absorptive region is due to 
the concomitant dropping out from the acoustical 
cycle of a group of energy states which have 
identical rates of adjustment, the effect is 
normally complete in a frequency range the 
upper limit of which is a hundred times the 
lower limit. If, for example, dispersion first 
becomes manifest at 10 ke the high frequency 


*3 Kneser, Ann. d. Physik 16, 337 (1933). 

4 This is made on the basis of Eq. (15) of page 740, 
Richards and Reid, J. Chem. Phys. 1 (1933). The approxi- 
mations are fo.=fi=---=0, (T/T’)'=1, and Bo/B,O 
= By /B,,'Q'. 
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terminal velocity will usually have been sensibly 
attained at 1000 kc at constant temperature and 
pressure. It will be noticed by consulting Fig. 5 
that although the velocity at log v/P=2.0 is 
sensibly greater than that at low frequencies, the 
dispersive region will not be complete at log 
v/P=4.0. Furthermore, a function of the type 
of the dotted line appears better suited to de- 
scribe the experimental points than the con- 
ventional S-shaped curve. This suggests that one 
of the series of energy states corresponding 
either to 1160 cm™ or to 1370 cm™ is somewhat 
less readily adjusted by collisions than those 
states corresponding to the 510 cm™ vibration. 
This is supported by the slight fall of the 
absorption coefficient measured by Grossmann at 
the beginning of the vibrational absorptive 
region. If this hypothesis is granted, the ab- 
sorption maximum predicted by Kneser must be 
shifted to higher frequencies at 20°C and 
atmospheric pressure, and the magnitude of the 
activation energy of collision which must be 
assumed is very greatly diminished. Neither the 
velocity nor the absorption measurements are of 
sufficient accuracy to define the matter more 
exactly, but there appears at present no con- 
vincing ground for believing them incompatible. 

Since the terminal velocity at high frequencies 
may be measured at 451 kc and 30°C only at 
pressures below 70 mm no quantitative statement 
of the transition probabilities of the various 
vibrational energy states may here be made. The 
predominant adjustment rate at 30°C is evidently 
about five times that of carbon disulfide, and 
about 15 times that of carbon dioxide. 


V. Discussion 


It is apparent that at least one dispersive 
region in the range of acoustical frequencies now 
experimentally available is to be expected with 
polyatomic gases in which the vibrational heat 
capacity is appreciable. It is of little interest, 
therefore, merely to multiply the number of 
examples of dispersive gases, and attention is 
rather to be directed to the information which 
they provide concerning the mechanism of the 
transformation of kinetic to vibrational energy. 
Unfortunately, however, attempts to correlate 
experiment and theory must at present be highly 
speculative, since the descriptions of these 
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phenomena in terms of quantum mechanical 
processes are so general that they do not lend 
themselves to the treatment of specific physical 
examples. In view of the predominantly experi- 
mental character of this study no theoretical 
discussion will here be attempted. The following 
paragraphs have been included to point out that 
our measurements are not at present unequivocal 
in their significance, and to indicate directions for 
further investigation. 

Carbon dioxide and carbon disulfide so closely 
resemble each other in structure that it is 
surprising to find an appreciable part of the heat 
capacity failing to participate in the velocity 
of sound at audible frequencies in the former, 
whereas in the latter all of the heat capacity 
remains fully active well into the ultrasonic 
range. It is natural, in considering carbon 
dioxide alone, to follow Eucken and Miicke in 
supposing that the symmetrical linear vibration 
is that with the extremely low rate of adjustment. 
In searching for a kinetic mechanism by which 
such a mode of vibration may be excited a 
fortunate triple-collision of the type A~OCOA 
suggests itself; any bond-softening process or 
other double-collision mechanism seems ineffi- 
cient in comparison. This supplies an excellent 
qualitative explanation for the behavior of carbon 
dioxide, but fails completely to account for that 
of carbon disulfide. 

An examination of the quantitative differences 
in the structure of the energy levels of the two 
molecules serves only to make their difference in 
behavior more puzzling. It is probable, for 
example, that the various energy states in the 
symmetrical linear vibration and the deformation 
vibration ‘‘communicate” during the time that 
two colliding carbon disulfide molecules are 
effectively in contact. The equality of the rates of 
adjustment of these states, which has experi- 
mentally been established, can reasonably be 
explained only on the supposition that vibrational 
equilibrium between them has essentially been 
reached on every effective collision. Under these 
conditions all the ‘‘communicating”’ vibrations 
will, of course, fall out of the acoustical cycle 
with the rate of adjustment of the most readily 
excited member of the group, and the steric and 
energetic factors which are measured characterize 
the most probable transition. (This is sharply in 


contrast to the behavior of nitrous oxide,” and 
may also differ from that of sulfur dioxide, where 
it appears that each of the vibrational modes is 
dropping out with a separate characteristic rate. ) 
It is further apparent that the probability of 
““communication”’ between two vibrational modes 
increases greatly as the number of excited states 
in each is increased, since the matrix component 
connecting two widely separated energy levels is 
necessarily small. Since carbon disulfide contains 
a greater number of excited vibrational states at 
any given temperature than carbon dioxide, and 
since the duration of a collision between the 
heavier molecules is also probably greater, it 
would be natural to attribute their difference in 
acoustical behavior to the effect of “‘communi- 
cation” alone. A _ serious difficulty to this 
interpretation results, however, from the pertur- 
bation of the first state of the linear vibration by 
one of the three levels corresponding to the 
second state of the transverse vibration. Per- 
turbation indicates that one of these states may 
transform itself into the other during the course 
of a few vibrations without the intervention of a 
collision. It is, therefore, difficult to believe that 
the linear vibration alone drops out of the 
acoustical cycle at low frequencies while the 
deformation states remain in equilibrium well 
into the ultrasonic range. 

It would be convenient, if it were possible, to 
discredit entirely the measurements at 9 ke on 
carbon dioxide, which are here reported, and 
obliterate the difference in behavior on a priori 
grounds. Unfortunately, however, they corrobo- 
rate a considerable body of independent evidence, 
and only the case of carbon dioxide is difficult to 
interpret. It is necessary, therefore, to confess 
that the low frequency dispersive region in 
carbon dioxide is not understood by the writers, 
and to leave the matter in this highly unsatis- 
factory condition. 

The situation is hardly more satisfactory when 
an explanation of the so-called activation energy 
of collision is sought. Three alternative hypothe- 
ses are at hand, and none of these may be 
eliminated with confidence. 

(1) A bond-softening process has been sug- 
gested as the mechanism by which the excited 
state is reached. In this case the collision energy 


25 Kneser and Ziihlke, Zeits. f. Physik 77, 649 (1932). 
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must necessarily be greater than the energy 
available for the bond-softening. process. This 
implies that the activation energy shall decrease 
roughly in proportion to the ratio between the 
energy of the state to be excited and the zero 
point energy. The experimental evidence fails to 
show this effect either with carbon dioxide or 
with carbon disulfide, since in neither gas is the 
symmetry of the S-shaped V*, log »/P curve 
appreciably disturbed. If, however, equilibrium 
between the vibrational states may be established 
during the time of a collision, as has been sup- 
posed above, this difference cannot be detected 
by experiment. 

(2) Heil has suggested that, because of the 
comparatively small rigidity of the bonds 
attaching the oxygen atoms to the carbon, a 
broadside collision of the type >--- : — would 
result in the transference of only about a third of 
the relative translational energy of the molecules 
into the vibration if the requirements of the law 
of conservation of linear momentum are to be 
satisfied. On this basis he has_ successfully 
predicted the order of magnitude of the activa- 
tion energy of collision and the effective cross 
section of carbon dioxide from Kneser’s measure- 
ments at atmospheric temperatures. A similar 
argument applies equally well to carbon disulfide: 
the available energy is smaller, but the rigidity of 
the bond is also less, and an activation energy of 
about the same magnitude is to be expected. The 
situation is somewhat complicated, however, by 
the rotation of the molecules during a collision. 
The activation energy represents, of course, 
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predominantly the minimum potential energy 
which must be available in order to de-excite a 
given state, and there appear to be several 
configurations of rotating molecules more favor- 
able for the process than that which Heil has 
described. An explanation based on a dynamic 
mechanism of this kind has, however, the 
advantage that it implies approximately equal 
activation energies for the various states of a 
vibrational mode, and this is in accord with 
experiment. 

(3) It is also possible that a molecule which is 
vibrating has an advantage over a molecule 
which contains only its zero point energy in 
bringing about transformations of the relative 
kinetic energy into vibrational energy by col- 
lision. If, for example, it is supposed that only the 
states 0 and 1 are present, and that a molecule in 
the state 1 is far more effective per collision in 
causing the transitions 1-0 and 0-1 than a 
molecule in the state 0, the transition probability 
fio will not be proportional to the total number of 
collisions but to the number of 1-1 collisions, 
The relation between the number of molecules in 
the 0 and in the 1 states is of course m; = mge~"o/** 
and consequently, for gases like carbon dioxide 
and carbon disulfide at ordinary temperatures, 
will increase exponentially with temperature, and 
an apparent activation energy will be observed 
which has no relation to the activation energy 
ordinarily discussed in chemical kinetics. This 
possibility may also seem remote, but it is 
difficult to devise an experiment by which it may 
convincingly be eliminated. 
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the Lower Vibrational States of Ethylene, Together with Some Observations 


The velocity of sound in pure ethylene has been deter- 
mined at various pressures, temperatures and frequencies. 
It is concluded that equilibrium is attained between the 
various states of the vibrational energy on every effective 
collision, but that a considerable activation energy of 
collision is required for the conversion of translational into 
vibrational energy. Collisions with argon, helium, and 
nitrogen molecules have no appreciable effect on the 
vibrational energy of the ethylene molecule. Collisions 
with hydrogen molecules are about ten times as effective 


Concerning the Excitation of Rotational Energy in Hydrogen 


T. RicHARDS AND JAMES A. Princeton University 
(Received February 7, 1934) 


as ethylene-ethylene collisions in producing transitions in 


the lower vibrational energy states of ethylene. Hydrogen 
shows a dispersive region at ultrasonic frequencies which 
is apparently due to the failure of the rotational energy 
transitions to follow the acoustical cycle. It has been 
found necessary to suppose that ethylene-hydrogen colli- 
sions are about twenty times as effective as hydrogen- 
hydrogen collisions in exciting the rotational energy of 
hydrogen molecules. 


INTRODUCTION 


NE of the outstanding unsolved problems in 
chemical kinetics is the mechanism of 
activation by collision. It has become increasingly 
apparent that great differences exist in the 
relative activation efficiencies of various types of 
molecules, but the theoretical approach to the 
problem is difficult, and it is apparently impos- 
sible at present to make quantitative predictions. 
The multiplication of suitably chosen experi- 
ments becomes therefore highly desirable. 


Measurements on the dispersion and ab-- 


sorption of sound provide ready access to data 
concerning a somewhat analogous process: the 
excitation by collision of the lower vibrational 
energy states. This differs from the chemical 
activation of the molecule in several respects, but 
is even harder to describe from a theoretical 
standpoint because the interaction of the mole- 
cules must here be attacked by quantum 
mechanics, whereas there is every indication that 
in chemical activation the problem will ulti- 
mately yield to a classical dynamical treatment. 
Any final description of the mechanism by which 
translational energy is converted into vibrational 
energy must, however, be sufficiently general to 
embrace transitions between both the lower and 
higher vibrational states. A study of the collision 
efficiencies of various molecules in adjusting the 
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heat capacity cannot fail, therefore, to have 
ultimate bearing upon the chemical problem. 

The first research of this character was 
undertaken by Abello,! who measured the 
absorption of ultrasonic energy in mixtures of 
carbon dioxide, hydrogen and nitrogen. The 
theoretical equipment necessary to interpret his 
results was provided shortly thereafter by 
Bourgin.? In 1932 Rogers* reported briefly con- 
cerning the effect of water vapor on the ab- 
sorption of sound in carbon dioxide, and an 
advance notice of the present work was pub- 
lished. During the past year Kneser has in- 
terpreted existing measurements on air and on 
oxygen-water-vapor mixtures,’ and Eucken and 
Becker® have published conclusions from a few 
measurements on the velocity of sound in binary 
mixtures containing as one constituent either 
carbon dioxide or chlorine. This appears to 
complete the meager list of measurements which 
have been published to date. 


1 Abello, Proc. Nat. Acad. Sci. 13, 699 (1927). 

* Bourgin, Nature 122, 133 (1928); Phil. Mag. 7, 821 
(1929); Phys. Rev. 34, 521 (1929). 

3 Rogers, Phys. Rev. 41, 369 (1932). 

4 Richards and Reid, Nature 130, 739 (1932). 

5 Kneser, Ann. d. Physik 16, 337 (1933); J. Acous. Soc. 
Am. 5, 122 (1933). 

6 Eucken and Becker, Zeits. f. physik. Chemie B20, 467 


(1933). 
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I. PurE ETHYLENE 


A. The purity of the ethylene gas employed 

The ethylene gas was obtained in compressed 
cylinders, which are supplied by the U. S. 
Industrial Chemical Company, and rated as 
about 99 percent pure. Four of these cylinders 
were employed, and of these three gave fairly 
reproducible results. The fourth gave values 
which varied from day to day, and all measure- 
ments on gas from this cylinder have been 
discarded. The cylinders were cooled to — 20°C to 
liquefy the ethylene; and then “blown off’’ to 
remove any traces of non-condensing gases. 
Ethylene was then evaporated through phos- 
phorus pentoxide into evacuated acoustical 
chambers of a design which has already been 
described.’ The apparatus was washed with gas 
several times before each measurement. 

The measurements, combined by means of 
smoothed curves, permit comparison for theo- 
retical purposes with an accuracy somewhat 
greater than 0.1 percent. It is not possible, 
however, to estimate here their absolute accu- 
racy. 


B. The low frequency terminal velocity 


A measurement of the velocity of sound at 9 
ke, 30°C, and 770 mm gave 332.8 m sec.~'. In 
order to idealize this quantity it was multiplied 
by a factor of 1.0035, which was obtained from 
the van der Waals constants a=4.47 10° cm® 
atm. mole! and b=57.14 cm* mole. The 
resulting figure is 333.9 m sec.~!, which disagrees 
gravely with 336.1 m sec. calculated in a 
similar manner from the empirical equation 
given in the International Critical Tables. We 
cannot account for this discrepancy. It should be 
noted, however, that the equation given in the 
International Critical Tables results from a combi- 
nation of the measurements of six different 
observers, who are in violent mutual disagree- 
ment. 

No wholly satisfactory data on the heat 
capacity appear to be available. The most recent 
study has been made by Haas and Stegeman,® 
who give values which are slightly lower than 
those reported by previous writers. Their results 


Richards and Reid, J. Chem. Phys. 1, 114 (1933). 
* Haas and Stegeman, J. Phys. Chem. 36, 2127 (1932). 
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are quoted without specification of the pressure, 
and are considered to vary linearly with tempera- 
ture from 3° to 67°C; both of these approxi- 
mations leave much to be desired from the 
standpoint of the present work. If their empirical 
equation for the temperature coefficient is valid 
(1/C,)(dC,/dT) =2.5X10-*, and if it is con- 
sidered that the empirical equation refers to the 
heat capacity of the gas at zero pressure Cy at 
30°C is 8.20 cal. mole“! deg.~!. The corresponding 
ideal velocity is 334.2 m sec.~', whereas the 
idealized experimental velocity at 9 ke is 333.9 
m sec.~'. The agreement is better than that to be 
expected considering the various possibilities for 
error, and it is concluded that no important part 
of the heat capacity of ethylene is lacking in the 
adiabatic compressibility at 9 kc, 760 mm 
pressure and 30°C. No band spectral analysis of 
sufficient scope to confirm or to modify this 
conclusion is known to us. 


C. The variation of the velocity of sound with 
pressure and frequency 


The results of a study of the variation of the 
velocity of sound in the dispersive region with 
pressure at 30°C are reported in Table I. A large 
number of determinations were made at 790 mm 
because this measurement served as a test of the 
purity of the ethylene. The spread of points is 
great, being +0.45 m sec.~! for the 94 kc and 
+0.75 m sec. for the 451 ke points. It is to be 
noted, however, that the difference V4;;— V9, has 
a maximum variation of +0.35 m sec.~, which is 
less than half that of the 451 kc measurements. 


TABLE I. The variation of the velocity of sound in ethylene 
with pressure and frequency at 30°C. 


At 451 kc 
472 236 83 60 
337.6 341.4 345.7 
1.004 1.002 1.001 


At 94 ke 

472 236 83 
333.9 334.5 339.5 
1.004 1.002 1.001 


Pinmm 790 
Vinmsec.-! 335.4* 
w 1.007 


60 
344.8 345.7 
1.001 1.001 


790 
V 333.1° 
w 1.007 


* These figures are the average of the following determi- 
nations each pair of which was made on a separate sample 
of gas. 


450 ke 
94 ke 


335.9 335.6 
333.4 


334.7 335.4 
332.9 332.7 


336.2 
333.3 


334.9 


335.8 
33 332.7 


333.6 333.1 


335.3 
333.1 
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1.200 


9 


Idealized 


0 2.0 3.0 40 
log,(v/P) 


Fic. 1. The variation of the velocity of sound in ethylene 
with frequency and pressure at 30°C. 9 kc measurements 
indicated by black dots, 94 ke by circles, and 451 kc by 
plus signs. The line connecting the points is drawn from 
theoretical considerations discussed in the text. 


Fig. 1 summarizes graphically the results of 
Table I by plotting, in the usual way, the square 
of the idealized velocity against the logarithm of 
the ratio of the acoustical frequency in Herzian 
cycles and the pressure in mm Hg equivalents. 
Several points in connection with this diagram 
should be noted. The S-curve has been drawn, 
according to the customary three-state approxi- 
mation, from the points corresponding to Co 
=8.20 cal. deg.' mole and C,=5.96 cal. 
deg.-! mole“. The experimental values in the 
neighborhood of the supposed high frequency 
terminal velocity were made at very low pres- 
sures, and, because of the small mass reaction 
of the gas on the oscillator under these conditions, 
were of conspicuously less accuracy than the 
other determinations. An oscillator of higher 
frequency was not available, and no search for a 
constant velocity in the upper part of the 
dispersive region could therefore be undertaken. 
The moment of inertia of the hydrogen atoms 
around the carbon-carbon axis is so small that 
ethylene may enter a rotational dispersive region, 
analogous to that exhibited by hydrogen,° before 


® Abello is apparently the discoverer of this dispersive 
region. Bourgin and Herzfeld and Rice (Phys. Rev. 31, 


the dispersive region due to the vibrations is 
wholly complete. The upper part of the S-curve 
has been drawn as a broken line to emphasize 
this uncertainty. No evidence for complexity of 
the dispersive region is, however, revealed by the 
measurements, since they fall uniformly on the 
S-curve within limit of experimental error. 


D. The variation of the velocity of sound with 
temperature 


A study of the temperature coefficient of the 
velocity at 94 and at 451 ke gave a series of 
measurements which are reported in Table II. 
The temperature coefficient was studied with 
somewhat more than the usual care, since by 
means of it certain conclusions may be reached 
which are directly serviceable in Section IT. 


TABLE II. The velocity of sound in ethylene at 790 mm 
pressure and various temperatures. 


At 451 ke 
150 15.3 225 250 WMO 3235 
V m sec. 328.3 328.7 331.8 333.5 334.2 335.4* 3306.3 
i 35.0 35.0 37.5 40.0 40.0 42.0 45.0 
V 336.8 337.1 338.6 340.1 340.0 340.5 341.7 
At 94 ke 


15.0 15.3 20.0 22.5 25.0 27.8 30.0 
325.5 325.7 328.1 329.4 330.7 332.0 333.1" 


35.0 35.0 40.0 40.0 42.0 45.0 45.0 
335.2 335.5 337.4 337.5 338.5 340.1 340.0 


* Each of these values is the average of the eight deter- 
minations on different samples of ethylene gas which are 
quoted in Table I. 


The continuous line connecting the 94 kc 
measurements is drawn from the data for the 
heat capacity quoted in B of this section, and 
from the equation Vo= {(RT/M)(1+R/C,)}?. It 
is to be noted that the disagreement of the 
points with the theoretical line is never greater 
than 0.3 m sec.~!, which is well within a reason- 
able limit of absolute error for the velocity 
determinations. The theoretical line is, however. 
almost rectilinear, whereas a slight consistent 
tendency for the points to fall under this is 
observable both below and above 30°C. A dotted 
line has been drawn between the 451 kc measure- 


691 (1928)) have attributed it to the failure of rotational 
energy adjustment; Kneser has recently supported this 
interpretation. 
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ments to give the most probable value of the 
velocity at this frequency at 790 mm pressure. 
The deviations of the experimental points from 
this average value is naturally somewhat greater 
than in the case of the 94 kc measurements. 


E. Conclusions from the velocity of sound in 
pure ethylene 

Owing to the lack of adequate band-spectral 
data it is not possible at present to describe the 
velocity of sound in the dispersive region in 
fundamental terms. The results may, however, be 
expressed in terms of numbers which will take on 
an exact meaning when the energy levels of the 
molecule are unequivocally defined. 

From the ideal velocity of sound at a given 
temperature and various pressures and_ fre- 
quencies a characteristic time, 3, may be com- 
puted by the relationship 


I= (Co/wC.,) (Vu? — Vo?) /(V..2— (1) 


in which Cp and C.,, are the heat capacities when 
=0 and =a, respectively, Vo and V,, the 
corresponding ideal velocities, and V. the 
velocity at a circular frequency w. This expression 
is approximate in several respects: notably a 
correction for the absorption has been neglected, 
and the time # is without quantitative meaning in 
terms of specific transitions. The variations of 3 
with pressure and temperature, nevertheless, 
permit certain conclusions. 

The conformity of experimental points to a 
symmetrical S-curve illustrated in Fig. 1 demon- 
strates, for example, that double collisions are the 
only appreciable means by which the various 
energy states which cause the dispersive region 
are excited. Similarly, since 3 has been found to 
be independent of frequency and directly de- 
pendent upon the numerical density, the trans- 
ition probabilities are sensibly equal from an 
acoustical standpoint for all of the contributing 
states. This strongly suggests that the vibrational 
energy attains equilibrium on every effective 
collision. The maximum acoustical absorption at 
760 mm and 30°C evidently occurs at about 
1.0 10° cycles sec.—!. At this temperature and 
pressure the time #@ has a value of 2.38 X10 sec. 

The activation energy of collision is evidently 
not great since the increase in V4;;— V4 is only 
1.2 m sec. between 45° and 15°C (Fig. 2). In 
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Fic. 2. The variation of the velocity of sound in ethylene 
with temperature at 780 mm. The line connecting the 94 kc 
points, which are indicated by circles, is the calculated ideal 
velocity. The dotted line connecting the 451 ke points, 
which are indicated by plus signs, assumes an activation 
energy of collision of 2800 cal. mole“. 


such cases it is probable that the variation of the 
effective molecular cross section with temperature 
must be taken into account in the kinetic 
analysis of the transition probabilities. However, 
an upper limit for the activation energy of 
collision A may be simply obtained by writing, at 
temperatures T and 7”, 


which involves the approximations that the 
cross section of the molecule is independent of the 
temperature, and that 
where x=hyo/kT. Since & is 1.93X10-" sec. at 
45°C and 2.73 X10~7 sec. at 15°C, the activation 
energy of collision is evidently about 2800 cal. 
mole. If this value is substituted into the usual 
kinetic expression, it is found that the least 
separation of ethylene molecules on an effective 
collision is about 2X 10-$ cm, which is only about 
a half of that found in ordinary kinetic processes. 
Since the measurements of the temperature 
coefficient of the velocity are unusually precise 
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and concordant, it is believed that the agreement 
of these two quantities is significant, and that no 
serious stearic difficulty can be responsible for the 
failure of the vibrational energy to maintain its 
equilibrium distribution in the sound wave at 
high frequencies. A line drawn on the suppo- 
sition that the activation energy of collision is 


v= 
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2800 cal. mole is identically the dotted line of 
Fig. 2. 
Il. Binary MIXTURES OF ETHYLENE WITII 


OTHER GASES 


The velocity of sound in a binary mixture of a 
gas A and a gas B will be approximately'’ 


pam+ pp(1—m) 


mixture, P, C and p refer to pressure, heat 
capacity, and density at P respectively, and m 
is the concentration of A expressed as a mole 
fraction. This expression is approximate in that it 
neglects absorption, and in that the charac- 
teristic time represents only roughly the reci- 
procal of the sum of the various transition 
probabilities involved in that part of the heat 
capacity (Cos — Coa) which causes the dispersive 
region. 

The velocity of sound in binary mixtures of 
ethylene with nitrogen, argon, helium and 
hydrogen has been measured, and the results and 
conclusions are reported below. The two gases 
were mixed under pressure in a reservoir and 
subsequently fed into the acoustical apparatus at 
a pressure of 780 mm. The gases were dried with 
phosphorus pentoxide, but not otherwise purified. 
The mixture was analyzed at the end of an 
experiment by absorbing the ethylene in fuming 
sulfuric acid, and measuring the residual volume 
of gas. This procedure was not of great accuracy, 
and a reproducibility of the measurements com- 
parable to that obtained in the previous section 
is not to be expected. The temperature was 
uniformly maintained at 30°C. 

Since the P, v, T relations of these mixtures 
have not been extensively studied, the procedure 
of idealizing the measured velocities of sound by 
the use of an acoustical equation which includes 
the second virial coefficient is not practicable. As 
an alternative the actual properties of the 
various gases, from data given in the International 
Critical Tables, were used for substitution into 
(3). They are summarized in Table III. 


{1 R Coam+Con(1 —m) +0 ] 


[Coam+Con(1—m) P+ 0??? 


where # is the characteristic time in the particular 


(3) 


TABLE III. The acoustical constants of several gases. 


Gas A CoH, He 


ping 1.649 1.166 1.651 8.407 4.357 
at 10-3 10-3 10-4 «105 x10" 
30°C and 
780 mm 


Co/R at 1.5 4.115 1.5 2.46 2.504 
30°C 


A. Ethylene-argon, nitrogen, and helium mix- 
tures 

The various measurements at 94 and 451 ke 
are summarized in Table IV. In the case of 
nitrogen the results may be numerically repro- 
duced by (3) only if 3 is taken as being identical 
with J obtained from the study of pure ethylene 
at a corresponding partial pressure. The agree- 
ment between the calculated and observed 
velocities is within experimental error, as a 


TABLE IV. The velocity of sound at 780 mm and 30°C in 
various binary gaseous mixtures containing 
ethylene as one constituent. 


Argon mixtures 


m 0.952 0.850 0.625 0.445 
Vas1 333.2 330.1 324.3 322.0 
Vos 331.3 327.0 319.9 316.8 


Nitrogen mixtures 
m 0.970 0.959 0.895 0.870 0.851 0.800 0.635 0.500 
Vas. 335.7 335.8 337.5 337.7 337.7 339.2 342.6 346.2 
Vous 333.8 334.1 335.2 335.3 335.8 343.1 339.5 — 


Helium mixtures 


m 0.940 0.900 0.852 0.760 
343.5 348.6 356.4 343.4 
Vo — 345.0 352.4 368.5 


- 10 This type of expression was first developed by Bourgin. 
The one here printed has been obtained from Einstein's 
treatment of the dispersion theory. See Richards, J. 
Chem. Phys. 1, 877 (1933). 
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V (m.sec-') 


332 
%CoH, 100 75 50 
N, 0 25 50 


Fic. 3. The velocity of sound in nitrogen-ethylene 
mixtures at 780 mm and 30°C. The lines are calculated 
from the data in Table III; black dots indicate 94 kc 
and circles 451 ke. 


glance at Fig. 3 demonstrates. With argon the 
agreement is not so satisfactory, since the 
experimental points lie increasingly above the 
calculated line with increasing concentrations of 
argon. A measurement on a s2mple of the argon 
employed gave a velocity of sound of 325.3 
m sec.-! both at 451 and at 94 kc. Since this is 
conspicuously higher than the value 324.1 m 
sec.-' obtained from (3) when m=O the dis- 
crepancy may safely be attributed to impurity of 
the argon. 

With the helium measurements it is apparent 
that an even greater error due to impurity is 
present. At m=0.75, V451 calculated using the 3 
obtained from pure ethylene, is 386.7 m sec.—! and 
Vs51 obs.y read from a smoothed curve through the 
experimental points, is 375.0 m sec.~'. Similarly, 
at m=0.75, Vos cate. 382.5 whereas ons. is 
only 370.3. It is to be noted, however, that 
(Visi— Voa)eate. iS 4.2 m sec.-' which agrees 
reasonably well with (V451— Vo1)ons. which is 4.7 
m sec.“'. Similarly at m=0.90 the calculated 
difference is 3.0 m sec.-! and the experimental 
3.6 m sec.-'. There can be little question, 
therefore, that the cause of error is the inclusion 
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Fic. 4. A comparison of the measured and calculated 
velocities of sound in ethylene-argon mixtures at 780 mm 
and 30°C, 


in the helium of a few percent of some denser 
impurity. Since a priori this is extremely likely, 
no serious concern is felt because of the lack of 
absolute agreement between the experimental 
and theoretical velocities. 

Within experimental error, therefore, collisions 
of argon, nitrogen, and helium molecules fail 
entirely to bring about vibrational energy 
transitions in ethylene. It may roughly be 
estimated that their collision efficiencies must be 
less than a tenth that of the ethylene-ethylene 
collision. 


B. Mixtures of hydrogen and ethylene 


Measurements on hydrogen-ethylene mixtures 
are reported in Table V. It should be noted that 
the determinations at 94 kc which result in a 
velocity of sound greater than 400 m sec.~ are of 


TABLE V. The velocity of sound in hydrogen-ethylene mixtures 
at 780 mm and 30°C. 


0.306 0 
1409 
1319 


m 0.950 0.940 
Van. 342.5 345.3 354.0 
341.6 344.7 353.1 


0.890 0.780 0.51 
376.4 462.9 584.2 
376.1 461.5 578.8 


212 WwW. 


considerably less accuracy than those below this 
figure. This is both because the tube-correction 
has not been determined for the apparatus at 
large velocities, and because the dimensions of 
the apparatus did not permit the detection of 
many points of maximum reaction when the 
wave-length was large. The agreement of the 
velocity at very low frequency calculated from 
(3) and that determined by experiment at 94 kc 
is in all cases satisfactory, however, and little 
hesitation is felt in basing approximate con- 
clusions on the entire range of measurements. 
It is at once apparent from an inspection of 
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Table V that two independent tendencies are at 
work, since A V(= V451— Vos) nearly vanishes in a 
mixture for which m=0.75, and _ increases 
markedly on either side of this composition. 
This behavior is plausibly explained if it is 
supposed that hydrogen-ethylene collisions are 
more efficient than ethylene-ethylene collisions in 
exciting the vibrational energy of ethylene, and 
that, conversely, ethylene-hydrogen collisions are 
more efficient than hydrogen-hydrogen collisions 
in exciting the rotational energy of hydrogen. A 
somewhat complicated situation results, which 
may most simply be described by the expression'! 


P Ria —wd 408) (Co— 49 + Rw? + 


p (Co— CI +0LC,,(84 +08) 


In (4) subscript A refers to ethylene and sub- 
script B refers to hydrogen, and 


Co= Coam+Con(1—m), 
C,am+C,,n(1 —m), 


the heat capacities being expressed per mole and 
m indicating the mole fraction of A. As before 
pam+pr(1—m)=p, pa and pz being, of course, 
the densities of A and of B at P. Finally 


vA —m) |, 
OB —m) ], 


where eg is the relative efficiency, at unit 
concentration, of the ethylene-hydrogen collision 
and the ethylene-ethylene collision, and €, is the 
relative efficiency, at unit concentration, of the 
ethylene-hydrogen and the hydrogen-hydrogen 
collision. #44 is the “relaxation time’’ of the 
ethylene-ethylene collision, and 3°, that of the 
hydrogen-hydrogen collision, both expressed at a 
numerical density corresponding to P. These are 
determined from the properties of the pure gases, 
and appropriate values of 34 and 0g found from 
gaseous mixtures containing small percentages of 
A and of B respectively. 

In order to relate 34 and #, to the relative 
efficiency of A—B and A—A or B—B collisions, 


This corresponds to Eq. (6.1) Richards, J. Chem. 
Phys. 1, 877 (1933). 


expressed per collision, it is further necessary to 
consider the velocity and cross section of the A 
and B molecules. Ordinary kinetic considerations 
show that, at equal numerical densities, the 
factor relating number of collisions to concen- 
where o indicates the kinetic diameter and J/ the 
molecular weight. By using o4 (ethylene) =4.3 


7.5 


ie) 
%CaH4 100 7S 50 25 
2s 50 75 


Fic. 5. A comparison of the measured difference in 
velocity of sound at 451 and at 94 ke. (AV) with that 
calculated on the assumption that an ethylene-hydrogen 
collision is twenty times as effective in causing transitions 
in the rotational energy of hydrogen as a_ hydrogen- 
hydrogen collision, and ten times as effective in causing 
transitions in the vibrational energy of ethylene as an 
ethylene-ethylene collision, 
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«10-8 cm and ox (hydrogen) it is 
found that a molecule suffers 1.53: times as many 
collisions in a given concentration of hydrogen as 
in a given concentration of ethylene. 

The results of the measurements at 780 mm 
and 30°C indicate that a hydrogen-ethylene 
collision is about 10 times as effective as an 
ethylene-ethylene collision in exciting the vi- 
brational energy of ethylene, and that an 
ethylene-hydrogen collision is about 20 times as 
effective in giving rotational energy to hydrogen 
as a hydrogen-hydrogen collision. The corre- 
spondence of these assumptions with experiment 
is indicated by Fig. 5, where the heavy line 
indicates AV calculated from (4) on this basis and 
the circles are the experimental points. The 
agreement is within experimental error, and 
justifies the assumption, which is implicit in (4), 
that the adjustment of rotational energy in 
hydrogen is effected by a double-collision mecha- 
nism. It remains to attempt to reach an under- 
standing of these findings. 


III. Discussion 


It has been concluded in the previous sections 
that collisions with hydrogen are extremely 
effective in bringing about transitions between 
the lower vibrational energy states of ethylene, 
whereas collisions with argon, helium and 
nitrogen are virtually without effect. This raises 
the question which has been discussed in con- 
nection with the maintenance of the unimolecular 
decomposition rates of the paraffins by hydrogen. 
What advantage has the hydrogen molecule over 
other molecules (with the exception of the 
vibrating molecule itself and. its close chemical 
analogs as, for example, its decomposition 
products) in effecting the conversion of trans- 
lational energy into vibrational energy? We 
cannot attempt to answer this question here. A 
very brief review of several methods of attack 
will, however, be given in order to emphasize the 
necessity for further experiment. 

It does not appear that the small mass of the 
hydrogen molecule can be of primary importance 
in this connection. If the laws of molecular 
impact are formulated in a purely classical 
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manner, after the method initiated by Olden- 
berg,” it is possible to imagine that a single 
hydrogen atom could give a part of its trans- 
lational energy to a single more massive atom of 
some other molecule, and that all of the energy so 
imparted should appear as vibrational energy 
after the termination of the collision. The 
necessity for the conservation of linear momentum 
allows hydrogen a maximum advantage of a 
factor of less than two over helium in such a 
process. The experimental data, which indicate 
that this factor is in excess of 20, clearly demand 
another explanation. 

The quantum mechanical approach to the 
problem leads to more promising conclusions. 
Zener™ and Rice!* have attempted solution by 
these means. Franck and Eucken" have reduced 
this treatment to a qualitative picture which is 
more suitable for discussion here. It is supposed 
that the forces which ultimately make themselves 
felt by atomic exchange reactions are, with low 
relative translational energies, capable of causing 
small displacements of the nuclei in the colliding 
molecules. The efficiencies of various collisions 
are thus in a sense conditioned by the chemical 
reactivity of the participating molecules. In this 
respect hydrogen differs profoundly from helium, 
and a qualitative explanation for its advantage is 
found. A somewhat similar variant of the 
quantum mechanical treatment was suggested in 
Part II'® of this series when a mutual ‘‘bond- 
softening’ of the colliding molecules was dis- 
cussed. The great disadvantage of these methods 
of treatment is that they cannot, at present, 
supply numbers which may be compared with 
those obtained from actual measurement. The 
necessity for estimating the duration of the 
collision by the formulation of some dynamical 
mechanism of approach and separation of the 
molecules is a source of difficulty. It is not 


2 Oldenberg, Phys. Rev. 37, 194 (1931); Heil, Zeits. f. 
Physik 14, 31 (1932). 

13 Zener, Phys. Rev. 38, 277, 556 (1931). 

14 Rice, Phys. Rev. 38, 1943 (1931); J. Am. Chem. Soc. 
54, 4558 (1932). 

% Franck and Eucken, Zeits. f. physik. Chemie B20, 
460 (1933). 

16 Richards, J. Chem. Phys. 1, 874 (1933). 


at 
la 
eS 
yn. J 
is 
ire 
in 
nd 
ire 
ns 

A 
ch 
(4) 
to 
A 
ns 
he 

’ 
he 
4.3 
in 
hat 
gen 
ions 
ren- 
ing 

an 


214 ma 


possible, therefore, to accept this type of 
explanation without reserve in spite of its great 
plausibility. 

Yet a third distinctive property of the hydro- 
gen molecule, which has not been greatly 
stressed in this connection, arises from the 
smallness of its moment of inertia. As a result the 
rotational quanta for hydrogen are immense, and 
it has failed to attain a classical distribution of 
rotational energy even at room temperature. 
The occurrence of a dispersive region due to the 
failure of the rotational transitions to follow the 
acoustical cycle is surely due to this property; 
Zener has, indeed, essentially predicted such 
behavior." It is suggestive that, since ethylene is 
conspicuously more successful than hydrogen in 
giving rotational energy to hydrogen, it must also 
be that hydrogen can give its rotational energy to 
ethylene with relative ease. It appears at least 
possible that the magnitude of the rotational 
quanta which are so exchanged may be connected 
with the conspicuous effectiveness of the hydro- 
gen molecule in exciting vibrational energy. If 
approximate ‘“‘matching” or resonance occurred 
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between a hydrocarbon vibrational quantum and 
a hydrogen rotational quantum a fresh basis on 
which to build an explanation would appear. 
Owing to the ill-defined condition of the 
present theoretical background it appears wiser 
to await an expansion of the very meager 
experimental material before attempting an 
explanation of the phenomena which have here 
been discussed. Comparative measurements of 
light and of heavy hydrogen at once suggest 
themselves as likely to provide desirable informa- 
tion. The chemical characteristics of the two are 
practically identical, and hence their behavior 
from the Franck-Eucken standpoint should be 
similar. The mass of the molecule of the heavy 
isotope is very near to that of helium, and any 
purely mass effects should, therefore, be readily 
separable. The larger moment of inertia of heavy 
hydrogen should on the other hand, profoundly 
alter any distinctive behavior which is due to 
rotational energy, and its smaller vibrational 
quantum should enhance any effects due to 
vibrational resonance. In this direction we hope 
for immediate progress. 


Al 
are 
cul 
opt 
thi: 
I 
cul 
are 
ion: 
obs 
tet! 
tha 
| low 
hav 
the 
solt 
| the 
and 
the 
404 
the 
of re 
tot 
stat 
v ob 
| for 
Silic 
| Pho 
Chi 
HC 
Eth: 
Ace: 
gase 
Whe 
solv 
of th 


APRIL, 


This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents, Contributions to 
this section must reach the office of the Managing Editor 


Raman Spectra of HCl 


In the course of study of the nature of dissolved mole- 
cules we have investigated the Raman effect of HCI dis- 
solved in various solvents. The results here communicated 
are for the vibration frequency of HCl in some non- 
ionizing solvents. Woodward! has reported his inability to 
observe the Raman spectrum of HCI dissolved in carbon 
tetrachloride, benzene and xylene, pointing out, however, 
that the negative result was probably to be ascribed to the 
low solubility of HCI in these solvents. This difficulty we 
have circumvented by the obvious device of examining 
the scattered spectra at a temperature at which the 
solubility of HCl in the various solvents was adequate, 
the mol. fraction of HCI in the solutions being about 0.2 
and the temperature —65°C. By »v in Table I is indicated 
the frequency difference between the exciting line (Hg 
4047) and the long wave side of the HCI band, which, if 
the latter is an unresolved Q branch with the same relation 
of rotational lines as in the gas, corresponds approximately 
to the vibrational frequency of the molecule in the lowest 
state. Av describes the difference between » in solution and 
v observed by Wood and Dieke and by Salant and Sandow 
for the gas (2886). 


TABLE I. 


Electric 
moment Diel. Const. v 
x 10'8 at —65° Ap 


0 2860 26 
0.90 10.3 2824 62 
1.03 7 2826 60 


1.08 8.8(—90°C) 2800 86 
Sulphur dioxide 1.63 25 2798 88 
Ethyl bromide 1.83 14 2797 +89 
Acetyl chloride re | 20.5 2804 82 


Solvent 


Silicon tetrachloride 

Phosphorus trichloride 

Chloroform 

HCl liquid (Salant and 
Sandow)? 


The electric moments listed are for the solvents in the 
gaseous state or in dilute solution in nonpolar solvents. 
When » is plotted against the electric moment of the 
solvents, a straight line is obtained, within the accuracy 
of the measurement of v, +5 cm™, up to the value of « for 
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not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


in Non-Ionizing Solvents 


ethyl bromide. It is not at all surprising that the shift in 
frequency should cease to be proportional to the electric 
moment of solvent molecules of high electric moment as 
measured in the attenuated state, for it is well known that 
the polarization of complicated molecules of large electric 
moment may diminish very markedly with concentration. 
The approximate proportionality of the frequency shift 
with the electric moment of the molecules of the solvent 
would seem to indicate that a large part of the shift is due 
to dipole interactions between solvent and solute mole- 
cules, and as the frequency of liquid HCI is of the same 
order as in solvents of the same electric moment as HCl, 
the shift between gas and liquid in this case can also be 
largely attributed to dipole action, as indeed is suggested 
by the fact that the frequency difference between gas and 
liquids for the nonpolar mols He, Ov, Ne, CH,y is much 
smaller than for polar molecules. 

The results have some bearing on the nature of the 
ionization process in HCI molecules. The diminution in 
vibration frequency with increasing dipole moment of the 
solvent indicates a diminution of the force constant for 
the vibration, i.e., a diminution of the strength of binding 
between the hydrogen and chlorine nuclei. The HCl 
molecule in a polar environment may be regarded in a 
state of incipient dissociation, and the diminution of 
binding strength is in a sense the effect recognized by 
Nernst and J. J. Thomson in their correlation between 
ionizing power and the dielectric constant of solvents; but 
it is also very obvious from our results that electrolytic 
dissociation by pure dipole interaction would require a 
solvent of enormously greater electric moment than any 
yet found. Where electrolytic dissociation of HCl actually 
occurs in the solvents in what it has been observed, many 
of which, like water, alcohol, ether, have electric moments 
lower than some of these investigated by us, the essential 


process is probably, as has often been stated, the attach- 


1 Woodward, Phys. Zeits. 32, 777 (1931). 

2 R. W. Wood and G. H. Dieke, Phys. Rev. 35, 1355 
(1930); E. O. Salant and A. Sandow, Phys. Rev. 37, 373 
(1931). 
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ment of a H* by valence forces to the oxygen atom of the 
solvent molecules. 

It is also of some interest to inquire if the increase of the 
electric moment of HCI from a value of 1.081078 e.s.u. 
in the gas to about 1.3X10-'§ in benzene, hexane and 
carbon tetrachloride, reported by F. Fairbrother* could 
be accounted for by the stretching of the molecule in 
solution. 

Assuming the constancy of the expression vr? where r 
is the internuclear distance, one finds that the greatest 
change in », in highly polar solvents, corresponds to an 
increase in r of some 3 percent and in nonpolar solvents it 
would be much less. An increase of 20 percent, to 30 
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percent in the electric moment of HCI dissolved in non- 
polar solvents cannot be attributed to the elongation of 
the molecule which takes place on solution. 

The investigation is being extended to the other hydro- 
gen halides both by the study of Raman and infrared 
spectra. 

Wa. WEst 
PauL ARTHUR 
Department of Chemistry, 
Washington Square College, 
New York University, 
February 20, 1934. 


3F, Fairbrother, J. Chem. Soc., p. 43 (1932). 


On the Calculation of the Probability of the Spontaneous Decomposition of Linear Triatomic Molecules 


In his paper Lifetimes of Unstable Molecules,! N. Rosen 
has calculated the probability of spontaneous decom- 
position of the ‘‘quasi-molecule’’ HO: into the molecule 
O, and atom H. 

His method of calculation is as follows: When we treat 
the molecule HO, as linear, the wave equation of the 
system may be represented in the following form: 


{_¥ 2 @ 


87? mor? Or? Me OY, Ore 


X(n, (1) 


where 7 =X2—%1, 72=X3—X2 (X1, X2 and x3 are the dis- 
tances of atoms m, m2 and m; from a certain fixed point), 
=m and Vi(r;), Vo(r2) 
are the potential energies of the two end-atoms with the 
middle one (the interaction between the end-atoms is 
neglected). 

The term, 


T = dr, = pips (2) 


is considered by Rosen as the perturbation energy; in this 
case the variables in the Eq. (1) can be separated and the 
probability of spontaneous decomposition of the molecule 
is defined in the usual way through the matrix elements of 
(2). 

The perturbation method can be applied, however, only 
in that case, when the matrix elements of the perturbation 
energy are sufficiently small compared with the energy 
differences of the states under consideration, i.e., when: 


(3) 


For the sake of simplicity we approximate the wave 
functions of the unperturbed problem referring to the 
discrete states by the wave functions of two harmonic 
oscillators (which is certainly valid for not too large values 
of the quantum numbers, and m2). The nonvanishing 


matrix elements of the perturbation energy are in this 
case those corresponding to transitions m;—>,+1, Non 
+1. We shall restrict ourselves to the case of two opposite 
jumps corresponding to the smallest values of the difference 
W,—W,+. We get in this case: 


ng, ng! =(1/my2) pri, ne! 
=hA[n,(n2+1)}, if not =n2+1 
=hA[n2(m, +1) }, if m4=n,+1, 


where A= vy, and v2 being the fre- 
quencies of the respective bonds. The condition (3) assumes 
the following form: 


Afni +1) PK | 
(3’) 
A[(m+1)m}PK | v1 


Now it can easily be seen that this condition is not fulfilled 
either for the ‘‘quasi-moleculé’’ HO. or for linear tri- 
atomic molecules known. 

For the ‘quasi-molecule’” HO. we have A =0.6X 
|vi--v2| =6X10" for the values of n(m,=4, 3, 2, 1; 
n2=16, 17, 19, 21) used by Rosen both sides of the in- 
equality (3’) are of the same order of magnitude. For the 
molecule we have: A >|v1—v2|, (A =1.63; |x — 
=1.36). For the symmetrical molecules CO. and 
CS, |v, —v2| =0, the term A[n(n+1)]} is, however, >0. 

It follows that by dropping in Eq. (1) the ‘‘perturbation 
energy” (2) we do not obtain any approximations to the 
actual stationary states of the molecules, so that Rosen's 
perturbation method cannot be applied. 

T. Kontrorova 
V. SOROKIN 
Leningrad, 
Institute of Chemical Physics, 
February 8, 1934. 


1N. Rosen, J. Chem. Phys. 1, 319 (1933). 
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